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INTRODUCTION  AND  SUMMARY 


This  report  gives  a quantitative  description  of  the  uplinks, 
crosslinks,  and  downlinks  of  a general  high  altitude  satellite  sys- 
tem. The  analysis  was  performed  at  intervals  during  1973  and  1974, 
in  support  of  the  Electronic  Systems  Division  (ESD)  responsibility 
to  define  a high  altitude,  crosslinked  system  of  satellites  for 
the  Air  Force  Satellite  Communications  (AFSATCOM)  II  program. 
Specific  quantitative  results  provided  as  a result  of  this  effort 
have  been  incorporated  in  previously  published  reports.  The 
computer  programs  have  also  been  employed  for  determining  test  sites 
having  low  elevation  angles  in  AFSATCOM  I testing.  The  present 
report  documents  the  detailed  analysis  and  associated  computer 
programs,  and  provides  further  examples. 

Efforts  to  define  an  AFSATCOM  II  system  have  required  analysis 
of  medium  altitude  (3,000  to  20,000  nmi) , high  altitude  (20,000  to 
250,000  nmi),  and  highly  eccentric  (eccentricity  >0. 7)  orbits. 
Various  system  concepts  also  employed  satellite  to  satellite  links 
or  crosslinks  as  shown  in  Figure  1.  While  there  were  some  computer 
programs  available  for  problems  in  each  class  of  orbit,  there  were 
severe  limitations  on  their  utility  in  context  of  the  AFSATCOM  sys- 
tem problem.  Close  analysis  of  such  a general  satellite  communica- 
tion system  (uplinks,  crosslinks,  and  downlinks)  has  been  going  on 
only  for  the  past  few  years. 

Principal,  significant  results  of  the  work  described  in  this 
report  are: 
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• Analytic  solution  for  high  altitude  lunar  perturbations, 
which  has  allowed  study  of  a wide  range  of  orbit  para- 
meters because  of  the  resulting  programs  are  less  costly 
to  run  than  previously  used  numerical  integration  pro- 
grams . 

This  led  to  a conclusion  in  a particular  case  that  12-day 
retrograde  orbits  have  stability  comparable  to  10-day  posigrade 
orbits. 

• Computer  programs,  which  offer  the  following  extensions  in 
comparison  to  other  programs  known  to  us: 

1.  Unified  uplink,  crosslink,  downlink  analysis  in  a 
single  program. 

2.  Programs  which  handle  high  eccentricity  orbits 
efficiently. 

3.  Allowance  for  incorporation  of  vector  antenna  radia- 
tion patterns  for  precise  received  power  at  the  ground 
terminals. 

4.  Low  central  processing  unit  (CPU)  time  requirements 
to  run  the  programs.  For  example,  less  than  30 
seconds  CPU  time  (IBM  370/158)  is  required  for  uplinks, 
crosslinks,  and  downlinks  for  two  ground  stations  and 
10  satellites. 

• Coverage  of  multiple  satellites. 


/• 


> 

, 


Part  I of  the  report  includes  a simple  two-body  analysis  (satel- 
lite, earth),  which  is  a good  basis  for  the  communications  engineer 
in  designing  orbital  links  of  3000  to  20,000  nautical  miles  (nrai) 
altitude.  General  Keplerian  orbits  are  analyzed  with  the  aid  of  a 
notably  efficient  iterative  solution  to  Kepler's  equation. 

Part  II  considers  orbital  altitudes  between  20,000  to  250,000 
nmi.  Analysis  of  higher  eccentricity  orbits  was  achieved  through 
use  of  a Taylor  series  expansion  of  eccentric  anomaly,  which 
allows  both  analysis  of  a higher  eccentricity  and  lower  CPU  time 
than  standard  methods. 

Analysis  of  high  altitude  orbits  had  been  previously  carried  out 
using  a numerical  integration  program.  ^ For  altitudes  greater 
than  approximately  20,000  nmi,  lunar  and  solar  perturbations  signi- 
ficantly disturb  the  satellite  orbit.*  Results  had  shown  that  retro- 
grade equatorial  satellite  orbits  are  more  stable  than  posigrade  orbits. 
A new  approach  was  adopted  of  combining  perturbed  orbital  elements 
to  provide  an  analytic  solution  to  orbit  perturbation  and  resulting 
stability. 

Derivation  of  the  rate  of  change  of  semi-major  axis  has  provided 
physical  insight  into  the  mechanism  which  causes  satellites  moving 
retrograde  with  respect  to  the  moon's  motion  to  be  more  stable  than 
those  moving  posigrade.  This  result  had  been  inexplicable  previously. 

As  noted  above,  the  low  running  cost  of  the  program  has  allowed  more 
extensive  analysis  of  orbits  for  stability.  As  in  Part  I,  the  uplink, 
crosslink,  and  downlink  analysis  is  reduced  to  a program. 
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Examples  of  the  use  of  Part  I and  II  programs  are  given  in  Part 
III.  Computer  programs  are  listed  and  described  in  Appendices  1 
through  9.  The  table  of  contents  identifies  the  author  of  each 
program.  A portion  of  the  above  programs  was  used  to  generate  a 
succinct  program  for  coverage  in  a multiple  satellite  system  (Example 
3 in  Part  III  and  Appendix  9). 


The  programs  are  intended  for  actual  communications  link  analy- 
sis. The  positions  of  all  stations  and  satellites  are  retained  in 
vector  form  in  the  programs  so  that  antenna  gain  patterns  can  be 
added  later  for  specific  link  calculations.  It  is  expected  that  these 
programs  will  be  useful  in  high  altitude  satellite  communications 
planning.  In  addition,  the  low  CPU  time  of  the  attached  programs 
may  allow  a spacecraft  with  limited  processing  capability  to  autonomously 
calculate  its  own  available  communications  links. 


PART  I 


MEDIUM  ALTITUDE  ORBITS 


Orbital  altitudes  between  approximately  3000  to  20,000  nmi 
(medium  altitude)  allow  the  possibility  of  an  especially  simple 
orbit  analysis.  A non-rotating  coordinate  system  is  chosen  a6  a 
basis  for  the  analysis  to  avoid  any  acceleration  problems  inherent 
in  the  coordinate  system  which  would  arise  if  Doppler  rates 
were  added  to  the  analysis  later.  The  coordinate  system  is  shown 
in  Figure  2;  it  is  a Cartesian  system  with  the  x axis  pointing  to 
the  vernal  equinox  and  the  origin  at  the  geocenter.  The  x - y 
plane  is  the  equatorial  plane.  This  coordinate  system  translates 
with  the  earth  as  it  moves  around  the  sun,  but  does  not  rotate. 

Aries  (the  direction  of  the  x-axis  and  the  vernal  equinox)  is  so 
far  away  that  movements  of  the  earth  around  the  sun  cause  insignifi- 
cant changes  in  direction  of  the  axes. 

In  Figure  2,  note  that  all  four  sites  are  moving.  The  ground 
stations  are  rotating  with  the  earth  at  15°/hr,  and  the  motion  of 
the  satellites  is  determined  by  Kepler's  laws. 

A convenient  starting  point  for  an  orbital  analysis  (time  = T - 
0 hrs)  can  occur  when  earth  coordinates  equal  celestial  coordinates. 
This  occurs  at  12  noon  Greenwich  mean  time  (GMT)  on  March  21.  Imme- 
diately after,  the  inertial  longitude  of  an  earth  station  will  be 
greater  than  the  earth  coordinates. 


ePRECfDirG  rlGS  BLANK. NOT  FILMED 
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ONE  WAY  TO  LINK  GROUND  STATIONS  I AND  2 IS  SHOWN.  THE 
COORDINATE  SYSTEM  IS  STATIONARY,  BUT  ALL  FOUR  SITES  ARE 
MOVING.  THE  ORBIT  OF  SATELLITE  SI  CAN  BE  COMPLETELY 
DESCRIBED  FOR  KEPLERIAN  ORBITS  BY 
W,  = RIGHT  ASCENSION,  DEGREES 

i,  = INCLINATION  WITH  RESPECT  TO  THE  EQUATORIAL  PLANE,  DEGREES 
Wp<  = ARGUMENT  OF  PERIGEE,  DEGREES 
TP|  = TIME  OF  PERIGEE,  HOURS 

a = SEMIMAJOR  AXIS,  NAUTICAL  MILES  (OR , IN  KM  BY  MULTIPLYING  THE 
SEMIMAJOR  AXIS  BY  1.852) 

• * ECCENTRICITY,  WHERE  0 < e <1.  AT  e ;0,  THE  ORBIT  IS  CIRCULAR 


Figure  2 GEOMETRY  OF  A TWO  SATELLITE,  TWO  GROUND  STATION  SYSTEM 


1.1  MOTION  OF  THE  SATELLITE 


R = P cos  0 
1+e 


(l-D 


The  parameters  are  shown  in  Figure  3. 

2.  In  the  coordinate  system  of  Figure  4 (called  the  prime 
system  in  the  remainder  of  the  discussion) , the  radius  vector  of 
each  satellite  sweeps  through  equal  areas  in  equal  times: 


R 0 » constant 


where  [l  * Gm 


(1-2) 


earth's  gravitational 
constant. 


3.  The  squares  of  the  periods  of  the  satellites  are  to  each 
other  as  the  cubes  of  the  semi-major  axes  of  their  respective  orbits: 


Tl2 

— 3 = — 3 
al  a2 


(1-3) 


After  extensive  analysis,  these  three  laws  can  be  interpreted  as 
a single  equation  (Kepler's  equation): 


Only  elliptical  orbits  are  considered  here,  because  a circular 
orbit  can  be  considered  a degenerate  ellipse.  Kepler's  laws  of 
planetary  motion  can  be  abbreviated  as : ^ 

1.  The  orbit  of  a satellite  is  an  ellipse,  with  the  earth's 
center  (geocenter)  at  one  of  the  foci: 
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THE  COORDINATE  SYSTEM  IS  AT  ONE  FOCUS  (GEOCENTER)  OF  THE  ELLIPTICAL 
ORBIT 

lure  3 RELATION  OF  TRUE  ANOMALY  6 TO  ECCENTRIC  ANOMALY  E 


M =■  E - e sin  E 


(1-4) 


where  M ■ mean  eccentric  anomaly  ■ n(t-tp) 
n = mean  angular  rate 
t ■ time,  hrs 
tp  3 time  at  perigee 
e = eccentricity 
and  E - eccentric  anomaly. 


The  relation  of  eccentric  anomaly  E to  true  anomaly  0 can  be 
seen  in  Figure  5.  The  argument  of  perigee  is  conveniently  located 
here  on  the  apsidal  axis.  The  true  anomaly  is  shown  as  the  angle 
from  perigee,  measured  from  the  focus  at  the  earth's  center.  Serious 
problems  arise  when  0 (t)  is  desired,  even  for  these  simple 
Keplerlan  orbits.  Many  extended  analyses  have  been  attempted  to  de- 
scribe the  progression  of  the  satellite  through  its  orbit  as  a func- 

[3] 

tion  of  time.  Moulton  developed  one  short  result  of  an  impor- 
tant analysis  which  is  useful  for  eccentricity  < 0.5.  Eccentricities 
greater  than  0.5  were  desired  for  this  analysis  and  for  most  of 
these  programs,  however. 


A direct  attack  on  Kepler's  equation  (Equation  (1-4)  has  been 
employed.  Although  the  Kepler  equation  is  transcendental,  an  ini- 
tial guess  at  eccentric  anomaly 


E^  = M + e sin  M 


(1-5) 


enables  one  to  take  a second,  much  more  accurate  estimate. 


M + e (sin  E^  - (e  cos  E^) 
1 - e cos  (E^) 
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IA  - 45,623 


<*'  = w + e 
T« 

Figure  5 MOTION  IN  THE  ORBITAL  PLANE 
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This  estimate  for  E2  is  the  result  of  linearizing  Kepler's 
equation  by  means  of  a Taylor  series  expansion.  All  terms  beyond 
the  first  power  in  E^  are  neglected.  More  estimates  for  eccen- 
tric anomaly  can  be  similarly  made  (e.g. , E3  as  a function  of  Ep. 
When  differences  of  successive  approximations  are  adequately  small, 

E,  can  be  accepted  as  the  solution  for  E.  The  true  anomaly  9 is 
J [2] 

related  to  the  eccentric  anomaly  by 


e 


cos 


-1 


/ cos  E - e 
( 1 - e cos  E 


(1-7) 


This  true  anomaly  will  be  needed  to  accurately  determine  the 
time-varying  position  of  the  satellite  on  its  prescribed  ellipse. 
The  range  from  the  geocenter  is  given  as 


R „ a (1  - a2) 

1 + e cos  9 


(1-8) 


These  relations  (Equations  (1-5)  to  (1-8))  are  implemented  in 
subroutine  ELLIP  of  the  computer  programs.  The  iterations  for  E are 
surprisingly  fast;  representatively,  they  are  as  fast  as  Moulton's 
esoteric  development,  which  was  meant  expressly  to  free  the  early 
twentieth  century  analyst  from  laborious  computations.  The  advan- 
tage over  Moulton's  development  is  that  the  programs  can  handle 
eccentricity  from  0 to  0.99.  The  subroutine  ELLIP  can  be  forced 

into  seven  iterations  at  e ■ 0.99.  This  iterative  technique  requires 

[4] 

less  CPU  time  than  a standard  procedure.  In  the  subroutine, 

M becomes 


2 7T,(t-t£j  _ 


z * 


r 


mean  angular  rate,  radians 


I 


Also,  <p  ’ “ Wp  + 9 since  9 is  measured  from  the  argument  of  perigee, 

s 

but  <f>  ' Is  measured  from  the  line  of  nodes  where  the  orbital  plane  In- 
s 

tersects  the  equatorial  plane. 

1.2  COORDINATE  TRANSFORMATIONS 

In  the  previous  section,  It  was  found  convenient  to  describe 
angular  motion  (true  anomaly  9)  in  the  plane  of  the  satellite.  How- 
ever, after  a position  (or  a velocity)  in  the  satellite  plane  has 
been  found,  it  should  be  converted  back  into  an  inertial  coordinate 
system,  which  is  required  for  correct  interstation  vectors. 

Fortunately,  a coordinate  transformation  that  compares  directly 
to  the  needs  here  has  been  extensively  analyzed.  ^ This  trans- 
formation makes  use  of  the  Euler  angles.  The  first  Euler  angle 
occurs  with  a rotation  about  the  z axis.  See  Figure  6,  where  the 
nomenclature  of  Goldstein  is  used.  This  first  rotation  through  W 
will  correspond  to  right  ascension  for  this  analysis. 

The  second  rotation  occurs  about  £ ; again,  this  rotation  is  coun- 
terclockwise. This  rotation  is  an  angle  i,  where  i here  corresponds 
to  orbital  inclination.  Finally,  a rotation  ($'  ) in  the  orbital 
plane  can  occur  which  is  measured  from  the  line  of  nodes.  See  Figure  4. 

The  three  rotations  define  a new  (PRIME)  coordinate  system 
(x'.y'jz').  For  the  purpose  of  this  analysis,  the  x'  axis  is  a 
radius  vector  from  the  geocenter  to  the  satellite. 

The  transformation  of  a vector  x'  in  the  prime  coordinate  sys- 
tem to  the  inertial  system  can  be  expressed  in  matrix  form  as 

x'  * A x 


■ 


: 
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45,621 


man 


and  A is  the  product  of  three  separate  matrices,  15,  C,  and  D.  Omit 
Goldstein’s  intermediate  steps, 

cos  <t>'  cos  W - cos  i sin  W sin  d>'  , cos  p'  sin  W + cos  i cos  W sin  <f> ' , 

T s s s s 

+ sin  W sin  i 

A = -sin  d> * cos  W - cos  i sin  W cos  4'  -sin  <f>'  sin  W + cos  i cos  W cos  d> ' 
— s s s 

cos  d>’  sin  i 
s 

sin  i sin  W,  - sin  i cos  W,  cos  i (1-9) 

and  the  inverse  transformation  from  satellite  coordinates  to  the 
inertial  frame  is  of  more  immediate  interest: 

-1  1 
x = A x 

(A  subroutine  UNPRIM  will  later  use  the  matrix  A directly  for  a 
local  satellite  coordinate  system.  UNPRIM  will  be  used  specifically 
for  crosslink  pointing  angles  in  program  SATLUNAE.) 

A ^ is  found  from  the  transpose  of  A because  the  magnitude  of 
A is  unity.  The  elements  of  A ^ are  then 

a,,  = cos  ' cos  W - cos  i sin  W sin  * 

11  s s 

a,  _ = -sin  $ ' cos  W - cos  i sin  W cos  p' 

12  s s 

a^  = sin  i sin  w 

Si  = cos  $ 1 sin  W + cos  i cos  W sin  <p  * 

21s  3 
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a0,  = -sin  <p  ' sin  W + cos  i cos  W cos 
t-t  s 

■ -sin  i cos  W 

a..  - sin  i sin  d> ' 

31  s 

a = sin  1 cos  d> ' 
o l s 

a33  =»  cos  i 

The  subscripts  refer  to  row  and  column. 


(1-10) 


All  of  these  elements  except  a13,  a23  and  a33  are  inserted  into 
the  attached  subroutine  PRIME.  They  are  omitted  because  there  is  no 
z component  for  satellite  coordinate. 

An  example  of  the  calling  of  subroutine  PRIME  can  be  useful.  If 
the  position  of  Satellite  //8  in  its  own  coordinate  system  is  spe- 
cified by  x'  = R8,  y'  =0.,  z'  -0.,  this  information  must  be  sent  to 
PRIME  along  with  the  three  angles  specifying  the  coordinate  transforma- 
tion. The  angle  of  the  satellite  from  the  line  of  nodes  is  FSP,  the 
right  ascension  is  W8,  and  the  inclination  is  xl8  for  the  purpose  of 
this  example.  The  attached  programs  call  PRIME  in  the  following  way. 

CALL  PRIME  (FSP,  W8,  xl8,  R8,  0.,  xs,  ys,  zs) . 

input  variables  output  variables 

The  output  variables  xs,  ys,  and  zs  represent  position  in  inertial 
space.  The  position  will  be  useful  in  getting  crosslink  and  downlink 
vectors. 
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A slightly  different  form  of  coordinate  conversion  will  be  re- 
quired for  range  rate.  When  it  is  discussed  later,  it  should  not  be 
confused  with  the  prime  coordinate  system  in  which  the  x'  axis  is 
pointed  through  the  satellite.  It  will  use  a stationary  x"  axis 
pointed  along  the  line  of  nodes. 

1.3  DOPPLER  ANALYSIS 

In  a central  force  field  such  as  that  assumed  for  Keplerian 
orbits,  the  velocity  of  a satellite  is  dependent  only  on  the  magnitude 
of  its  radius  vector  once  its  initial  conditions  have  been  determined. 
Similarly,  it  is  also  determined  by  9 (its  angle  measured  from  perigee) 

since  R(e)  is  determined  by  the  Keplerian  ellipse.  Expressions  directly 

[2] 

related  to  velocity  are 


i ; 


[■ 


R = 


(2ttA  e ) 
t v/l-e2 


Sin  9 


and 


0 = 


2 tt  (1  + e cos  9)  ‘ 


t (1-e  ) 


2\  3/2 


(1-11) 


(1-12) 


The  geometry  is  given  by  Figure  6, 

where  $ ' = 9 + W , W = argument  of  perigee 

s P P 

and  ($')  = 9 . 

The  x"  axis  is  stationary  and  lies  in  the  equatorial  plane. 


:: 
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Further,  the  derivatives  of  the  position  vector  (with  the  x"  axis 
now  lying  along  the  line  of  nodes  and  the  y"  axis  90°  ahead  in  the 
orbit  plane)  can  give  the  x"  component  of  velocity  and  the  y"  compon- 
ent of  velocity.  The  symbol  0 ' is  later  changed  to  FSP  so  it  can  be 

s 

more  directly  related  to  the  FORTRAN  programs. 

x"  component  = ^qq)  (R  cos  (0'g)  - R 9 sin  (0’s))  nmi/sec  (1-13) 

and 

1 

y"  component  = (j^Qg) (R  sin  Wfl)  + R 9 cos  (0'g))  nmi/sec  (1-14) 

To  convert  these  velocity  components  to  components  in  the  inertial 
system,  the  general  coodinate  conversion  (PRIME)  can  be  called,  but 
with  0 ' =0,  because  the  x"  axis  does  not  rotate  with  the  satellite 

in  this  case.  Velocity  components  are  found  in  the  subroutine  DOPE  of 
the  attached  programs.  This  stationary  x"  axis  was  not  required  for 
range  rate  but  will  be  desirable  if  R is  analyzed. 

The  velocities  of  the  earth  stations  are  straightforward. 

*1  ~3^  Rfe  Sin  91  8in  *1 
h = 3^0  Re  sin  91  cos  *1 


where 


15°/hr 

57.296 


earth  rotation  rate, 


rad/hr  * 


*Note  that  the  distinction  between  solar  hours  and  sidereal  hours  is 
not  retained  here  because  only  three  place  accuracy  is  desired  in  the 
answers. 
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z^  = 0 because  rotation  is  about  the  z axis 


The  background  for  velocitites  of  satellites  and  earth  stations 
is  now  concluded.  The  relation  of  velocities  to  uplink  and  crosslink 
Doppler  shift  is  of  more  interest  to  the  communications  engineer  than 
velocities  per  se. 

The  relation  of  a Doppler  shifted  frequency  to  the  transmitted 
frequency  is  shown  by  Jackson ^ to  be 


(1  - — cos  0j) 


(l-15a) 


Where  v is  relative  velocity  between  transmitter  and  re- 

K 

ceiver  and  0^  is  the  angle  between  the  relative  velocity  vector  and 
the  pointing  vector  between  transmitter  and  receiver. 

This  expression  includes  the  (relativistic)  transverse  Doppler 
shift  which  may  be  of  interest  for  crosslink  Doppler.  For  vR  « c. 
Equation  (l-15a)  becomes 


v 2 v 

F (1  + 1/2  -y-)  (1  - y-  cos  0j) 


(l-15b) 


A F = Doppler  shift  = F (1/2  ) 

2 


c (vR  C0S  V 


(l-15c) 


The  two  components  are  called  transverse  Doppler  shift  and  ordi- 
nary Doppler  shift,  respectively.  At  synchronous  altitude,  the  trans- 
verse Doppler  shift  can  be  of  the  order  of  (first  term  of  Equation 
(l-15c)) 

(35  x 109  Hz)  (I/2A-*65-5  nml/8fc \ ) - i.82  Hz  at  K band. 

\l.6198  x 10  nmi/sec J 

This  is  so  small,  even  with  a millimeter  wave  carrier,  that 
tranverse  Doppler  is  dropped  for  the  remainder  of  the  discussion 
and  in  the  computer  programs  (Appendices  1-9). 


The  ordinary  Doppler  shift  is 


F - - - (VR  cos  9j) 


In  order  to  get  (cos  6^),  the  relative  velocity  vector  v 
is  found  by  taking  differences  of  velocity  components  in  the  inertial 
f rame : 


'r  “ ^*2  ” *i>  i + ($2  “ J + ^*2  " *i^  k 


By  \ or  analysis, 
vr  (cos  Oj)  “ vR  . p~ 


where 


Pp  ■ unit  vector  which  points  from  transmitter  to  receiver. 
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(x2  - x^  i + (y2  - yx)  j + (z2  " z{>  k 

Rl 


VR  (cos  0j)  - vR  . P12 


(x2-x1)  (x2-x1)  + (y2~y j)  ^-Vj)  + (z2_zi)  (‘7"fzl> 


(1-16) 

and  the  one-way  Doppler  shift  can  be  found  by  substituting  this  result. 
Equation  (1-16) , into 


AF  = — (vR  cos  Gj). 

This  calculation  is  done  in  the  attached  programs  after  the  sub- 
routine DOPE  is  called. 


1.4  ELEVATION  ANGLE 

The  elevation  angle  for  the  ground  station  antenna  is  important 

to  the  communications  engineer  in  many  ways:  (1)  it  serves  as  a check 

for  satellite  visibility;  (2)  it  can  give  an  estimate  of  mean  atmospheric 
attenuation  of  signal  strength  and  (3)  its  derivative  is  useful  in  an 
estimate  of  antenna  slewing  rate.  Only  the  elevation  angle  is  treated 
here,  and  not  its  derivative. 

Elevation  angle  is  found  by  noticing  the  following  geometry,  which 
is  in  the  plane  of  the  ground  station  position  vector  and  the  satellite 
position  vector  (Figure  7). 


The  geocentric  angle  A. , can  be  found  by  taking  the  dot  product 


R1  * R1SJ 
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43,«13 


r 


A1S  = cos 


-1  / <xlxs  + ylyS  + Zl2s) 


R R 
e s 


with  the  aid  of  the  law  of  sines,  the  elevation  angle 


(1-17) 


D1S  ~ sin  1 / (sin  Alg)  — 

\ R1S 


(1-18) 


is  found.  A sign  check  for  D^g  is  still  necessary  because  the  obtuse 
angle  (D^g  + tt/2)  in  Figure  7 is  read  as  (rr/2  - D^)  by  the  computer. 

It  cannot  distinguish  between  sines  of  the  first  and  second  quadrants. 


Azimuth  angle  calculations  are  also  included  in  the  program  AZ1 
(Appendix  8). 

1.5  CROSSLINK  POINTING  ANGLES 

The  angles  are  found  for  a vector  Rg^  in  a coordinate  system 
centered  at  Satellite  // 8.  The  coordinate  system  has  axes  parallel 
to  the  inertial  coordinate  system.  This  kind  of  coordinate  system 
makes  sense  only  if  the  satellite  has  to  know  the  directions  of 
the  inertial  axes  anyway.  Otherwise,  the  coordinate  system  is  a 
local  satellite  coordinate  system  (shown  as  angles  TU12  and  FU12 
of  the  program  SATLUNAE  in  Appendix  6) . 


The  angles  are  found  in  terms  of  local  longitude  and  colatitude. 


Colatitude  0g^  g£  = cos 


-1  ( 2S2  ~2S1 
\RS1,  S2  . 


(1-19) 


Longitude 


*S1,S2  = t3n 


-1  ( yS2  ~ySl 
\XS2  ”XS1 


(1-20) 
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These  angles  are  found  near  the  end  of  the  main  part  of  the 
programs.  Again,  some  sign  tests  must  occur  because  the  computer 
has  trouble  with  correct  quadrants.  Another  problem,  which  is  merely 
inconvenient,  remains:  Longitudes  less  than  -360°  sometimes  occur. 

If  this  is  too  inconvenient  for  the  user,  a two-line  logical  check 
can  be  done  as  on  lines  662  and  664  of  the  program  AZ1  (Appendix  8). 

Crosslink  visibility  is  also  calculated  in  the  attached  programs. 
It  is  found  by  checking  the  angle  between  2S  and  (-R^g)  * ^ 
this  angle  is  greater  than  half  the  subtended  earth  angle,  the  cross- 
link is  declared  visible  (CVIS  * 1 in  the  programs). 

1.6  MEDIUM  ALTITUDE  PROGRAMS 

Four  FORTRAN  IV  computer  programs,  all  based  upon  the  flowchart 
listed  on  Blocks  1 and  2 (Figures  8a  and  8b)  have  been  written  to 
examine  various  relations  regarding  satellite  links  for  a maximum  of 
ten  satellites  and  two  ground  stations.  The  four  programs  may  be  di- 
vided into  two  categories  according  to  the  manner  in  which  input  data 
is  supplied  to  the  program:  those  for  which  input  data  are  supplied 

by  answering  input  prompting  questions  during  program  execution  at 
a Time  Sharing  Option  (TSO)  terminal  and  those  for  which  input  data  are 
supplied  by  changing  the  data  statements  within  the  program  before  exe- 
cution. 

All  four  programs  require  the  following  input  data: 

1.  The  semi-major  axis  in  nautical  miles  for  each  satellite. 

2.  The  eccentricity  for  each  satellite. 

3.  The  right  ascension  in  degrees  for  each  satellite. 

4.  The  argument  of  perigee  in  degrees  for  each  satellite. 

5.  The  time  of  perigee  in  hours  for  each  satellite. 
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Figure  8 a FLOWCHART,  BLOCK  I 
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\SET  VALUE  OF  Ml  AMD  M2 
\ M2* 


7 


, Ml  AMD  M2  AHE  2 OF  THE  10 
SATELLITES  FOR  WHICH  THE 
CROSSLIN*  VALUES  WILL  SE 
COMPUTED 


TAU  * (S.9S7  ■ 10'*  ) A(Ml)S/2 
TAU2  * (S.GBT  ■ I0_«  >A(M2)V2 

00  )001  *1,25 

J 

I TAU- PERIOD  OF  ROTATION 
OF  SATELLITE  Ml 
TAU2  - PERlOO  Of  ROTATION 
OF  SATELLITE  MS 


i Rt -HOURS 


CALL  ELLIP 


DETERMINE  R, RANGE  FROM  GEOCENTER  TO 
SATELLITE  Ml  S fSP,  ANGLE  MEASURED  IN 
ORBIT  PLANE  [REFERRED  TO  NOOE  OF  ORBIT 
PLANE  ANO  EQUATORIAL  PLANE] 


I 


CALL  ELLIP 

/ DETERMINE  R2,  RANGE  FROM  QEOCENTER  TO' 


SATELLITE  M2  B F2P,  ANGlE  MEASUREO  IN 
ORBIT  PLANE  [REFERRED  TO  NODE  OF  ORBIT 
^ PLANE  AND  EOUATORjAL  PLANE] 


/oF 

CAf 

Lmi 


TERMINE(XS,VS,23),THE  INERTIAL 
CARTESIAN  COORDINATES  OF  SATELLITE 


CALL  »RIME 


DETERMINE  (X2.Y2.22),  THE  INERT lAl 
CARTES. AN  COORDINATES  OF  SATELLITE 


|Rl2»lU2-X5)*+(Y2-YS)*  + (22 -ZS)*)''* 

T 


,ftl?  - RANGE  BETWEEN  SATELLITES 
Ml  B M2 


FI2  ( a T 12 ) - POINTING  ANGLES 
BETWEEN  SATELLITES  Ml  AND 
M2  (THE  LOCAL  COORDINATE  SYSTEM 
CENTERED  AT  Ml , WITH  AXES  PARALLEL 
TO  TME  GEOCENTRIC  INERTIAL  SYSTEM) 


, T)2  F12)  - POINTING  ANGLES 

BETWEEN  SATELLITES  Ml  AND  M2 


T52*  37. 29®  [ARC COS  ((22-  2S)/Rl2)]  | 

x - Z 


f FOR  SATE 

(XV,  f 2 


SATELLITE  Ml,  DETERMINE 

V),  WHERE  ( XV,  YV,2V)*(x's,  YS.2S 


X 


D 


FOR  SATELLITE  M2,  DETERMINE 
(XV2.YV2.2V2),  WHERE  (XV2.YV2. 2V2)  * 
(X2.V2.22) 


T 


[dot*  [(XV2-XV)(X2-XS)-MYV2- YV)(Y2-VS)4(2V2-2V)(22-2S)]/RI2  [ 


; DOPI2  - DOPPLER  SHIFT  EXPERIENCED 
AT  SATELLITE  M2  DUE  TO  SIGNAL 
TRANSMITTED  FROM  SATELLITE  Ml 


OOP  !2-£-OOT»  !)•■  10*)] /i  ,10* 


Figure  8b  FLOWCHART,  BLOCK  2 


33 


M 


r 


6.  The  inclination  in  degrees  for  each  ground  station. 

7.  The  colatitude  in  degrees  for  each  ground  station. 

8.  The  longitude  in  degrees  for  each  ground  station. 

9.  The  uplink  frequency  in  Hz. 

10.  The  downlink  frequency  in  Hz. 

11.  The  crosslink  frequency  in  Hz. 

All  four  programs,  using  the  above-mentioned  input  data,  compute  the 
following: 


1.  The  range  in  nautical  miles  between  each  satellite  and  each 
ground  station. 

2.  The  elevation  angle  in  degrees  from  each  ground  station  to 
each  satellite. 

3.  The  uplink  Doppler  shift  in  Hz  for  each  satellite  and  each 
ground  station. 

4.  The  downlink  Doppler  shift  in  Hz  for  each  satellite  and  each 
ground  station. 

5.  The  9ubsatellite  colatitude  in  degrees  for  each  satellite. 

6.  The  subsatellite  longitude  in  degrees  for  each  satellite. 

7.  The  range  in  nautical  miles  between  two  satellites  for  all 
pairs  of  satellites. 

8.  For  all  pairs  of  satellites,  the  pointing  angles  in  degrees 
between  a first  satellite  and  a second  satellite  (with  the 
local  coordinate  system  centered  at  the  first  satellite  and 
having  axes  parallel  to  the  geocentric  inertial  system). 

9.  For  all  pairs  of  satellites,  the  one-way  crosslink  Doppler 
shift  in  Hz  experienced  at  a second  satellite  due  to  a signal 
transmitted  from  a first  satellite. 

10.  The  crosslink  visibility  between  two  satellites  for  all  pairs 
of  satellites. 
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1.6.1  SATE  (Appendix  1) 


The  first  program  listed  has  been  entered  on  TSO,  and  may  be 
executed  either  in  the  foreground*  with  the  output  printed  at  the 
terminal  or  in  the  background  with  the  output  printed  on  the 
high  speed  printer.  Input  data  is  entered  in  this  program  by 
changing  the  data  statements  within  the  program  before  execution. 

1.6.2  SAID  (Appendix  2) 

The  second  program  listed  is  simply  the  double  precision 
version  of  SATE;  that  is  to  say,  SATD  has  the  same  attributes 
as  SATE  except  for  the  fact  that  it  has  a precision  of  approximately 
16  decimal  digits  as  opposed  to  the  precision  of  approximately 
seven  decimal  digits  of  SATE.  Because  of  this  increased  accuracy, 
SATD  has  a running  time  longer  than  that  of  SATE. 

1.6.3  SATV1ZE  (Appendix  3) 

The  third  program  listed  has  also  been  entered  on  TSO  but,  un- 
like the  two  above-mentioned  programs,  is  on-line  foreground  execu- 
table only  with  output  being  printed  at  the  terminal.  Also,  the 
necessary  input  data  in  this  program  are  entered  during  program  exe- 
cution merely  by  answering  the  input  prompting  questions  supplied  by 
the  program.  For  the  casual  user,  default  values  and  points  for 
stopping  various  parts  of  program  execution  have  been  incorporated 
within  the  program. 


*To  run  this  program  in  the  foreground,  lines  10  to  40  inclusive 
roust  be  deleted  first. 
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1.6.4  SATVIZD  (Appendix  4) 

The  fourth  program  listed  Is  merely  the  double  precision 
version  of  SATVIZE.  Because  of  its  increased  precision,  SATVIZD, 
like  SATD,  has  a running  time  longer  than  that  of  its  single 
precision  counterpart. 


PART  II 


HIGH  ALTITUDE  SATELLITES  WITH  LUNAR  PERTURBATIONS 


Conventional  three-body  (earth-moon-satellite)  perturbation 
analysis  can  consume  great  amounts  of  computer  time  with  concomitant 
expense.  For  example,  the  carefully  planned  and  very  extensive 
Lincoln  Laboratory  Planetary  Ephemeris  Program  (PEP)  can  perform 
numerical  integrations  of  satellite  motion,  including  the  perturba- 
tions of  the  satellite  orbit  by  many  non-terrestrial  bodies.  The 
implicit  penalty  of  this  program  is  its  expense.  The  disadvantages 
of  numerical  integration  become  even  more  marked  for  a multiple 
satellite  communications  system.  For  these  reasons,  and  to  account 
for  long  term  lunar  perturbations,  programs  incorporate  existing 
short  analytic  (i.e.,  closed-form)  results.  Dr.  M.  Ash's  elegant 
and  useful  analytic  results  ^ for  variation  of  right  ascension, 
argument  of  perigee,  and  eccentricity  are  combined  with  a short 
approximate  result  for  the  time  rate  of  change  of  a semi-major 
axis  which  we  developed.  A changing  semi-major  axis  can  account 
for  the  difference  in  stability  of  posigrade  and  retrograde 
orbits. 

The  programs  discussed  and  listed  in  this  section  (Appendices 
5-7)  are  the  only  ones  known  to  us  that  can  be  used  to  estimate 
orbital  stability  with  modest  computation  time,  while  noting  the 
difference  between  posigrade  and  retrograde  orbits. 

2.1  ANALYSIS 

One  way  to  avoid  extensive  computations  for  the  position  of  a 
satellite  in  a non-Keplerian  orbit  is  to  describe  the  fundamental 
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motion  at  any  time  as  Keplerian  with  defined  orbital  elements.  How- 
ever, these  orbital  elements  can  be  allowed  to  change  as  a function 
of  the  perturbing  forces.  When  the  Keplerian  motion  is  elliptical, 
the  corresponding  orbital  elements  for  perturbed  motion  are  called 
"osculating"  elliptic  orbital  elements. 

With  a few  changes  to  the  nomenclature  of  Ash, 

a = semi-major  axis 

e = eccentricity 

I = inclination  with  respect  to  lunar  plane 

£>  = right  ascension  of  ascending  node  on  lunar  plane 
Wp  = argument  of  perigee 

M = mean  anomaly 

^ = gravitational  constant  times  lunar  mass 

^ = gravitational  constant  times  mass  of  earth 
1/2  -3/2 

n = U a = mean  motion,  rad/hr 


= true  anomaly 
= semi-latus  rectum 


= radius  from  geocenter 


1 + e cos  \[> 
lunar  radius 
0 + W 


Further,  if  R,  S,  and  W are  the  three  orthogonal  perturbation  force 
components  (R  along  radius  vector  from  geocenter,  S close  to  the 
velocity  vector,  and  W completing  the  right  handed  coordinate  system) 
for  a lunar  mass  spread  into  a torus  at  lunar  altitude.  Ash  derived 
the  relations  for  R,  S,  W given  by  Equations  (2-1),  (2-2)  and  (2-3). 


j, 


K>|3 


m 

2 


2*  (f)2*'1  yj  (-i)k+z 

= 1 k = 0 


1*3*5  (2k  + 21  - 1) 

(H-k) ! 2l+k  k! 


(cos  I)2k~2m 
m!  (k-m) ! 

5 


(cos  n)2m  (sin  n)2k_2m 


2 (r)2£_1S  (_1)k+* 

l = 1 P k = 1 


1*3-5  (2k  + 21  -1) 

U-k)!  2X'+kk! 


Icos  i)2k~2nl  £2k-2m)  (cos  n)2m+1  (sin  n)2^21""1 
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-2m(cos  n) 2m_1  (sin  n)2k"2m+1] 


U sin  I 
m 


2 
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CO 
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(i- k)!  25,+k  k! 


(2k  - 2m) 
m!  (k-m)! 
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When  Equations  (2-1),  (2-2),  and  (2-3)  are  substituted  into  Gauss 
form  of  the  equation  of  the  osculating  elements  and  extensive  opera- 
tions are  performed.  Ash  finds  the  change  of  right  ascension,  eccen- 
tricity, and  argument  of  perigee.  Change  in  semi-major  axis  and 
inclination  were  deemed  negligible. 

The  changes  in  right  ascension,  eccentricity,  and  argument  of 
perigee  per  orbit  were  found  to  be: 


AW  — 2tt  (— ) (a)3  cos  I r + (-) 2 
— 4 p 

U P 


135  . 315  2 

128  128  COS  1 
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The  changes  represented  by  Equations  (2-4),  (2-5)  and  (2-6)  can 
be  very  useful  for  long  term  (secular)  changes  in  orbits  as  a function 
of  lunar  perturbations.  The  computational  order  which  should  be 
followed  is: 


1.  4#^  should  be  computed  from  Equation  (2-6)  from  an 

initial  value  of  argument  of  perigee  W . A new  value 

o 

of  W is  then  found  from  (W  + W ).  Sometimes,  how- 
P op' 

ever, is  identically  zero  and  a stable  argument  of 
perigee  exists;  this  can  happen  at  high  inclination 

angle  (XI)  and  large  semi-major  axis  (a).  Program 
PERTP,  Appendix  5,  calculates  the  stability  of 
argument  of  perigee  in  a subroutine  ARGPER. 

2.  Wp  is  then  substituted  into  Equation  (2-5)  to  find  A e 
per  orbit.  This  is  done  in  PERTP  on  line  630  (App.  5). 

3.  The  change  of  right  ascension  is  calculated  from  Equation 
(2-4).  This  is  done  in  PERTP  on  line  620  (App.  5). 

PERTP  is  also  more  comprehensive;  this  will  be  shown  in  the 
following  estimates  of  change  in  semi -major  axis. 
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2.2  ESTIMATED  CHANGE  IN  SEMI-MAJOR  AXIS 

A surprising  difference  in  long  term  stability  of  posigrade 
compared  to  retrograde  orbits  was  identified  by  Ash^.  For 
example,  12-day  retrograde  orbits  appear  approximately  as  stable  as 
10-day  posigrade  orbits  over  a two-year  interval.  This  difference 
in  stability  does  not  appear  in  Equations  (2-1) , (2-5) , and  (2-6) 
after  the  development  of  perturbations  due  to  a torus  of  lunar 
matter. 

Clearly,  at  least  one  major  physical  concern  has  been  deleted 
from  the  development  for  Equations  (2-4),  (2-5),  and  (2-6).  This 
is  the  interaction  time  between  the  satellite  and  the  moon,  which 
would  show  a difference  between  posigrade  and  retrograde  orbits. 
Figure  9 gives  a planar  diagram  of  lunar  perturbation. 

If  the  satellite  is  in  a circular  orbit  about  the  earth,  it 
becomes  convenient  to  consider  the  tangential  impulse  and  the  radial 
impulse  imparted  to  the  satellite  from  the  moon.  By  symmetry,  it  can 
be  seen  that  the  satellite  loses  about  as  much  tangential  impulse  as 
it  passes  the  moon  as  when  it  approaches.  However,  the  impulse 
along  the  radial  direction  is  always  positive. 


Radial  force  = p = GMm  (cos 

per  unit  mass  r m 


(2-7a) 


where 


r^  = distance  between  satellite  and  lunacenter 

G = gravitational  constant 

R ~ 

9 . s for  9 < 0.15  rad 

m e (R  -R  ) 6 

m s 


(2- 7b) 


M = lunar  mass 
m 
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W.  RAO/SEC 


9 GEOMETRY  FOR  LUNAR  PERTURBATIONS 


=.3 


Further 


R -R 

r . . can  be  noted  to  be  m a 

cos  6 


m 


for  e 


e <0.15  rad,  and  Equation  (2-7a)  becomes 


GM 


F i ™ 

r (R  -R.) 

m s 


cos3  f 9e  .(  — 


R -R  ) 
m s 


(2-8) 


The  first  order  change  in  radial  velocity  is  related  to  radial 
impulse  by 


total  time 


Av  = 


' / 


F dt 
r 


(2-9) 


Equations  (2-8)  and  (2-9)  can  be  integrated,  yielding  a closed 
form  solution,  and  this  opportunity  to  save  CPU  time  will  be  adopted. 

If  an  angle  6eo  is  chosen  as  appropriate  to  the  position  at  which 
noticeable  lunar  perturbations  begin,  the  time  until  the  earth, 
satellite,  and  moon  are  lined  up  (maximum  lunar  perturbation  force) 
is  9eo(w  -uj.  Using  the  small  angle  approximation.  Equation  (2-7b), 
and  substituting  Equation  (2-8)  into  Equation  (2-9)  yields: 


1 


2 GM 


m 


^ Vr  ut  , (R  -R  ) 
rel  ms 


1/3  sin  (0 


eo 


s ) 


R -R 
m s 


2 

cos  (e 


eo 


s ) + 2 


R -R 
m s 


(2-10) 


where  = relative  angular  velocity  of  satellite  and  moon.  This 

is  shown  explicitly  later  in  Equation  (2-16). 
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Equation  (2-10)  has  been  implemented  in  subroutine  MOON  of  PERTP 
(Appendix  5).  The  limit  of  integration,  0eo,  is  entered  as  TE  » 0.3 
radians  on  line  1210.  This  is  felt  to  give  an  underestimate  of  lunar 
perturbations,  because  a tidal  displacement  actually  occurs  during 
the  radial  impulse.  The  radial  velocity  of  Equation  (2-10)  must  be 
related  to  orbital  energy,  and  hence  to  semi— major  axis  (a)  as  soon 
as  the  satellite  leaves  the  influence  of  the  moon. 


The  vis  viva  integral  is  the  relation  between  orbital  energy 
and  semi-major  axis  for  a simple  central  force  field  and  has  the 
form 


where 


for  elliptic  motion 


V = total  velocity 


H-  Gm 

e 

m - mass  of  earth 
e 


(2-11) 


R = geocentric  distance. 


When  Eauation  (2-11)  is  rewritten 

vi  - 4 - 4 (-  I \ 

2 r 2 \ a / , 

the  left  side  is  recognized  as  total  orbital  energy,  or 


f 


45 


- 77  (total  energy)  = — 
m a 


After  rearranging  Equation  (2-13)  and  differentiating  the  expression 
of  (a) , a function  of  (total  energy) , one  finds 


= , £ ) d (total  energy) 

2 (total  energy) 


(2-14) 


Changing  Equation  (2-14)  to  an  incremental  form  and  substituting 
2 

AV r for  change  in  total  energy  gives 


2 2 

/)a  = CAVr>  a 


(2-15) 


Since  is  proportional  to  interaction  time,  the  square  of 
in  Equation  (2-15)  implies  that  the  instability  of  the  semi- 
major axis  increases  in  a nonlinear  way  with  the  lunar  interaction 
time.  The  actual  interaction  time  depends  on  a vector  difference  of 
angular  velocities,  ( CJ  - <j  ).  The  relative  angular  velocity  can  be 
approximated  by 


2 2 

U>  = (u  cos  I - <j  ) + (cl>  sin  I) 

rel  s ms 


(2-16) 


where  again  I is  the  inclination  angle  of  the  satellite  with  respect 
to  the  lunar  plane. 

The  remaining  considerations  which  have  gone  into  the  subroutine 
MOON  are: 
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1.  The  total  region  of  possible  Interaction  of  the  satellite 
with  the  moon  is  approximately  limited  (to  ANG)  if  the 
inclination  angle  is  large  and. 


2.  Change  of  (a)  per  year  is  proportional  to  the  number  of 
times  per  year  in  which  interaction  happens,  or /la/ 

yr 

a(Jrel  * 

Even  with  the  second  item  included  for  completeness,  however, 
a distinct  difference  between  the  stabilities  of  posigrade  and  retro- 
grade orbits  is  seen  in  only  a few  seconds  of  CPU  time.  When  the 
effects  of  subroutine  MOON  are  included  in  PERTP  and  a variable  to 
represent  instability 


a (1  + e)/a  (1  + e ) 
o o 


(2-17) 


is  generated  to  represent  normalized  apogee  distance,  a 10-day  orbit 
at  inclination  0°  has  approximately  the  same  value  for  V5  (1.003) 
as  a 12-day  orbit  at  inclination  180°. 

It  appears  that  the  gross  secular  results  for  orbital  insta- 
bility^ have  been  qualitatively  reproduced  with  very  modest  CPU 
requirements. 

The  reader  will  notice  that  a running  change  in  units  has 
occurred  in  these  programs.  At  the  time  these  lunar  perturbation 
programs  were  generated,  a new  Federal  interest  was  expressed  in 
metric  units.  Semi-major  axis  (a)  is  expressed  in  the  program  of 
Example  2 and  Appendices  5-7  in  kilometers.  The  relation  between 
nautical  miles  and  kilometers  is: 

1.852  x (length  in  nautical  miles)  ■ length  in  kilometers. 
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2.3  PROGRAM  PERTP  (Appendix  5) 

Program  PERTP  calculates  a new  semi-major  axis  (a,  km), 
eccentricity  (E) , right  ascension  (W) , and  normalized  apogee 
distance  V^.  It  does  this  from  estimated  lunar  secular  perturba- 
tions. Iterations  in  orbital  elements  are  performed  every  10  days 
for  a total  of  400  days. 

CALCOMP  plots  are  generated  (e.g.,  vs  time  in  days). 

This  program  gives  an  estimate  of  orbital  stability  implicit 
in  the  variable  V^.  When  V,.  increases  to  approximately  1.01  in 
400  days,  the  orbit's  stability  is  very  questionable.  This  stability 
has  qualitative  agreement  with  Ash's  numerical  results. 

2.4  PROGRAM  SATLUNAE  (Appendix  6) 

Program  SATLUNAE  combines  a variation  of  SATE  and  a subroutine 
LUNA.  To  save  CPU  time  for  orbits  which  are  to  be  perturbed  by 
the  moon  for  a long  time  (say,  T > 3 years),  integrated  versions 
of  Equations  (2-4),  (2-5),  and  (2-6)  were  prepared  for  subroutine 
LUNA.  The  equation  for  argument  of  perigee  offered  two  choices: 
either  (1)  it  kept  rotating,  or  (2)  it  was  assigned  a stationary 
value  immediately  as  in  the  ARGPER  subroutine  of  PERTP.  If  it 
kept  rotating,  an  average  value  of  angular  velocity  was  assigned 
by  using  si^  u * 1/2. 

The  order  of  calculations  in  LUNA  was  as  follows: 


1.  A linear  estimate  of  growth  in  semi-major  axis  (a) 
was  done,  so  that  A (AB  in  the  program)  could  be  used 
for  later  calculations . 

2.  Argument  of  perigee  was  estimated  as  either  stationary 
or  rotating,  which  gave  the  possibility  of  two  integrals 
for  the  eccentricity  calculation. 


3.  If  a stationary  existed,  a simple  exponential 

growth  in  eccentricity  resulted.  This  is  similar  to 


the  result  of  the  1974  Ash  report 


A less  simple 


result  occurred  (on  line  3020)  if  u>  kept  rotating. 

P 

Lines  3020  and  3030  are  deceptive  in  their  simplicity, 
and  to  see  the  complexity  one  should  substitute  the 
various  c’s.  Many  of  the  c's  are  not  constants,  but 
are  functions  of  the  orbital  elements. 


Right  ascension  and  normalized  apogee  distance  are  finally  esti- 
mated. The  results  show  a satisfying  but  not  exact  comparison  to 
an  iterated  (not  integrated)  change  in  orbital  elements.  A pro- 
gram which  iterates  the  orbital  elements  appears  as  SCOREE  (Appendix  7), 
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PART  III 
EXAMPLES 

3.1  EXAMPLE  1 

Figure  10  gives  an  illustration  of  the  10  satellites  chosen  for 
this  example;  they  are  a running  example  throughout  this  section. 

These  10  satellites  are  examined  here  with  program  SATE.  They  appear 
again  in  Example  2,  where  lunar  perturbations  are  examined,  and  in  a 
coverage  problem  of  Example  3.  The  10  satellites  give  an  idea  of 
the  capability  of  these  programs  and  are  not  intended  as  a specifically 
useful  satellite  communications  system.  Satellites  #1  and  #2  are 
Aerospace  Type  satellites  selected  to  give  polar  coverage. 

These  12-hour  orbits  differ  in  right  ascension  (W)  by  90°  in  order 
to  give  the  same  ground  trace.  Satellites  //3,  #4,  #5,  and  #6  are 
synchronous  equatorial  satellites.  Satellites  // 7,  //8,  #9,  and 
//10  are  in  4-day  polar  orbits.  The  orbital  elements  are  entered 
in  program  SATE  on  lines  160  to  240  (Table  1).  The  ground  station 
coordinates  are  entered  on  line  250;  Lexington,  Massachusetts 
and  San  Diego,  California,  are  entered  as  ground  stations  1 and  2„ 

Lines  160  and  170  list  the  semi-major  axes  of  the  satellites  in 
nautical  miles,  in  order  of  increasing  altitude.  Line  180  gives  the 
respective  eccentricities.  Lines  190  and  200  give  the  right  ascen- 
sions in  degrees.  Line  210  gives  arguments  of  perigee  in  degrees. 

Lines  220  and  230  give  the  time  at  perigee  in  hours,  and  line  240 
gives  orbital  inclination  in  degrees.  Line  250  gives  colatitudes 
of  Lexington,  Massachusetts,  and  San  Diego,  California;  then,  the 
longitudes  of  the  same  sites. 
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TABLE  1 


ORBITAL  ELEMENTS  ENTERED  IN  SATE 


00150 

Dimension  A(10) , E(10),  W(10) , WP(10),  TP(10), 

XI  (10), 

B(4) 

00160 

Data  A/14342.,  14342.,  22767.,  22767.,  22767., 

22767., 

00170 

3 57369.2,  57369.2,  57369.2,  56369.2/ 

00180 

Data  E/ .725,  .725,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

0 ./ 

00190 

Data  W/o.,  270.,  0.,  0.,  0.,  0., 

00200 

• 

o 

• 

o 

• 

o 

•* 

• 

o 

00210 

Data  WP/-90 . , -90.,  0.,  0.,  0.,  0.,  0.,  0.,  0., 

0./ 

00220 

Data  TP/0.,  -6.,  -3.,  -9.,  -15.,  -21., 

00230 

1 0.,  -24.,  -48.,  -72./ 

00240 

Data  xI/63.435,  63.435,  0.,  0.,  0.,  0.,  90.,  90.,  90., 

90./ 

00250 

Data  B/47.54,  57.23,  288.73,  242.8/ 

As  an  important  aside,  it  should  be  noted  that  the  communica- 
tions designer  does  not  really  care  as  much  about  orbital  elements 
as  he  does  about  the  subsatellite  traces  which  he  will  require  in 
order  to  get  adequate  satellite  visibility.  The  designer  can  tell 
that  Satellite  ill  was  at  -90°  longitude  (90°  West  longitude)  at 
T =*  0 hrs  because  TP  * 0,  WP  = -90°,  and  W » 0.  The  subsatellite 
"starting  point"  of  a satellite  with  TP  ^ 0 requires  some  calculation. 
Satellite  # 2,  with  TP  **  -6  hrs,  was  at  perigee  6 hours  before  the 
computation  began.  It  is  therefore  halfway  through  its  12-hour 
orbit  at  T * 0 hrs,  and  is  at  apogee.  Since  the  right  ascension 
is  270°  and  the  argument  of  perigee  is  -90°  for  Satellite  #2,  the 
perigee  occurred  at  180°  longitude  in  celestial  coordinates. 

However,  perigee  occurred  6 hours  previously  for  Satellite  it 2; 
the  earth  advanced  1/4  revolution  in  that  period,  or  perigee 
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occurred  over  270°  East  longitude  on  the  earth.  Since  this  is 
involved,  a simple  check  for  the  communications  engineer  is 
available  on  the  computer  output  which  gives  the  subsatellite 
trace. 


The  computer  output  of  SATE  gives  the  time  history  of  all  combi- 
nations of  links.  It  is  not  a graphical  output,  but  a compensatory 
feature  is  its  28-second  CPU  time  on  a 370/158  for  25  hours  of  data. 

A sample  link  from  Satellite  #1  to  Lexington,  Massachusetts,  is 
shown  on  Table  2.  From  left  to  right,  the  columns  are  satellite 
number,  time  in  hours,  range  in  nautical  miles,  elevation  angle  in 
degrees,  uplink  Doppler  shift  (Hz)  for  a 300  MHz  signal,  downlink 
Doppler  shift  (Hz)  for  a 245  MHz  signal,  subsatellite  longitude, 
subsatellite  colatitude,  and  ground  station  number.  It  is  seen 
that  range  increases  until  a maximum  is  reached  at  6 hours  for  the 
first  apogee. 

No  earth  harmonics  for  the  potential  field  have  been  intro- 
duced for  Table  2,  so  this  near  earth  orbit  should  not  be 
examined  beyond  two  orbits  or  24  hours. 

The  55  possible  crosslinks  are  also  listed  in  the  computer  output. 

A crosslink  from  Satellite  //I  to  Satellite  //3  is  examined  in  Table  3. 
This  represents  a highly  eccentric  Aerospace  Type  A orbit  crosslinking 
to  a synchronous  satellite.  From  left  to  right,  time  is  given  in  hours, 
range  in  nautical  miles,  L12  and  C12  are  the  longitude  and  colatitude 
of  the  crosslink  pointing  vector  (in  the  inertial  coordinate  system) , 
crosslink  Doppler  (Hz)  is  given  for  the  60  GHz  crosslink1  , and  CVTS 
gives  a visibility  check  on  the  crosslink.  CVIS  = 1 means  the  cross- 
link exists;  CVIS  = 0 if  the  earth  blocks  the  crosslink.  The  first 
maximum  for  crosslink  range  is  seen  to  occur  at  6 hours  and  27822 
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nautical  miles.  Crosslink  Doppler  peaks  at  1.39  MHz  at  12  hours. 

A more  convenient  and  intuitive  crosslink  pointing  system  will  be 
given  later  in  program  SATLUNAE. 

These  results  were  given  in  single  precision.  Because  elevation 
angle  was  required  only  within  0. 3°  and  earth  harmonic  terms  were 
omitted,  greater  accuracy  was  not  felt  to  be  justified  here.  Single  pre- 
cision can  lead  to  a strange  result  for  crosslink  Doppler  between  two 
coplanar  satellites  in  circular  orbit;  the  required  zero  Doppler  shift 
is  calculated  to  be  a fraction  of  a Hz  because  it  is  the  result  of  a 
subtraction  of  two  large  numbers.  Another  embarrassment  for  single 
precision  can  occur  in  crosslink  pointing  angles.  The  tangent  of 
angles  near  90°  may  be  read  as  tan  (90°).  If  these  Inconveniences  are 
to  be  avoided,  the  double  precision  program  SATD  (Appendix  2)  can  be 
used. 
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3.2 


EXAMPLE  2 


The  orbital  configuration  of  Figure  10  is  used  again  in  Example  2, 
but  this  time  the  stability  of  the  orbits  relative  to  lunar  perturba- 
tions is  to  be  checked.  The  orbital  elements  of  Figure  10  are  entered 
into  program  SATLUNAE  (Appendix  6) , but  unlike  the  previous  example 
the  semi-major  axis  is  entered  in  km.  The  initial  semi-major  axes 
are  26561  km,  42164  km,  and  106247  km  respectively  for  the  12-hr, 

24-hr,  and  4-day  orbits. 

The  first  part  of  the  output  of  SATLUNAE  is  an  estimate  of  new 
orbital  elements  after  a long  period  of  lunar  perturbation.  Although 
SATLUNAE  is  meant  for  very  long  periods,  like  five  years,  365  days  is 
the  period  used  for  this  example.  The  estimate  for  the  new 
orbital  elements  after  one  year  is  shown  in  Table  4a. 

The  semi -major  axis  for  the  12-hour  satellites  is  seen  to  be  the 
same  as  the  initial  value,  but  the  eccentricity  has  increased  from 
0.725  to  0.745.  This  jump  in  eccentricity  is  due  to  the  entry  of  the 
fourth  column,  WPS.  WPS  is  a stationary  argument  of  perigee  which  was 
set  immediately  in  the  computer  program  SATLUNAE.  This  stationary 
argument  of  perigee  allows  the  eccentricity  from  the  program  to  start 
growing  immediately,  although  real  satellite  eccentricity  does  not 
immediately  grow.  Therefore  the  new  eccentricity  El  should  be  con- 
sidered an  upper  bound.  The  new  right  ascension  W is  only  slightly 
changed  from  the  original  0°  for  Satellite  #1.  The  normalized  apogee 
distance  has  increased  1.2%  in  one  year. 
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TABLE  4a 


FUTURE  SEMI-MAJOR  AXIS,  ECCENTRICITY,  RIGHT  ASCENSION, 
ARGUMENT  OF  PERIGEE,  AND  NORMALIZED  APOGEE  DISTANCE  V5 


ABIG 

El 

W 

WPS 

V5 

26561.0 

0.745798409 

-0.360600 

44.963928 

1.0120573 

26561.0 

0.745798409 

269.639160 

44.963928 

1.0120573 

42164.2 

0.0 

-1.655902 

0.0 

1.0000038 

42164.2 

0.0 

-1.655902 

0.0 

1.0000038 

42164.2 

0.0 

-1.655902 

0.0 

1.0000038 

42164.2 

0.0 

-1.655902 

0.0 

1.0000038 

106249.1 

0.0 

-0.000002 

39.585815 

1.0000191 

106249.1 

0.0 

-0.000002 

39.585815 

1.0000191 

106249.1 

0.0 

-0.000002 

39.585815 

1.0000191 

106249.1 

0.0 

-0.000002 

39.585815 

1.0000191 

Surprising  stability  appears  in  the  V,.  column  for  Satellites 
#3  through  #10.  This  is  because  the  initial  eccentricity  was  0. 

A more  reasonable  estimate  of  stability  can  be  found  if  a non-zero 
value  for  initial  eccentricity  is  used.  Table  4b  lets  the  initial 
eccentricity  value  be  0.05. 

The  synchronous  Satellites  #3  through  #6  are  still  stable  since 
retains  a value  close  to  unity.  However,  it  is  seen  that  the  4-day 
polar  satellites  show  a 1.4%  increase  in  normalized  apogee  distance. 
Questionable  stability  is  indicated. 
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TABLE  4b 


SATELLITES  #3  THROUGH  #10,  INITIAL  ECCENTRICITY  =*  0.05 


ABIG 

El 

W 

WPS 

V5 

26561.0 

0.745798409 

-0.360600 

44.963928 

1.0120573 

26561.0 

0.745798409 

269.639160 

44.963928 

1.0120573 

42164.2 

0.050000001 

-1.655902 

0.0 

1.0000038 

42164.2 

0.050000001 

-1.655902 

0.0 

1.0000038 

42164.2 

0.050000001 

-1.655902 

0.0 

1.0000038 

42164.2 

0.050000001 

-1.655902 

0.0 

1.0000038 

106249.1 

0.065072298 

-0.000002 

39.585815 

1.0143747 

106249.1 

0.065072298 

-0.000002 

39.585815 

1.0143747 

106249.1 

0.065072298 

-0.000002 

39.585815 

1.0143747 

106249.1 

0.065072298 

-0.000002 

39.585815 

1.0143747 

SATLUNAE  then  goes  through  link  calculations  with  the  new  orbital 
elements.  These  are  largely  omitted  here,  since  they  are  reminiscent 
of  the  output  from  SATE.  One  other  crosslink  calculation  does  occur 
here,  however,  which  makes  it  worthwhile  to  include  one  crosslink. 
Crosslink  angle  in  a local  satellite  coordinate  system  is  shown  in 
Figure  11  and  it  is  included  as  TU12  and  FU12  in  Table  5;  this  table 
is  part  of  one  output  table  of  SATLUNAE. 

Except  for  the  columns  FU12  and  TU12,  Table  5 was  discussed  in 
Example  1.  TU12  and  FU12  are  discussed  in  Figure  11.  Again,  TU12 
is  a local  co-elevation  angle  and  FU12  is  a local  azimuth  angle. 
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»A-  45,124 


THE  LOCAL  ZENITH  ANGLE  9ai  „ IS  COMPUTED  IN  SUBROUTINE  UNPRIM 
AND  PRINTED  OUT  AS  TU,*  IN  PROGRAM  SATLUNAE.  THE  LOCAL  AZIMUTH 
ANGLE  (PROJECTION  IN  THE  XUYUPLANE,  MEASURED  FROM  THE  Xu  AXIS) 
IS  PRINTED  AS  FU,t 


n 

Figure  II  LOCAL  COORDINATE  SYSTEM  FOR  SATELLITE  CROSSLINK 


CROSSLINK  RELATIONS  FROM  PROGRAM  SATLUNAE 
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Program  SATLUNAE  can  also  be  run  repeatedly  to  get  entire  regions 
of  stability.  An  example  of  the  stability  analysis  appears  In 
Figure  12,  which  gives  a topology  of  stability  that  results  from  an 
initial  eccentricity  equal  to  0.005.  The  curves  should  be  read  as  a 
contour  map  with  the  orbital  instability  increasing  to  the  right. 

Each  curve  represents  a constant  stability  line.  The  curve  labelled 
0.3%  change  in  normalized  apogee  distance  per  year  shows  that  10-day 
posigrade  (I  ■ 0°)  orbits  approach  the  same  stability  as  8-day  polar 
orbits  and  12-day  retrograde  (I  » 180°)  orbits.  On  Figure  13  it  is 
seen  that  initial  eccentricity  is  very  important  to  the  orbital 
stability  of  satellites  outside  the  lunar  plane.  The  initial 
eccentricity  (0.05)  of  Figure  13  gives  greatly  increased  instability 
at  I = 90°  as  compared  to  Figure  12.  Since  stability  of  these  orbits 
is  a strong  function  of  initial  eccentricity,  launch  techniques  pro- 
ducing low  initial  eccentricity  should  be  employed  for  high  altitude 
orbits. 
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ESTIMATED  AVERAGE  CHANGE  IN  NORMALIZED  APOGEE  DISTANCE  PER  YEAR 

Lo,(i+.,)  J 

(AVERAGED  OVER  5 YEARS)  AS  A FUNCTION  OF  LUNAR  PERTURBATIONS  e0=  .005 


3.3  EXAMPLE  3 

Two  of  the  initial  problems  which  face  a satellite  communications 
designer  are:  (1)  Given  a set  of  satellites,  at  what  points  on  the 

earth  can  an  observer  see  a satellite?  (2)  If  many  satellites  are 
involved,  at  what  spots  on  earth  are  no  satellites  visible?  Preferably, 
the  designer  should  see  a graphical  output  of  covered/uncovered  areas. 

In  seeking  an  answer  to  (1),  one  might  get  swamped  by  information 
from  the  180°  x 260°  projection  which  forms  the  Mercator  projection. 

The  really  interesting  locations  in  a good  satellite  system  might  be 
those  which  can  see  no  satellite  at  all.  Identification  of  these 
uncovered  locations  yields  a fast,  comprehensible  output.  It  is 
with  this  idea  of  limited  output  that  program  NOLINKE  (Appendix  9) 
was  generated.  With  10  satellites  in  Keplerian  orbits,  it  examines 
a 12  x 24  grid  of  possible  ground  stations  (covering  the  180°  x 
360°  Mercator  projection  in  15°  increments)  for  elevation  angle  to 
all  ten  satellites.  It  performs  a printout  at  a given  coordinate 
only  if  elevation  angle  to  each  of  the  10  satellites  is  less  than 
some  minimum  acceptable  elevation  angle  (EM) . 


* 

. 
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The  lines  10  to  70  of  NOLINKE  (Appendix  9)  represent  the 
orbital  elements  of  10  satellites.  The  four  DO  LOOPS,  progressing 
from  outside  to  inside,  are  time  (RI),  colatitude  (T),  longitude 
(F) , and  satellite  number  (L) . So,  at  each  (T,  F)  all  satellites 
are  considered.  Elevation  angle  (D)  to  each  satellite  is  generated 
at  each  location.  A print  statement  is  allowed  only  if  all 
satellites  have  been  considered  and  if  no  satellite  is  visible 
(SATNO  = 0). 


This  limited  output  can  be  very  helpful  to  the  designer  of  a 
many  satellite  system. 


i 
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The  program  NOLINKE  runs  in  the  foreground  on  TSO.  The  limited 
output  allows  this  to  be  done  conveniently.  For  example,  the  orbital 
elements  implied  by  Figure  10  can  be  typed  into  lines  20  to  70,  a 
minimum  elevation  angle  * 35°  typed  into  line  210,  and  the  resultant 
output  (at  T - 0 hrs,  1 hr)  is  shown  on  Table  6.  T1  and  FI 
are  ground  station  colatitude  and  longitude  at  which  no  satellite  is 
observed.  No.  is  the  number  of  satellites  observed  at  these  loca- 
tions. T is  time  in  hours. 

TABLE  6 

AN  OUTPUT  OF  NOLINKE 


T1 

FI 

NO. 

T 

60.000 

90.000 

0.0 

0.0 

60.000 

270.000 

0.0 

0.0 

120.000 

90.000 

0.0 

0.0 

120.000 

270.000 

0.0 

0.0 

60.000 

90.000 

0.0 

1.000 

60.000 

270.000 

0.0 

1.000 

120.000 

90.000 

0.0 

1.000 

120.000 

270.000 

0.0 

1.000 

Although  only  four  locations  appear  at  each  hour,  12  x 24  ■ 288 
calculations  were  done  worldwide  to  see  if  elevation  angle  require- 
ments were  met.  The  locations  and  times  of  Table  6 were  the  only 
places  which  did  not  meet  the  required  35°  elevation  angle. 
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No  output  occurs  for  NOLINKE  for  a 6-hour  period  if  a 25°  eleva- 
tion angle  is  required.  Therefore,  this  10-satellite  system  meets  a 
25°  minimum  elevation  angle  requirement. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


An  analysis  has  been  presented  which  gives  fast,  convenient  com- 
putations for  range,  Doppler,  crosslink  data,  and  pointing  angles  for 
a two  ground-station,  10- satellite  system.  With  the  programs  listed, 
data  on  all  ways  to  link  one  ground  station  with  the  other  via  the 
satellite  are  given  on  one  computer  output.  The  writers  are  not  aware 
of  the  existence  of  other  complete  link  solutions  on  one  output. 

Three  place  accuracy  is  representatively  given.  Completely  general 
Keplerian  orbits  are  allowable  with  0 < eccentricity  < 0. 99. 

The  coverage  of  a 10- satellite  system  is  quickly  found  in  a 
program  with  very  limited  output. 

In  continuing  the  emphasis  on  low  computation  time,  analytic 
results  for  orbital  stability  as  a function  of  lunar  perturbations 
were  used  to  generate  a very  efficient  program.  It  can  be  used  as  a 
first  economical  estimate  of  high  altitude  orbital  parameters  after 
many  years  of  lunar  perturbations.  Figures  12  and  13  were  also 
included  to  indicate  stability  regions. 

None  of  the  listed  FORTRAN  programs  requires  more  than  one  minute 
CPU  time  on  an  IBM  370/158  computer. 

In  the  future,  it  is  recommended  that  real  signal  propagation  be 
interposed  between  the  stations.  This  could  include  ideal  antenna  gain 

patterns  and  regions  of  ionospheric  disturbance  to  yield  estimates  of 
the  received  signal/noise  ratio.  A "point  ahead"  capability  should 
also  be  added  to  the  crosslinks;  this  can  be  done  conveniently  with 
the  existing  velocity  calculations. 
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NOTATION 


Symbol 


Definition 


Comment 


semi-major  axis;  Invariant  for 
Keplerian  orbits 


Nautical  miles, 
except  in  pro- 
grams for  lunar 
perturbations , 
where  a is  en- 
tered in  km. 


average  semi -major  axis 


initial  semi-major  axis  (be- 
fore lunar  perturbations) 


1 nmi  - 1.852  km 

Subscript  1 is 
used  to  denote 
satellite  1,  etc. 

Average  taken  over 
a period  of  lunar 
perturbation. 


all"a33 


elements  of  Euler  rotation 
matrix 


future  semi-major  axis,  subject 
to  lunar  perturbations 

Euler  rotation  matrix 


inverse  of  A 

geocentric  angle  between  sub- 
satellite point  and  ground 
station  1 

FORTRAN  program  which  includes 
azimuth  calculations  (App.  8) 

semi-minor  axis 


Not  used  here. 


velocity  of  light  = 1.6187  x 
105  nmi /sec 


. r PKECjQIrG  ' BLANK- NOT  FILMED 


elevation  angle  from  ground 
station  to  satellite,  degrees 


denotes  small  changes  in  orbi- 
tal elements  (e.g.  4e=small 
change  in  eccentricity) 


eccentricity;  for  the  ellipti- 
cal Keplerian  orbits  considered 
here  0 < e < 1. 


The  programs  with 
the  iterative  sub- 
routine ELL IP  work 
for  0<e<0.99 


future  eccentricity,  subject  to 
lunar  perturbations 

eccentric  anomaly  E^,  E ...  used 

for  successive 
approximations  to 
eccentric  anomaly. 

carrier  frequency,  Hz 

Doppler  shifted  carrier  fre- 
quency, Hz 

lunar  radial  force/unit  mass  of 
satellite  (scalar) 

angle  measured  from  the  line  of 
nodes  in  the  satellite  plane 

same  as  FSP  = + ^ 


gravitational  constant  = 6.6732 
x 10”H  newton  meter  2/kg^ 


orbital  inclination  angle  with 
respect  to  equatorial  plane 


Becomes  XI  in 
FORTRAN  programs. 
Numbered  subscripts 
refer  to  satellite 
number. 


unit  vectors  in  x,  y,  and  z 
directions,  respectively 

orbital  inclination  with  re- 
spect to  lunar  plane 


M 

M 

e 

M 

m 

^m 

nautical 
mile  (nml) 

n 

NOLINKE 


P 

PERTP 

^12 

R 

R 

R 

e 

*12 


mean  anomaly 
mas 8 of  earth 
mass  of  moon 
GM 

e 

(31 

m 

1.852  km  (exact) 


mean  angular  rate,  rad/hr  M 

a coverage  program  of  very 
limited  output  for  ten  satel- 
lites (App.  9) 

angle  in  satellite  plane  measured 
from  line  of  nodes  in  the  lunar 
plane 

semi-latus  rectum 

FORTRAN  program  which  cal- 
culates effects  of  lunar  per- 
turbations on  orbital  elements 
(App.  5) 

unit  vector  from  transmitter 
to  receiver 

geocentric  distance  of  satel- 
lite at  a given  instant  of  time 
(nmi  for  programs  without  lunar 
perturbations,  km  for  those 
with) 

d (R) 
dt 

earth  radius  ■ 3440  nmi 

distance  from  transmitter  to 
receiver 


- n(t-tp) 
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£ r 

- /-I 

R1S 

distance  from  ground  trans- 
mitter 1 to  satellite 

RS1,  S2 

distance  from  Satellite  //I  to 
Satellite  \\2 

1- 

RI 

time  index  (hrs)  in  programs  Note  that  solar 

(T  in  analysis)  hours  are  not 

distinguished 
from  sidereal 
hours  because  only 
three-place  accuracy 
is  desired  in  re- 
sults. 

R,  S,  W 

lunar  perturbation  components 
in  Ash  analysis 

SATD 

FORTRAN  program  which  calculates 
ground-satellite  relations 
(e<0.5)  in  double  precision 

(App.  2)  r 

SATE 

FORTRAN  program  which  cal- 
culates ground-satellite  rela- 
tions (0<e<0.99)  in  single  pre- 
cision (App.  1) 

SATLUNAE 

like  SATE,  but  includes  long  term 
estimates  (>5  years)  of  lunar  per- 
turbations (App.  6) 

f i 

SATVIZD 

like  SATD,  but  interactive  with 
the  user  (App.  4) 

V 

SATVIZE 

like  SATE,  but  interactive  with 
the  user  (App.  3) 

SCOREE 

like  SATE,  but  includes  effects 
of  lunar  perturbations  on  orbi- 
tal elements  for  periods  up  to 
a few  years  (App.  7) 

T 

time  in  hours  RI  in  programs 

TP 

time  at  perigee,  hours 

Lu 

74  1 

J 

r period  for  Keplerian  orbit 

(subscripts  are  added  for 
particular  satellites) 

TU12,  FU12  crosslink  angles  in  a local 

satellite  coordinate  system 

6 true  anomaly 


6 


d ( ) 
dt 


e 

SI,  S2; 

$ 

SI,  S2 


9 

eo 


6 

m 


angle  between  relative  velocity  vector 
and  pointing  vector 

spherical  coordinates  associated 
with  ground  station  1 


inertial  spherical  coordinates 
associated  with  crosslink  vec- 
tor (called  L12,  C12  in  com- 
puter output) 


TU12,  FU12  are 
usually  more  con- 
venient than 
0 • d) 

SI,  S2’  SI,  S2. 


angle  of  satellite  with  respect 
to  earth-moon  vector 


bounding  angle  for  lunar  pertur- 
bations (representatively,  0.3 
radians) 


angle  of  satellite  with  respect 
to  moon-earth  vector 


relative  velocity  between  trans- 
mitter and  receiver 


V 


r 

V 


satellite  radial  velocity 

potential  energy  in  v is  viva  Integral 

normalized  apogee  distance  » final 
apogee/ initial  apogee 
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vis  viva  integral  gives  relation  between  potential  and 

kinetic  energy  for  Keplerian  orbits 


W 


right  ascension  for  Keplerian  Subscripts  are  used 
orbit,  degrees  for  particular 

satellites. 


W argument  of  perigee  for  Keplerian 

orbit,  degrees 

WPS  stationary  argument  of  perigee, 

subject  to  lunar  perturbations 

(j  earth  rotation  rate  - 15°/solar 

hour  = 15.04°/sidereal  hour 


U 

s 


CJ 

m 


u>  . 
rel 


W 


X,Y,Z, 


X'.Y'.Z', 


angular  rotation  rate  of  satellite 
about  earth,  rad/hr 

angular  rotation  rate  of  moon 
about  earth,  rad /hr 

vector  difference  of  u , u u , determines 

8 m rel 

the  time  of  in- 
teraction of  satel- 
lite and  moon. 


one  of  three  lunar  perturbation 
components  used  by  Ash  (others, 
R,  S) 

inertial  Cartesian  coordinate 
system,  with  x pointed  toward 
Aries  (X,Y,Z,  associated  with 
ground  station  1) 

a moving  coordinate  system  in 
the  satellite  plane 
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X",Y",Z",  a stationary  coordinate  system  The  double  prime 

In  the  satellite  plane  system  will  be  con- 

venient if  Doppler 
rates  are  later  added 
to  the  programs. 

XS,YS,ZS  coordinates  of  satellite  in 

inertial  space  axis  of  a re- 
lated coordinate  system  which 
lies  along  Z axis 
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APPENDIX  1 


PROGRAM  SATE 


l> 


This  background  version  of  SATVIZE  is  usually  more  valuable  than 
SATVIZE  for  more  than  two  satellites  or  two  ground  stations.  The  high 
speed  printer  allows  the  extensive  output  to  be  more  quickly  and  neatly 
finished  than  the  terminal  output  of  SATVIZE. 


Orbital  elements  are  entered  on  lines  160-240.  The  ground  sta- 
tion latitudes  are  entered  consecutively  for  the  ground  stations, 
then  the  longitudes  (East).  For  example,  the  locations  of  Lincoln, 
Massachusetts  (47.54°  colatitude,  288.73°  E.  Long.)  and  San  Diego, 
California  (57.34°  colatitude,  242.8°  E.  Long.)  are  entered  as: 


DATA  B/47.54,  57.23,  288.73,  242.8/ 


240 
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ePRECiQrrn  rSl-'r  n. 

■•LAADC.NL'T’  jf 


non  o n o o 


r 


//TS042CA  JOB  (6? so , D°1 ,DFSK) , • CHRISTOPHER  P' , NOTI?Y=TSC420,  00000010 
//  TIME=1  00000020 
//  EXEC  FORTGCG  00000030 
//FORT.SYSIN  DD  * 00000040 
C THIS  S A" ELL  IT"  VISIBILITY  PROGRAM  IS  WRITTEN  FOE  A BA  XI  MUM  OF  00000050 
E TEN  SAT"LLIT"S  AND  TWO  OpprjND  STATIONS.  00000060 
C THE  APFAYS  CONTAIN  THE  FOLLOWING  INFORMATION  FOR  EACH  SATELLITE:  000OOO70 
C A — SEMIMAJOB  AXIS  IN  NAUTICAL  MILES  00000C80 
C E — ECCENTRICITY  00000C90 
C W--RIGHT  ASCENSION  IN  REGRESS  00000100 
C WP — ARGUMENT  OF  PERIGEE  IN  DEGREES  000CQ110 
C TP--TI ME  OF  PERIGEE  IN  HOURS  00000120 
C XI--INCLINATION  IN  DEG?  EES  00000130 
C APRAY  E CONTAINS  TH"  COLATITUDF  AND  LONGITUD"  FOR  EACH  GROUND  STATION. 00000 1 40 


DIMENSION  A (10)  ,"(1r'),W<1P),WP(10),TP(1C),XI(10),B(4) 
TATA  A/14  342. ,14342. ,22767. ,22767. , 2276  7.  ,2  27  67., 

3 5716°.  2,  5136°.  2,  5'  36  9. 2,  5736°.  2/ 

DATA  E/. 7 2 5,.  12 5,0.  ,0  . ,0.  ,C.,0.  ,0.,0.,0./ 

DATA  W/C.,  210. ,243.733,243.733,0.  ,0.  , 

4 0. ,0. ,0. ,°./ 

DAT*  HP/- 90.  ,-90.  ,0.  ,0.  ,0.  ,0.  ,0.,0.,n.,0./ 

DATA  TP/n.  , -6. ,0., -6.  ,-16. ,-21., 

1 0.  ,-24. ,-4B. ,-72./ 

DATA  ni/6  3. 435, 63. 435, 23. 4, 23. 4,0.,°., 90. ,9°. ,90. ,90./ 
DA-A  B/153.  436,  26. 565,  27C.  , 0./ 

NUBS  IS  THE  NUMBER  OF  SATELLITES  TO  BE  CONSIDERED. 

NUMG  IS  THE  NUMBER  OF  GROUND  STATIONS  TO  BE  CONSIDFRED. 

NUMS=4 

NUMG=2 

P=3 . ’ 41 5926 


00000150 
00000160 
00000170 
000001 80 
00000190 
00000200 
0000021 C 
00000220 
00000230 
00000240 
00000250 
00000260 
00000210 
00000280 
00000290 
00000300 

c E=  34  40 . 00000310 

RTD  = 57. 2°R7-iori  00000320 

WRAD=15 .0/Ri'D  00000330 

C Cl  IS  THE  V-LOCITY  OF  LIGHT  IN  NAUTICAL  MILTS  PER  SECOND.  000n0340 

CL=1. 61084*10. **5  00000350 

"R  IS  TH"  UPLINK  FREQUENCY  IN  HZ.  00000360 

F®1  IS  THE  DOWNLINK  FREQUENCY  IN  HZ.  00000370 

R?=200.0*10.**6  00000380 

rP1=245.O*10.**6  00000390 

CON=6. 9R7*  1 0 .**  (-6)  00000400 

00  9C0  J= 1 , NUBS  00000410 

4 (J)  = W (J)  /RTD  00000420 

W®  (J)  = WP  (J)  /PTD  00000430 

XI  ( J)  =Xi  ( J)  /RTD  00000440 

TAU  IS  THE  PERIOD  0"  ROTATION  OF  SATELLITE  J . 00000450 

TAU  = CON*A  (J)  **  1 . 5 00000460 

DO  100  K= 1 , NUMG  00000470 

WRI  TE  (6 , 1 0 j 00000480 

10  FORM  AT(  1H  , ' SAT*  , 5 X , ' HP  S • ,12X,' RANGE*  ,9X, 'ANGLE*  ,8  X,  ' UPDOPPLER0C0004  90 

1 • ,PX ,'DNDOPPLSR' ,°X, • SU FL • , 9X , ' SUBC* ,6X , • STATION' ,//)  00000500 

T=8(K)/PTD  00000510 

G = B (K  *2)  /RTD  00000520 

I IS  THE  HOUR.  00000530 

70  *00  1=1  ,25  0000054 0 

?I=T  00000550 

SUBROUTINE  "LLIP  COMPUTES  THE  RANGE  FROM  GEOCENTER  TO  SATELLITE  J AND  00000560 
THE  ANGIE  MEASURED  IN  ORBIT  PL  AN r ; R EFEP RED  TO  NODE  OF  ORBIT  PLANE  00000570 

AND  EQUATORIAL  PLANE.  00000580 

CALL  ELLIP(FI  ,E  (J)  , WP  (J)  ,TP  (J)  , A (J)  ,TAU,FSP,R)  00000590 

E SUP  ROUTT  NT  PRIME  COMPUTES  THE  INF'TIAL  CARTESIAN  COORDINATES  OF  THE  0C000600 

C S AT  ELL  IT  * J (WITH  3EOCENTEP  AS  THE  ORIGIN).  00000610 
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CALL  PR  1ST  (ESP,  W (J)  , XI  ( J)  ,R  , 0 . , X£ , Y S, ZS) 

C ES  IS  THE  SUBSATELLTTF  LONGITUDE  IN  DEGREES. 

”S=  (ATAN  (YS/XS)-WRAD*RI)*RTD 
IF(XS.LT.O.)  FS  = c’S  + P*PTD 
11  IF  (FS.LT.  -360  . ) FS=FS*360. 

IFCS.LT.-Se".  )50  TO  11 

C TS  IS  THE  SU  BSATELL  IT?  COLATITUPF  IN  DEGRESS. 

TS~RTD*AFC9S ( ZS/R) 

F=G  + MR  AD*  R I 

C (X,Y,Z)  ARE  "HE  INERTIAL  CAPTESIAN  COORDINATES  OF  GROUND  STATION  K. 
X = RF*SIN(T)  *COS  (F) 

Y*RF*SIN  (T)  *SIN  (?) 

Z = ?E*COS  (") 

C ACC  IS  THF  ANGLE  PETWFFN  GROUND  STATION  K AND  SATELLITE  J FROM  THE 
C CENT  EF  OF  THE  EADTH. 

C RGF  IS  TH*  PANGS  EPOM  GROUND  S*A”ION  K TO  SATELLITE  J IN  NAUTICAL 
C WILES. 

ACC=APCOS ( (X*XS*Y*YS*Z*ZS)/ (RE*P) ) 

”G?=(  (XS-X)  **2*  (YS-Y)  **2*  (ZS-Z)  **2)  **.5 
ARG=SI N (ACC) *R/RGE 
IF  (AFG. GT. 1 .COOOCCO) A P3* 1 . 0 000" 0"0 
C D IS  THE  ELEVATION  ANGLE  IN  DEGRESS  nFOI!  GROUND  STATION  K TO 
C SATELLITE  J. 


D=  (ARSIN  (A»G)-P/2  . ) *RTD 
RT=SQFT  (R*R-PE*PE) 

IF  (RGF.  IE.  RT)  D=-D 

C SUBROUTINE  DOPE  COMPU~ES  THF  COMPONENTS  OF  SATELLITE  VELOCITY. 

CALL  DOPE  (R  ,E  (J)  , A (J)  ,TAU,ESP,WP  (J)  ,W  (J)  , XI  ( J)  , XV,YV ,ZV) 

C (XT , YT)  ARr  THE  COMPONtnTS  OP  THE  RELATIVE  VELOCITY  OF  GROUND 
C STATION  K . 

XT  = -W R A D*P  E *STN  (T)  *STN  (F)  /1600. 

YT=NPAD*Pr*SIN  (T)  *COS  (F)  /I600. 

DOT*  (X  V-XT) * (YS-X)  ♦ (YV-YT)  * (YS-Y)  *ZV*(ZS-Z) 

DOT*-DOT/RG" 

C UPDOP  IS  THE  UPLINK  ROPPLEP  SHIFT  IN  HZ. 

C INOOP  TS  THE  DOWNLINK  DOPPLER  SH"??  IN  HZ. 

UPDOP=DOT*F?/CI 
DNDOP=DOT*FP1/CL 

WRITE  (6, 20)  J,R  I,  RGF,  D,  U PDOP,  DNDOP,  ES,TS,K 

FOPS  AT  ( 1H  ,T3,EP.1,P17.3,n1U.3,F17.3,E’7.3,F13.3,F13.3,I10) 
CONTINUE 
WP ITS  (6  ,25) 

FORMAT  ( 1 H ,//) 

CONTINUE 
CONTINUE 
WR IIP (6 , 30) 

FORMAT  (1  H ,/////) 

I = NUM  S-1 

IF  (L  . FO  . 0)  GO  TO  3 SO 
T FP  IS  THE  CROSSLINK  FREQUENCY  IN  HZ. 

C Ml  IS  THE  FIRST  SATFLLITE . 

C M2  IS  THF  SFCOND  SATELLITE. 

DO  350  HI =1  ,L 
TAU=CON*A  (H  1)  **1.  5 
N=H1*1 

DO  250  M2=N,NUWS 
WR  TTr  ( 6 , 40) 

POPMAT  (1H  , 2X  ,' HFS  • ,12X, 'RANGE' , 1 1X, 'L12' , 12X,' 02' , 1 1X, 

2 • CFDOPPLFR'  ,6X  , • CVIS • , 6 X , • S AT  1 • , 6X  , 'S  AT 2 • , //) 

TA  U2  = CON*  A (12)  * * 1 . 5 


2C 

500 

25 

"00 
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ir 
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00000870 
00000880 
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00000900 
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00000940 
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00000970 
00000980 
00000990 
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00001050 
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00001120 
00001130 
00001140 
0 CO0 1 1 50 
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00001170 
00001180 
00001190 
00001200 
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DC  300  1=1,05 
PI  = I 

C NOVIZ  DETERMINES  CROSSLINK  VISIBILITY  H HI  PR  0 M FANS  NO  VISIBILITY 
C AND  1 MFANS  VISIBILITY. 

NOVIZ= 1 

CALL  ELLIP(PI,F  (HI)  ,WP  (Ml)  , T P (M  1)  , A (1 1)  ,TAU,FSP,R) 

CALL  EL  LIP  (RI,F(M2)  , WP(M2)  ,TP(K2)  ,A(M2)  , T AO  2,  F2  P,  H2) 

CALL  PRIMF(rSP,W(M1)  ,XT  (HI)  ,R,°.,XS,YS,ZS) 

CALL  PRIME  ( F2 F , W ( M2)  ,X1(M2)  ,P2,0.,X2,Y2,Z2) 

TCI =X2-XS 
TC2  = Y2-YS 
TC3=Z2-ZS 

C PI  2 IS  THF  RANGE  BETWEEN  TWO  SATELLITES  IN  NA'JTICAL  MILES. 

C Fl  2 AND  T 1 2 APE  THE  POINTING  ANGLES  BETWEEN  TWO  SATELLITES  IN  DEGREES 
E12=(TC1**2*TC2**2*TC3**2) «*. c 
FI  2 = PTD*ATA  N (TC2/TC  1 ) 

IF(”C1.LT.r. ) F 1 2=r 1 2*P*  PTD 
VIZ=(-XS*TC1-YS*mC2-ZS*TC3)  / (R12*P) 

I17  (VI2.GT.  .99°99P99)  VIZ  = .9°999999 
DEL  = A P CO c (VIZ) 

DELM=APSIN  ( CF/P) 

IF  (DEL.  LE.DELM)  NOVIZ  = C 
T ’2  = AS  COS  (TC  3/P  1 2)  * P TD 

CALL  DOPE  (P,  E (Ml)  , A (M  1)  , TA  U , FSP  , WP  (Ml  ) , W (Ml ) , XI  (Ml ) , XV  , YV  , ZV ) 
CALL  DOPE  (R2,E  (M2)  , A.  (M2)  ,TA'J2,  F2P,WP(M2)  , W (M2)  ,XI  (M2)  ,S2,02,V2) 
T°T  = (S2-XV)  * TCI  ♦ (U2-YV)  *TC2  + (V2-ZV)  *TC3 
DOT=DOT/R 1 2 

C rOP’2  IS  THE  CROSSLINK  POPPLE"  SHIFT  IN  HZ. 

DOP1  2 = - DOT*  "P/CL 

WPI"’E(6,50)PI,D12,F12,T12,DOP12,NOVIZ,M1,M2 
c°  FORMAT  (1H  , '='5.  1 , F 17  . 3 , F 1 4 . 3 , El  5 . 3 , F 2D  . 1,  19 , 1 1 0,  II  0) 

300  CON  TIN  OF 

W R ITS  ( 6 , 6 °) 

60  FORMAT  ( 1 H ,//) 

250  CONTINUE 

350  CON'INl'F 

END 

SUBPOUT  INF  PRIME  (FSP  , WS  , XIS  , X PS  , Y PS  , X S,  YS , Z S) 

A1 1 = COS  (ESP)  *COS  (WS)  - COS  (XIS)  *SIN  (WS)  *SIN  (FSP) 

A12  = -SIN(FSP)  »CCS  (WS)  -COS  (XIS)  *SIN  (WS)  *COS  (FSP) 

A2 1 =COS  (FSP)  *SI  N (WS)  *COS  (X  IS  ) *COS  (W  S)  *S  IN  ( FSP) 

A22=-STN  (FSP)  * SIN  (WS)  *CO  S (X  IS)  *COS  (WS)  *COS  (FSP) 

>31  = STN(XIS)  * S I N ( F S P) 

A 32  =SI  N (XIS)*  COS  (FSP) 

XS  = A11*XPS*A12*YPS 
YS=A21*XPS»A22*YPS 
ZS  = A.  31  *XPS*A32*  YPS 
p^URN 
END 

SUBROUTIN  F "LLTP(T,,',NP,TP,A,TAU,FSP,  ") 

P=3.  1 41 5°  26 
7 = 2 . * P*  (T-TP)  /'"An 
P2=2.*P 

2 IF(Z.GT.P2)  Z = Z-P2 

IF  (Z.GT.  P2)  GO  TO  2 
El  =7 *E* SIN  (Z) 

F2=  (Z*F*(FIN  (FI)  ) - (F*COS(F1)  ) *F  1)  / (I  ,-F*COS  (El)  ) 

0*o. 

u r‘1  = (Z  ♦ F*  (SIN  (r2)  ) - (r*C0S  ("  2 ) ) • = 2 ) / ( 1.  - E*COS  ( E2)  ) 

0=0*1. 

DF=F3-E2 


C 0001 2 1 0 
00001220 
00001230 
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DE2=DE**2 

E2=F3 

IF  (DE2 . GT.  • 00000001)  GO  TO  4 
TH=»PCOS  ( (COS  (E2)  -E)  / (1 . - E*COS  (E2)  ) ) 
IF(Z.GT.P)?H*2.*P-TH 
FSP=WP*TH 

F=  (A*  (1  .-F»*2)  )/(1.  ***COS  (FSP-WP)  ) 

PFTOPN 

END 

SOBPOOTINE  DOPE(P,F,A,TAO,FSP,ilP,IIS,XIS,-D,  TD,ZD) 

P=3.  1415926 

T=FSP-HP 

FD=  (A*E*2.*P/(TAU*  (1.-E**2)**.5)  )*SIN(T) 

TD=  (2.  *P/TA0)  *(  (1.-E**2)  **(-1.5)  ) * (1  .♦E*COS  (T)  ) ** 2 
XPC=FD*COS  (FSP)  -P*TD*SIN(FS  P) 

XPC=XPC/3600. 

YPC=PD*SIN (FSP)  *P*TD*COS(FSP) 

YPC=YPC/?609. 

CALL  PP  IN  E (9  . ,WS,XIS,XPC,YPC,XD,Yr,ZD) 

PETOPN 

END 
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APPENDIX  2 


S’ 


PROGRAM  SATD 


Unlike  SATE,  SATD  runs  in  double  precision.  It  is  rarely  used 
except  to  check  crosslink  Doppler  for  low-eccentricity  satellites 
within  a given  constellation.  Differences  of  large  velocities  can 
be  found  accurately  with  SATD,  and  accurate  crosslink  Dopplers  result. 
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//TS042  o A JOB  (6’60  ,D°1  .DESK)  , 'CHPISTOPHER  P' , NOTIFT=TS04  20  , 00000010 

//  TTME=3  00000020 

//  EX  EC  EORTGCG  00000030 

//PORT . SYS!  N DD  * 00000040 

C THIS  SATELLITE  VISIBILITY  PROGRAM  IS  WITTEN  FOR  A MAXIMUM  OP  00000050 

C TEN  SATELLITES  AND  TWO  GPCUND  STATIONS.  00000060 

C THE  ARRAYS  CONTAIN  THF  FOLLOW  I NG  INFORM  AT  ION  POR  EACH  SATELLITE:  00000070 

C A — SFHIMAJOR  AXIS  IN  NAOTICAL  MILES  00000080 

C E — ECCENTRICITY  00000090 

C W — FIGHT  ASCENSION  IN  DEGREES  00000100 

C WP — AFGOMENT  Or  PERIGEE  IN  DEGREES  C0000110 

C TP — TIME  OF  PERIGEE  IN  HOURS  00000120 

C XI--INCLINATION  IN  DFGPEES  0000013C 

C ARRAY  E CONTAINS  THE  COLATITUDE  AND  LONGITUDE  POR  EACH  GROUND  STA TION . 000001 40 
IMPL ICTT  F*AL»8  (A-H.O-Z)  00000145 

DIMENSION  A (1C)  ,E  (10)  ,N  (10)  , Wp  (10)  ,?P  (10)  ,XI  (10)  ,B  (4)  00000150 

DATA  A/22’67. ,22767. ,22767. ,22767., 36140. 36,36140.36,  00000160 

3 36140. 36, 36 140. 36, ’4342. 3,14342 .3/  00000170 

DATA  E/.  1 ,.  1 ,.  1 ,.  1 ,.  1 ,.  1,.  1, . 1,  . 65,  .65/  00000180 

DATA  W/243. 733, 24 3. 733, 24 3. 733, 24 3. 73 3, 243. 733, 243. 733,  00000190 

4 243.733,  243.733,243.733,243.733/  00000200 

DATA  WP/C. ,7. ,P.  ,0.  ,C. ,0. ,0. ,0. , -90. 0,-90. n/  00000210 

DATA  TP/0. ,-6 . ,- 12 . ,- IP . ,0. ,-12. ,-24. ,-36. ,-.75,-6.75/  00000220 

DATA  XI/23. 4,23.4,23.4,23.4,113.4,113.4,113.4,113.4,63.4,63.4/00000230 
DATA  B/4  7 . su,S7.23,2PR.7’,242.  8/  00000240 

C NUMS  IS  the  NUMBER  OP  SATELLITES  TO  BE  CONSIDERED.  0000025C 

C NUHG  IS  THr  NUMBER  OP  GROUND  STATIONS  TO  BE  CONSIDERED.  0000C260 

NUMS-17  00000270 

NUNG=  2 00000280 

P=3.  1415926  00000290 

op=3U40 . 00000300 

RTD=57. 24577061  000C03 1 0 

WR AD= 15. 0/RTD  00000320 

C CL  IS  THE  VELOCITY  Of  LIGHT  IN  NAUTICAL  MILES  PER  SECOND.  00000330 

CL=1.61°84*1C.**5  00000340 

C FR  IS  the  UPLINK  FREQUENCY  IN  HZ.  00000350 

C PP 1 IS  THE  DOWNLINK  FREQUENCY  IN  HZ.  CC000360 

FR  = 300 . 0» 1 0 .»*  6 00000370 

'■R1  = 24C.0*1D.**6  00000380 

CON=6 . Q87* 1 0 .«» (-6)  00000390 

DO  900  J = 1 , NUN  S 00000400 

W(J)  =W  (J)  /PTD  00000410 

WP(J)=WP(J)/PTD  O00C042  0 

XI  (J)  =XI(J)  /FTD  30000430 

P TAU  IS  TH*  PERIOD  OF  ROTATION  OF  SATrLLITE  J.  00000U40 

TAU=CON«a  (J)«*l . 5 00000450 

DO  700  K = 1 , NUM G 00000460 

WRITE  (6,10)  00000470 

10  FORMAT ( 1 H , 'SAT1 ,5X , 'HRS ' , 1 2X , ' RANGE' , 9X , • ANGLE* , 8X , ' UPDOPPLER00000480 

1 ' , 8X,' DN DOPPLER • , °X, • SUBL* ,9X,' SUBC • ,6X, 'STATION' ,//)  00000490 

7=  B (K)  /RTD  00000500 

G = B (K*2  ) /PTD  00000510 

C I IS  THF  HOUR.  00000520 

DD  500  1=1,  25  00000530 

PI=I  00000540 

C SUBROUTINE  ELLIP  COMPUTES  T Hr  PA NGE  FROM  GEOCFNTER  TO  SATELLITE  J AND  00000550 

C THE  ANGLE  MEASURED  IN  ORBIT  PLA NE ; REFERRED  TO  NODE  OP  ORBIT  PLANE  00000560 

C AND  EQOATQPTAL  PLANE.  00000570 

CALL  ELL  It  ,RI,E(J)  ,WP(J)  ,TP(J),A(J)  ,TAU,FSP,R)  00000580 

C SUBROUTINE  PPIME  COMPUTES  tHe  INERTIAL  CARTESIAN  COORDINATES  OF  THE  00000590 
C S AT  ELI  IT r J (WITH  GEOCENTER  AS  THE  ORIGIN).  00000600 


t 


CALL  PPIME(E*P,  W(J)  ,XI  (J)  ,R  ,0.,  XS,YS,ZS) 

C FS  IS  THE  SUBSATELLITF  LONGITUDE  IN  DEGREES. 

*S  = (DATAN  (YS/XS)  -NRAD*RI)  *RTD 
IFfXS.LT.O.)  FS=FS*P*PTD 

C TS  IS  THE  SUBSATELLITE  COLATITODF  IN  DEGREES. 

TS=RTD*  DA  RCCS  (ZS/P) 

P=  G*WR  AD*R  I 

C (X,Y,Z)  APE  THE  INERTIAL  CARTESIAN  COORDINATES  OF  GROUND  STATION  K. 
X=RE*DSIN  (T)  *DCOS  (F) 

Y = P E*  DS  IN  (T)  * DSTN  (F) 

Z=PE*DCOS(T) 

C ACC  IS  THE  ANGLE  BETWEEN  GROUND  STATION  K AND  SATELLITE  J FROH  THE 
C CFNTEP  OF  THE  EARTH. 

C RGE  IS  THE  RANGE  FROM  GROUND  STATION  K TO  SATELLITE  J IN  NAUTICAL 
C RILES. 

ACC=DARC05  ( ( X* XS*Y* YS *Z* ZS)  /(PE*R) ) 

RGE  = ( (XS-X)  **?*  (YS-Y)  **2*  (ZS-Z)  **2)  *•  .5 
ARG  = DSI N (ACC) * R/RGE 
IF(ARG.GT. 1.0000000)ARG=1. OOOCOCOO 
C D IS  THE  ELEVATION  ANGLE  IN  DEGREES  FPOH  GROUND  STATION  K TO 
C SATELLITE  J. 

D=  ( PAR  S IN  <ARG)-P/2.)*PTD 
RT  = DSQRT(R*7-RE*RE) 

IF  (PGF.  LE.  RT)  D=-D 

C SUBROUTINE  DOPF  COMPUTS  7HF  CONPONENTS  OF  SATELLITE  VELOCITY. 

CALL  DOPE(R,E(J)  , A ( J)  ,TAU,ESP,WP  (J)  ,W  (J)  , XI  (J)  , XV,YV,ZV) 

C (XT , YT)  ARE  THE  CONPONENTS  OF  THE  RELATIVE  VELOCITY  OF  GROUND 
C STATION  K . 

XT  = -WRAB*P”*RSTN(T)  *DSIN  (*)  /3600  . 

YT=WPAD*RE*DSIN(T)*DCOS  (E)/3600. 

DOT=(XV-X^)  * (Yc-X)  ♦ (TV-  YT)*  (YS-Y)  *ZV*  (ZS-Z) 

DOT=-DOT/RGF 

C IJPDOP  IS  THE  UPLINE  DOPPLEP  SHIFT  IN  HZ. 

C DNDOP  IS  THE  DOWNLINK  DOPPLER  SHEET  IN  HZ. 

UPDOP=DOT*FR/CL 

DNDOP=DOT*FR1/CL 

WRITE  (6, 20)  J,PI,RGE,D,UPDOP,DNDOP,FS,TS,K 
20  FORMAT  ( EH  , I 3,  F<>.  1 , F17 . 3,  FI  9.  3,  T1  7.  3 , *17 . 3,  FI  3.  3 ,F1  3 . 3 , 1 10) 

500  CONTINUE 

WRITE  (6, 25) 

25  *ORMAT(1H  ,//) 

700  CONTINUE 

900  CONTINU* 

WRITE<6, 30) 

30  FORMAT  (1H  ,/////) 

L = NMN  S- 1 

IF  (L.  FQ.  0)  GO  TO  350 
C FP  IS  TH*  CROSSLINK  EPEQUFNCY  IN  HZ. 

ER  = 2. 83*10.**  13 
C Ml  IS  THE  FIRST  S A T ELL I TF . 

C M2  IS  THE  SECOND  SATELLITE. 

DO  350  M1  = 1,L 
TAD=CON»A  (Ml)  **1.5 
N=  HI  ♦ 1 

DO  250  H 2 = N , NON  S 
WRITP(6,9P) 

96  FORMAT  (1  H , 2 X , • HRS  ' , 1 2X  , • R ANG*«  , 1 1X  , ■ L 1 2*  , 1 2X  , • C 1 2'  . 1 1X  , 

2 ' CR  DOPPLER'  , 6 X , ' C VI  «•  ,6  X,  • SAT1  • . 6X  , • S AT2  • ,//) 

T AU  2xCON  *A  (M2)  **1.5 
DO  300  1=1,25 
PI=I 
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NOVIZ  DSTEPBTNSS  CROSSLINK  VISIFILTTY  WHFPE  0 BEAMS  NO  VISIBILITY 
AND  1 BEANS  VISIBILITY. 

NOVI Z=1 

CALL  ELLIP  (rI  (B1)  ,VP(A1)  ,TP(B1)  , A (B  1)  .TAU.FS",  P) 

CALL  ELLIP  (PT,F  (H2)  ,HP  (B2)  , TP(B2)  ,A(B2),TAU2,F2P,P2) 

CALL  PR  IMF  (ESP,  W(H  1)  , XI  (B  1)  , R , 0 . , XS  , Y S , Z S) 

CALL  PPIHE  (F2P,  W(B2)  ,XI  (N2)  , R2,  0.  ,X  2,  Y2,Z2) 

TC1=X2-XS 
TC2  = Y 2-YS 
TC3=Z2-ZS 

C PI  2 IS  THE  RANGE  BETWEEN  TWO  SATELLITES  IN  NAOTICAL  BILES. 

C F 1 2 AND  T 1 2 APE  THE  POINTING  ANGLES  BETWEEN  TWO  SATELLITES  IN  DEGREES 
R12=(TC1**2*TC2**2  + ?C3**2)**.5 
FI  2=RTD*DATAN  (TC2/TC1 ) 

IF(TC1.LT.1.)  E12=F12*P*PTD 
VIZ=(-XS*^C1-YS»TC2-ZS*TC3)/(P12»P) 

IF(VIZ.GT..9q99R9RO)VIZ  = .9°9  9c>999 
DEL  = DAP  COS  (VIZ) 

DELB=DAFSIN  (FE/P) 

IF  (DEL.  LE.  DELB)  NOVIZ=0 
T12=DAPCOS  (TC3/P12)*PTD 

CALL  DOPE  (P,E  (B  i)  , A (B1)  , TAU,  FSP,  W?  (HI ) , W (B1)  , XI  (B1)  ,XVrYV  ,ZV) 
CALL  DOPF  (P2  , " (12)  , A (12)  ,TAU2,  ”2P,WP(N2)  , W ( IB 2)  , XI  (H2)  , S2.U2,  V2) 
DOT=  (S2-XV)*rT'C1  + (U2-YV)  *TC2*  ( V 2 — Z V ) *T  C 3 
DOT  = DOT/P  1 2 

C DOP 1 2 IS  THE  CROSSLINK  DOPPLEP  SHIFT  IN  HZ. 

DO  PI  2 = - DOT*  FP/CL 

WFITE(S,SO)HT,P12,E12,T12,DOP12,NOVIZ,H1,H2 
SC  F0PEAT(1H  ,r5.  1 ,Er>.  3,F1U.  3,  F15.  3.F20.3, 19,11  C,  I1F) 

300  CONTINUE 

WP  ITE  ( 6 , Sn) 

PC  FORBAT(1H  ,//) 

250  CONTINUE 

350  CONTINUE 

END 

SUBFOUT  IN  r PPIBE(FSP,WS,XTS,XPS,YPS,XS,YS,ZS) 

IHPLICIT  REAL*P  (A-H,0-Z) 

All  =DCOS(Fsr)  *DCOS  (WS)  - DC  OS  (X  IS)  *PSIN  (WS)  *DSIN  (FSP) 

M2  = -DSTN  (FSP)*DCOS  (WS)-DCOS(XIS)*DSIN(WS)  *DCOS(FSP) 

A21  = ECOS(  FSP)  *DSIN  (WS)  *DCOS  (XIS)  »DCOS  (WS)  *DSI  N (FS  P) 

A22=-DSIN  (FSP)  *DSIN(WS)  »DCOS  (XIS)  *DCOS  (WS)  *DCOS  (FSP) 

A31  = TSIN  (XIS)*DSIN  (FSP) 

A3 2=  DS IN  (XIS)  * DCOS  (ESP) 

X S=  A 1 1*XPS* A 1 2*YPS 
YS  = A.  21*VPS*A22*YPS 
ZS=A31»XPS*A32*YPS 
PFTUPN 
ENT 

SUBFO'JT INF  FLLIP  (T,t',WP,TP,A,TAU,FSP,R) 

IBPLICI"  PrAL*P  (A-H.O-Z) 

P=  3 . 1 9 1 5°26 
Z=2.*P* (T-TP) /TAU 
SSM  = DSIN(5.  *Z) 

S6H=DSI N (6  . * Z ) 

STH  = DSIN  C.»Z) 

C5N  = DCOS(5.»Z) 

C6H=DCOS  (6.*Z) 

C^H  = DCOS  C7 . * Z) 

SH=DSI N ( Z) 

S2B  = DSIN(2.*Z) 

S 3H  = DS  IN  ( 3.*Z) 


0C001220 
C 0001 23  0 
00001290 
00001250 
00001260 
00001270 
00001280 
00001290 
00001300 
00001310 
C0001  320 
00001330 
00001390 
00001350 
00001360 
00001370 
00001380 
00001  390 
00001900 
0000191  0 
00001920 
00001930 
00001990 
0C00195C 
00001960 
0.0001970 
00001980 
00001490 
00001500 
000C1510 
00001520 
00001530 
00001 590 
0C0C1550 
00001560 
0F001 570 
00001575 
00001580 
00001590 
0000160C 

000  0 1 6 1C 
00001620 
00001630 
00001690 
00001650 
00001 660 
00001670 

00001  680 
00001  690 
00001695 
00001700 
00001710 
0000 17  20 
00001730 
00001790 

0000  1750 
00001760 
00001710 
00001780 
00001790 

00001  800 
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00001810 
00001820 
00001  8j  0 

00001840 
00001850 
00001860 
00001 820 
00001880 
00001 890 
00001900 

00001910  y 

000  01 92  C 
00001925 
00001930 

00001  940 
00001950 
00001960 
00001970 
00001980 
0000 1990 
00002000 
00002010 
0000202C 
CC902030 


S4B  = DSIN(4.»Z) 

CB=DCPS  (2) 

PSP*Z*2.*E*SB+1 . 25*  (F**2)  *S2B*  ( (E»*3)  /12.)  *(13.*S3B-3.*SB) 
F?P=FSP*WP*  ( (E**4)  /96.)  *(103.»S4B-9U.*S2B1 
f'=  ( (■=■-•» 5)  /960. ) ‘(loop.  *s  59-645.  *S3b*50.*SB) 

P6-  ( (»*»6)  /960.  )*(1223.*S6B-902.*S4B*85.*S2B) 
r7*  ( (E**7)  /32256.)  * (47273. *S7b-41699.*S5B*5985.*S3B*74°.*CB) 
FSP=PSP*F5*F6*F"' 

P«  (A*  (1  . -E**2)  )/  (1.  ♦?*DCOS  (FSP-WP)  ) 

PFTOPN 

END 

SOBFOUTIN  F DOPE  (P  ,F,A  ,TA0  , FS  P,  4 P,  US , X IS  , X D,  Y D,  ZD) 

IB  PLTCIT  PEAL*8  (A-H.O-Z) 

P=3. 1415926 
T=FSP-WP 

PD=  (A*E*2.  *P/(TAO*(  1 .-E**2)  **.t-))*DSTH(T) 

TD  = (2.*P/TAU)  *(  ( 1.  -E*»2)  *»(-1.5))*(1.  ♦E*DCOS  (T)  ) * * 2 
XPC  = BD*DCOS  (FS  P)  -P*TD*  DSIB  (FSP) 

XPC=XPC/3600. 

TPC=PD*DSIN  (FS  P)  ♦P*TD*DCOS(FSP) 

YPC=YPC/3600. 

CALL  PRIRF(P.,8S,XI<!,XPC,  YPC  , XD  , YD  , ZD) 

PT TURP 
END 


APPENDIX  3 


s' 


PROGRAM  SATVIZE 


SATVIZE  is  an  interactive  FORTRAN  program  which  queries  the  TSO 
user  for  the  number  of  satellites,  number  of  ground  stations,  orbi- 
tal elements  of  each  satellite,  location  of  each  ground  station,  and 
frequencies.  Default  values  are  included  in  the  program  to  allow 
the  undecided  user  a chance  to  study  the  form  of  the  output. 

Arbitrary  Keplerian  elements  can  be  entered,  with  0 < eccen- 
tricity < 0.99. 


95 


IT 

115 

120 

’21 
1U2 
1 1r 

145 

1 15 
1=5 

150 

Of 

1 40 
605 


DIMENSION  1(H)  ,"(11)  , « (11)  ,«P  (1 1)  , IP  (1 1)  ,XI  (11)  ,C(3)  ,RL(3)  0 0000010 

J1=C  0 0000020 

DATA  A/0. ,22767. ,22767.  ,227  67.,  22767.,  36140. 36,  36140. 36,  0 000003  0 

1 36140.  36,  36140.36,  14342.3,14  342.3/  00000040 

DATA  E/0.,. 1, . 1,. 1, . 1,. 1, .1, .1, .1 ,.1 ,.1/  00000050 

DATA  N/24  7. 133,24  3. 13  3, 24  3.  7 33,  24  3.  7 3 3,  24  3.  73  3,  24  3. 7 33,  243.  7330  0000060 
2,  24  3. 73  3,  243.  733,  24  3. 733,  243. 73  3/  00000070 

DA-A  WP/O.  , 0.  ,0.  ,0.  ,0.,  0.  ,C.,C.  ,0.  ,-90.  '',-90.0/  00000080 

DATA  TP/C.  ,i . ,-6.  ,-12.  ,-18.  , 0. , -12.  ,-24.  ,-36. ,-  .75,-6.75/  00000090 

DATA  XT/2  3.  4,  23. 4, 2’. 4, 23. 4, 23. 4, 1 1 3 . 4 , 1 1 3 . 4 , 1 1 3 . 4 , 1 1 3.  4 , 63.  4 , OOOC0 1 00 
163.4/  00000110 

DATA  C /O. ,47. 54,57. 23/  000C0120 

DATA  HL/O., 288.73, 242. e/  C00C0130 

WRITE (6, 1 10)  00000140 

FORMAT ( 1 H ,///)  OnOOQ1 50 

WRITE  ( 6 , 115)  00000160 

FOP MA T ( 1 H , ’"HIS  S ATFLL ITr  VISIBILITY  PROGRAM  IS  WFITTEN  FOP  A"0000170 

I MAXI  MOM  Or  T SA^ELLITR'S  AND  2 GPOTJND  STATIONS.',//)  00000180 

WRITR(6, 120)  00000190 

FORMAT  (1H  , 'FOP  EACH  SATELLITE  THE  TLLOWING  DATA  MOST  BE  SUPP09000200 
’LIED:',/,'  SPMI  MAJnP  AXIS',/,'  SCCEMT  RICITY  • , /, ' RIGHT  A SCENSION  ' , C 00002 1 0 
V,'  ARGUMENT  0"  PERIGEE', /,'  TIME  O'  PEFIGT  • , /,  ' I N CL  IN  ATION  • , //)  0 0000  22  0 
WPI”F(6,  125)  00000230 

Ri.R"f?(1H  ,"’0?  Each  GPO'JNr  STATION  THE  FOLLOWING  DATA  MOST  B00900240 

II  SDPDLIED:'  ,/,  ' COLATITUDE'  ,/,  ' LONGITUDE',//)  C0000250 

»?r'(6,14i)  00000260 

FORMAT (1H  , 'ADDITIONAL  DATA  WHICH  MUST  PE  SUPPLIED:',/,'  UPLIN0000027C 

If  FREQUiNCY' ,/,  • DOWNLINK  FREQUENCY ',/,  • CROSSLINK  FREQUENCY • , //)  00000280 
WP  ITE  ( 6,  111)  00000290 

FORMAT  ( 1 H , ••’’H'7  USrR  IS  ASKED  TO  ANSWP  SEVERAL  QUESTIONS.  •,/, 00000300 
1*  Ir  THE  "5  FD  DOFS  NOT  KNOW  THE  VALUE  OR  AN  IN^U?  V APIA BL ",  • , /, • H0CO0031  n 
IE  SHOULD  EN " RR  THR  NUM0RR  400.0.  A Dp  R A ULT  VALUE  WILL  THEN  BE  US ER 0000 320 
1".  ’,/,'  THE  USEP  IS  RXPECTED,  HOWEV  R®,TO  SPECIFY  THE  NUMBER  OF  SATE00000330 
1LLITRS  ANT  GFOUND  STATIC'S  TO  3r  CONSIDERED.  ',////)  0000034C 

WFITE(6,  ia')  00000350 

FOP»AT(1H  , ' HOW  MANY  S A"  rLL  ITES  APE  TO  BE  CONSIDERED?  ENTER  INO"000  360 


IN  COLUMN 


') 


00000370 
00000380 
00000390 
00000400 
00000410 
NT00000420 
00000430 
00000440 
00000450 
n0000460 
00000470 
00000480 
00000490 

EN"3F  FOLLOWING  DA"A  IN  DECIMAL  FCRM  ANYWHEPE  IN  TO00005O0 


1 TFGER  WITH  RIGHTMOST  DTGI 
READ  150,  M 
WPI  TE  (6  ,1  35 ) N 
FOFM  AT ( 1 H ,13) 

WRITE  (6,  1 e 5) 

cOPM A T ( 1 H ,'HOW  MANY  GFOUND  STATIONS  ARr  ”0  BE  CONSIDERED? 

1 EP  INTEGER  IN  COLUMN  2.') 

READ  150,  N 
WPI  TE  (6,135)  N 
FOPMAT  (12) 

W°I"R  (6,  25) 

WPI  TE  (6 ,25) 

WRITE  (6,0R) 

FOFM  AT ( 1 H , 

1 HE  FIRST  2n  CnLUMNS  . ' ,//)  00000510 

WRI"R  (6 ,1  40)  00000520 

FOFM  AT  ( 1 H ,'WHAT  IS  THE  SEMINAJOP  A XI  F VALUE  IN  NAUTICAL  MI  L ES000  00  53C 

1 'OF  SATELLITE  1?')  00000540 

DO  101  1=1, M 00000550 

K=  T* 1 00000560 

R'AD  170,  A(I)  00000570 

IF  (A  (I)  . NE.  400. 7)  GO  TO  605  00000580 

A ( I)  = A ( K)  0000D590 

.11  = 1 00000600 

WR  ITE  ( 6 , 30 5)  A ( I)  00000610 
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305  F0RNAT(1H  ,F20.8) 

IF  (I.  EQ.  M)  GO  TO  1C  1 
PRINT  3 10  , K 

310  pOP!HT  ( 1 H , 'FOR  SATELLITE’, 

101  CONTINUE 

WR  ITF  (6 , 3 1 5) 

315  FORMAT ( 1H  , /, • WHAT  IS  THE 

DO  103  1=1, M 
K=  1 + 1 

DEAD  170,  F ( I) 

IF  (F  (I)  . NF.  «n0.  ")  GO  TO  61" 
1(1)  =MK) 

J1=1 

610  RRITF  (6 ,305  ) F (I) 

IF  (I.  SQ.M)  GO  TO  103 
PRI  NT  310,  K 
1n3  CONTINUE 

WRITE(6, 320) 

320  FORMAT ( 1H  , /, • WHAT  IS  THE 

1 R SATELLITE  1?'  ) 

DO  10a  1=1, M 
K=  I ♦ 1 

READ  I^O,  W(I) 

IF(W(I)  .NE.  UOO.OJGO  TO  615 
H(I)=W(K) 

J1  = 1 

6 IF  WRIT"(6,  30F)  w (I) 

IF(I.',Q.H)GO  TO  1C4 
PRINT  310,  K 
10  4 CONTINUE 

WRITE  (6  ,325) 

325  PORMA’-JIH  ,/,'  WHAT  IS  TH71 

’S  NOR  SATELLITE  1?') 

DO  106  1=1, M 
K=  I ♦ 1 

READ  INC,  WP(I) 

IF  (WP(I)  .NE.UOO.O)  GO  TO  620 
WP  (I)  =WP  (K) 

J1=1 

620  WRITE(6,  305)  WP(I) 

IF  (I.FQ.  M)  GO  TO  106 
PRINT  310,  K 
106  CONTI  NUE 

WRITE (6,330) 

330  FORM  AT  ( 1H  , /,  • WHAT  IS  THE 

1ITN  1?') 

DO  107  1 = 1 , M 
K = !♦  1 

READ  170,  TP  ( I) 

IF  (TP  (I)  . NF.  400.0)  GO  TO  625 
TP  (I)  =TP  (K) 

J1  = 1 

625  WRITE  (6,305)TP(I) 

IP(I.  EQ.M)  GO  TO  107 
PRINT  310,  K 
1 C 7 CONTINUE 

WRITE(6,  335) 

335  PORMAT ( 1H  , /, ' WHAT  IS  THE 

1 TELLITE  17'  ) 

DO  ICR  T=1,M 
K«I*  1 


000C062C 

CC00063C 

00000640 

13,  * ?' ) 000O0650 

00000660 

00000670 

ECCENTRICITY  OP  SATELLITE  1?')  00000680 

00000690 

00000700 

00000710 

00000720 

00000130 

00000740 

00000750 

00000760 

00000770 

00000780 

00000790 

RIGHT  ASCENSION  VALUE  IN  DEGREES  F000000800 

00000810 

00000820 

00000830 

00000840 

00000850 

00000860 

00000870 

00000880 

00000890 

0C0C0900 

00000910 

00000920 

ARGUMENT  OF  PERIGEE  VALUE  IN  DEGRFE0C0C0930 

00000040 
00000950 
000O0960 
00000970 
00000980 
00000990 
00001  000 
00001010 
000C1C20 
00001  030 
00001040 
00001C50 

TIHF  OF  PERIGE"  IN  HOURS  FOR  SATELL00001 060 

0000 1070 

0 CO  01 080 

00001  090 
00001100 
000C1110 
00C0112C 
00001130 
00001140 
00001150 
00001160 
00001170 
00001180 

INCLINATION  VALUE  IN  DEGREES  EOF  SA0000H90 

000C1200 

00001210 

00001220 
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READ  in,  X I ( X) 

IF(XI(I)  .NE. 400.0)00  TO  6 30 
X I ( X)  =X  T (K) 

J1  = 1 

630  WRITE  (6 ,305)  XI  (I) 

IF (I. FQ.M) GO  TO  ice 
PRINT  310,  K 
118  CONTINUE 

WRITE  (6,  341) 

3Ui  FORMA  T (1  H ,//, 

1 TATTON  1?') 

DO  102  1=1, N 
K = I + 1 

READ  HO,  C(T) 

I F ( C ( I)  .NE.  400  .)  GO  TO  635 
C (I)  =C  (K) 

J1=1 


000C1230 
00001  280 
00001250 
00001260 
00001270 
00001280 
00001290 
00001300 
00001310 

WHAT  IS  THE  COLATITUDF  IN  DEGREES  FOR  GROUND  S00001320 

00001330 
00001390 
00001350 
00001360 
00001370 
0C0C1380 
10001  390 


635 


20  5 
102 

345 


WRITE  (6, 305)  C (I) 

IF  (I.  EQ.  N)  GO  TO  1A2 
PRINT  205  , K 

FORM  AT ( 1 H ,'FOP  GROUND  STATION ' ,1 3 ,'?' ) 
CONTINUE 
WP  ITF  (6,345) 


00001400 
00001410 
00001  420 
00001430 
00001440 
00001450 


FORMAT  ( 1H  ,/,'  WHAT  IS  THE  LONGITUDE  IN  DEGREES  FOR  GROUND  STA000C1460 
1TION  1?')  00001470 

DO  109  1=1, N 00001480 

K= I ♦ 1 00001490 

READ  171,  Pl(T)  00001500 

171  FORMAT  (F20. 8)  00001510 

IF(RL  (I)  .NE.40C.C)  GO  TO  641  00001520 

RL(T)=?1(K)  00001530 

J1=1  00001540 

641  WRITF  (6, 305)  PL  (I)  00001550 

IF  (I.EQ.  N)GO  TO  10°  00001560 

PRINT  205,  K 00001570 

10°  CONTINUE  00001580 

WRITE  (6  ,841)  00001590 

840  FORMAT(1H  ,////,'  ENTER  THE  FOLLOWING  DATA  IN  EXPONENTIAL  FQRM00001600 

2 (EG.  300 .0*10  . **6  IS  WRITTEN  AS  300  . OS  6)  • , / , • WITH  THE  RIGHTMOST  0000  1610 
2DIGIT  APPEARING  IN  COLUMN  10.',//)  00001620 

WP  ITE  (6  ,845)  00001630 

845  FORM  AT ( 1 H ,'WHAT  IS  iHE  UPLINK  FREQUENCY  IN  HZ?')  00001640 

PFAD  850,  F P 1 00001650 

981  FORMAT (E10. 2)  00001660 

IF (FR 1 . =Q . 400. C) Fc 1 =300. 056  00001670 

WPITE(6,855)  FP1  00001680 

855  FORMA"(1H  ,E20.C)  00001690 

WPI”E(6, 860)  00001700 

960  FORMAT ( 1H  , /,  ' WHAT  IS  THE  DOWNLINK  FREQUENCY  IN  HZ?')  00001710 

READ  861,  F*2  00001720 

IF  (FR2.  EQ  . 400. 0)  FR  2 = 2 4 5.  OS  6 00001730 

WFI”= (6, 855) FR2  00001740 

WRITE  (6  ,855  ) 00001750 

865  FORM  AT  ( 1 H , /,  ' WHAT  IS  TH=  CROSSLINK  FREQUENCY  IN  HZ?')  00001760 

PEAD  850,  FR3  00001770 

IF (FR3.FQ.410.l)FP3=38.OE9  00001780 

WRI"E(6, 855) FP 3 00001790 

IE  ( J1 . EQ.O)  GO  TO  660  00001800 

WRITE  (6  ,6i0)  00001810 


670 


FORM  AT ( 1 H ,//,'  YOU  HAVE  NOT  USED  A COMPLETE  DATA  SET.',/,'  TH00001  820 


1 EPEPORF , "HE  RESULTS  APE  FOP  DEMONSTRATION  PURPOSES  ONLY.') 


00001830 


98 


ff 


660 

665 

675 

1"5 


190 


GO  TO  6 i 5 00001840 

HFITE(6  ,665)  00001850 

FORM  AT ( 1 H ,//,'  YOUR  DATA  SET  IS  COMPLETE.')  00001860 

MR!  TE  (6 , 105)  00001870 

FORMAT  (1  H ,///)  00001880 

P=3.  1415926  00001890 

p==3440.  00001900 

PTD=57 .29577951  00001910 

HRAD=  15.0/PTP  00001920 

CL=1 . 61984»10. **5  00001930 

CON=6. 087*10.**  (-6)  00001940 

HRITE(6,199)  00001950 

FORM»T(1H  .'P^SULTS  LISTED  IN  ORDER  0"  APPEARANCE  ARE:',/,'  RA00001960 

2NGE  BET8EEN  GROUND  STATION  AND  SATELLITE  IN  NAUTICAL  MILES',/,'  ELC000197C 

2EVATION  ANGLF  IN  DEGREES  FROM  GROUND  STATION  TO  SATELLITE',/,'  DPL00001980 
2INK  DOPPLER  SHIFT  IN  HZ',/,'  DOHNLINK  DOPPLFR  SHIFT  IN  HZ',/,'  SUPC0001990 
2 S A T ELL  ITE  COLATITUDE  IN  DEGREES',/,'  SUBSATFLLITF  LONGITUDE  IN  DEG00002000 


2REFS' , ///) 

DO  890  I N=  1 , M 
H (IN)  =W  (IN)  /PTD 
HP  (IN)  =WP  ( IN)  /RTD 
XI  (IN)  = XT  (IN)  /FTP 
890  CONTINUE 

TO  900  J= 1 , M 
TA  U=CON*A  (J)  **1 . c 
DO  7f0  K=  1 , N 
HRITE(6, 10) 


0 000  20 1C 
9000  20  2C 
C0002030 
00002040 
00002050 
0C002060 
00002070 
0 000208C 
CD002C90 
00002100 


ID  FORMAT  (IH  ,'SAT',5t,'H?S',12X,'RANGE',°X,'ANGLE',8X,' UPDOPPLEPOOOC21 10 

1 • ,»X, ' DN DOPPLER' ,9 X , • SnBC • , 9X , • SOBL ' , 6X , • ST  AT  ION • , //)  C 0002 120 

T=C  (K)  /RTD  00002130 

G = RL(K)/RTD  00002140 

DO  500  1=1,25  "0002150 

=1=1  00002160 

CALL  ELLI  P ( =1 , E ( J)  ,HP(J)  ,TP(J),  A (J)  ,TAU,FSP,R)  000C217C 

CALL  PRIME  (ESP.  H (J)  ,XI  (J)  , R,0.  ,XS,YS,ZS)  00002180 

= S=(ATAN  (YS/XS)  - HPAD*RI)  *RTD  00002190 

Tp(XS.LT.O.)  FS=FS  + P*RTD  00002200 

TS  = RTD*  APCOS  (ZS/P)  00002210 

F=G*HPAD«RI  00002220 

X = RE*SIN  (T)  *COS  (F)  00002230 

Y = PE*SIN  (T)*SIN(P)  00002240 

Z = R S*COS  (*")  00002250 

ACC=  A RCOS ( (X*XS  + Y*YS*Z*ZS) / (RF*R) ) 0000  2260 

RGE=  ((XS-X)**2*(YS-Y)**2*(ZS-Z)**2)**.5  00002270 

ARG=SIN  (ACC)  »R/»GE  00002280 

IF  (APG.GT.  1 .OODPPCD)  ARG=  1 . "00''0009  0 00  0 2290 

D=(APSIN  (APG)-P/2.)*RTD  0C002300 

RT  = SQRT  (R*R-RE»RF)  00002310 

IF  (RGE.  LE.  RT)  D=-D  00002320 

CALL  DOPE  (R,p(J)  ,A  (J)  ,TAU,FSP,HP(J)  ,H  (J)  , XI  ( J)  ,XV,YV,ZV)  00002330 

XT=-HPAP*RE*SIN  (T)  *SI N (F ) /3  600  . 0000234  0 

YT  = HRAD*RE*SIN(T)  *COS  (F)  /3600.  0000  2350 

DOT=  (XV-XT)«  (XS-X)  * (YV-YT)  * (YS-Y)  *ZV*(ZS-Z)  00002360 

DOT=DOT/RGE  00002370 

DPDOP= -DOT*  FR 1/CL  00002380 

DNDOP  = -DOT*  FR2/CL  00002390 

HR ITE (6 , 20) J,RI,PG?,D,UPDOP,DNDOP,TS,FS,K  00002400 

20  FORMAT(1H  ,13  ,F9.  1 ,F17.  3,  F14.  3,  F17.  3,  F17. 3,  F13. 3.F13.3, 110)  0000241O 

500  CONTINUE  00002420 

HP  ITE  (6, 25)  00002430 
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680 


FO  PM  AT  ( 1 H ,//)  00002440 

IF  (J1  . E Q.  c ) GO  TO  7°u  00002450 

WPITF(6  ,689)  CC0C2460 

"08MAT(1H  , • DO  YOU  WANT  THIS  PAPT  OF  "HE  PPOGPAM  TO  CONT INU E? • 00002470 
• 'NTFS  AN  INTEGER  0 IF  THE  ANSWER  IS  NO  OR  AN  INTEGER  1 IF  THE00002U80 


ANSWER  IS  YES  IN  COLUMN  2.’) 


00002490 


READ  150,  NO 

00002500 

IF(NO.EQ.O)  GO  TO  685 

00002510 

WRITE  (6,25) 

00002520 

"OO 

CONTINUE 

00002530 

900 

CONTI NUE 

00002540 

685 

WRITE  (6  ,30) 

00002550 

30 

FORM  AT  ( 1 H ,/////) 

00002560 

L=  M-  1 

00002570 

I"  (L.  EQ.O)  GO  TO  690 

000 02 5 80 

WR  IT  5(6,260) 

00002590 

260 

FORMAT  ( 1H  , 'THIS  PAP"  OF  THE  PROGRAM 

COMPUTES  CROSSLINK  VALUES00002600 

7 FOR  ALL  PAIRS  OF  S AT  EL LIT ES .’,///,  ■ RESULTS  LISTED  IN  ORDEF  OF  AP00002610 
1PEARANCE  AR F: 1 , /,  1 RANGE  3ETWEEN  TWO  SATELLITES  IN  NAUTICAL  MILES' 0C00262C 
1,/,'  POINTING  ANGLES  BETWEEN  SA  TE LLI TE S (C 1 2 &L  1 2 ) IN  DEGREES',/','  D00002630 
20PPLER  SHIFT  EXPERIENCED  AT  SECOND  SATELLITE  DOE  TO  SIGNAL  TRANSMI000O264C 
2TTED  FROM  THE  FIRST  SATELLITE  IN  HZ.',/,*  CROSSLINK  VISIBILITY  WHE00002650 
2RE  0 MEANS  THEPF  IS  NO  VISIBILITY  BETWEEN  SATELLITES  AND  1 MEANS  T00002660 
2HERE  IS  VISIBILITY',///)  00002670 

DO  350  M 1 =1 , L 00002680 

TAU=CON«A  (M  1)  **1. 5 00002690 

N1=M1*1  00002700 

DO  250  M 2 = N 1 , M 0000271  0 

WR I"E (6 , 40)  00002720 

40  "OF  M A T (1H  ,2X , ' HRS ' , 12X , ' RANGE' , 1 1X,  • Cl  2' , 12X, ' L12' , 1 1X,  00002730 

2'CF  DOPPLER'  , 6X  , ' CVTS'  ,6X,  'SAT1',6X,  'SAT  2',  //)  0 00  027  4 0 

TAU2=CON*A  <M2)  *»1 . 5 00002750 

DO  300  1=1,25  00002760 

RI=I  00002770 

NO  V TZ= 1 00002780 

CALL  ELLIP(RI,"(M1)  ,WP(M1)  ,TP(M1)  ,A  (Ml)  ,TAU,FSP,F)  00002790 

CALL  ELLI?(PT,E(M2)  ,WP(H2),TP(M2)  , A(M2)  ,TAU2,F2P,R2)  00002800 

CALL  PRIME  (FSP,  W (Ml)  ,XI  (Ml)  ,R,o.  , XS  ,YS,ZS)  00002810 

CALL  PRIM"  ("2P,  W (M2)  , XT  (M2)  ,R2,0.  ,X2,  Y2  , Z2)  00002820 

TC 1 = X2- XS  00002830 

TC2=Y2- YS  00002840 

TC 3=Z  2-ZS  00002850 

PI  2=  (TC1**2*TC2**2*TC3**2) **.  5 0 0002860 

F12=RTD*ATAN  (TC2/TC1)  00002870 

IE(TC1.LT.O.)  F12="12*P«RTD  00002880 

T1 2= A "COS (TC3/P1 2 ) *"TD  00002890 

VIZ=  (-XS«TC1- YS*TC2-ZS*TC3)/  (R12*®)  0 0002  90C 

I"  (VIZ . GT.  . 0990  07°)  VI  Z = . 999  099  9 00002910 

DE L=A  PCOS  (VI  Z)  00002920 

DELM  = ARSTN  (PE/P)  00002930 

IP(DFL.LE.DELM) NOVIZ=0  00002940 

CALL  POP"  (R,E  (Ml)  , A ( H 1 ) ,TAU,FSP,WP(N1),W(H),XI(M1),XV,YV,ZV)  000029  50 
CALL  DOPE  (R2  ,E  (M2)  , A (M2)  ,TAU2,  F2P,  WP(M2)  , W (M2)  , XI  (M2)  , S2,  U2,  V2)  00002  960 
DOT"  (S 2-  XV)  *TC1*  ( U2-  YV)  *TC2*  (V2-ZV)  «TC?  0000  297  0 

DOT=DOT/P1 2 00002980 

DOP1 2 = - DOT*  FR  3/CL  00002900 

WRITE  (6,5°)  RI,P  12,T12,"  1 2 .OOP  1 2 ,NOVI  Z , Ml  , M2  00003000 

5°  FORMAT  ( 7 H ,F5.  1 ,R17. 3.F14, 3,  E15. 3.E20. 3,  19, 110,  110)  00003010 

300  CONTINUE  00003020 

WR ITE  (6, 60)  00003030 
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r 


60 

FORMAT! 1H  .//) 

00003040 

IF  (J1  .EQ.OJGO  TO  257 

000030  50 

WRITE  (6, 680) 

00003060 

READ  150,  RO 

00003070 

IF (NO.EQ. 0) GO  TO  690 

00003080 

WpITE  (6,25) 

00003090 

210 

CONTINUE 

00003100 

350 

CONTI  ROE 

00003110 

690 

CONTINUE 

000P3120 

END 

00003130 

SUBROUTINE  ORIME  (FS  P,  NS  ,XIS  , X PS  , T PS , XS,  TS,  ZS) 

00003140 

All  =COS  (PSP)  *C0S  (NS  (-COS  (XIS)  *SIR(NS)  *SIN  (FSP) 

00003150 

A12  — SIN  (FSP)  *CQS(WS)-C05  (XIS)  *SIN  (NS)  *COS  (PSP) 

00003160 

A21  =COS  (FSP)  *51 N ( NS)  *COS  (XIS)  *COS  (NS)  «SIH  (FSP) 

00003170 

A22  — SIN  (FS  P)  *SI  N (NS)  *COS  (XIS)*COS  (NS)  • COS  (FSP) 

00003180 

A31-STN  (XIS)  *STN  (FSP) 

00003190 

A32  = SIN(XIS)  *COS(FSP) 

00003200 

XS=A11*XPS*A12*YPS 

00003210 

YS=A21*XPS*A22*YPS 

00003220 

ZS=A31*XPS+A32*YPS 

00003230 

RETURN 

00003240 

FN  P 

00003250 

SUBROUTINE  FILIP (T, E, NP, TP , A , TA 0 , FSP ,R) 

00003260 

P=3 . 141592  6 

00003270 

Z = 2.*P*  (T-TP)  /TAU 

00003280 

P2=  2.  *P 

00003290 

IF(Z.G?.P2)  Z=Z-P2 

00003300 

IF  (Z . GT.  P 2)  GO  TO  2 

0000331 P 

E 1 =Z+?*SI  N (Z) 

00003320 

E2  = (Z*E*  (SIN  (El ))  - (E*COS  (El)  ) *E1>  / ( 1. -E*COS  (El)  ) 

00003330 

Q-0 . 

00003340 

E3=(Z*r*(SIN(F2) ) -<E*COS(E2)  ) *52)  / (1 . -E  *COS  (E2) ) 

00003350 

Q*Q* 1 . 

00003360 

DE-E3-  E2 

00003370 

DE2=DE**2 

0000  3300 

E2-E3 

00003390 

IF(EE2.GT..0000D001)GO  TO  « 

00003400 

”H  = ARCOS(  (COS  ( p2)  -E)  /(  1 . - F*COS  (F2)  ) ) 

00003410 

IF(Z.GT.P)TH*2.*P-TH 

00003420 

FSP-HPFTH 

00003430 

R=(A*(1.-E**2) ) /( 1.*'*COS(FSP-WP) ) 

00003440 

RETURN 

00003450 

END 

00003460 

SUBROUTINE  DOPE(R,E, A, "AO, FSP , NP , NS , XIS , XD , ID , ZD) 

00003470 

P-3.1415928 

00003480 

T-ESP-WP 

00003490 

RD  = (A*E*2.*P/(TA0*(1.-E**2)  **.5))  *SIN  (T) 

00003500 

TD*  (2  • *P/TA  0)  * ( (1 . -S**2) **  (-1.5)  )*(1.*E*COS(T)  ) *•  2 

00003510 

XPC-RD* COS  (FSP)  -R*TD*  SIN  (»SP( 

00003520 

XPC-XPC/3600. 

00003530 

YPC-RD*SIN  (FSP)  *R*TD*COS  (FSP) 

00003540 

YPC-YPC/3  600. 

00003550 

CALL  PRIMS(0.  ,NS,XI.S,XPC,YPC,XD,YD,Zr) 

00003560 

RETURN 

00003570 

END 

00003580 

i 
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APPENDIX  4 


PROGRAM  SATVIZD 

This  program  is  similar  to  SATVIZE  but  it  gives  results  in  double 
precision.  Its  CPU  time  requirements  can  be  more  than  50%  greater 
than  SATVIZE,  so  its  use  has  been  infrequent  and  specialized. 

The  subroutine  ELLIP  listed  here  should  not  be  confused  with 

the  ELLIP  routines  of  the  other  appendices.  It  contains  the  expan- 

r 3 1 

sion  derived  by  F.  R.  Moulton  for  true  anomaly,  and  is  accurate 
only  for  eccentricity  <0.5. 
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*PRECf  QING  r 4£§.  BIANK^NOT  FILMED 


DATA 

WP/O. 

, 0 . , 0 . , 

0..0 

.,0, 

DATA 

TP/0  . 

, o . , — 6 . 

,-12 

• 1 " 

DATA 

X 1/2  3 

.4,23.4 

,23. 

“,2: 

I? 3. 4/ 

DATA 

C /D. 

,47.54, 

E7.2 

3/ 

DATA 

PL/0. 

,288.73 

,242 

.8/ 

W °I  T' 

(6,11 

0) 

110 

FOFHA 

T (1  H 

,///) 

WRITE 

•(6,11 

5) 

115 

"ORHAT ( 7 H 

, 'THIS 

SATFLLI' 

1 F1AXT«U!1 

0?  1° 

SATELLITES 

AN! 

WPTT? 

■(6,12") 

12"' 

p C R "I  A. 

T(1H 

, » rO  ? E 

*CH 

sat: 

IMPLICIT  REAL*S (A-H,0-Z)  00000005 

DIHENSION  »(11)  (I1)  ,«  (11)  ,TP(11)  ,1TI(  11J  ,C(3J  ,RL<3)  0 0000010 

J 1 = 0 00000020 

DATA  A /O. , 227fi.,2276i. , 22767.  ,2276  7.,  3 6140. 36,36140.  38,  00000030 

1 36140.36,36140. 36 ,14?42. 3,  14342. 3/  OOOOOOUC 

DATA  E/P. , . 1, . i, . 1, . 1 , . 1, . 1 , .1  , .1  1 , .1 / 00000050 

DATA  W/2U3.733,24i.  733,243.  733,  2«  3.  733,  243.  733,  243.133,243.73300000060 
2, 243.  "’3  3, 243. 733, 24  3.  733, 24  3.  031/  0000007C 

,0.,0.,0. ,0., -90. 0,-9 0.0/  00000080 

°.  ,0.  ,-12.  ,-24.  ,-36.  ,-.75,-6.75/  00000090 

. 4, 23. 4,  11  3. 4,  11  3. 4,  113.4,  113.  4,63.4,00000100 

00000110 
00000120 
00000130 
000001 40 
00000150 
00000160 

E VISIBILITY  PROG  RAH  IS  WRITTEN  FOP  A00000170 
2 GROUND  STATIONS.*,//)  00000180 

000n01 90 

LL ITF  THE  FOLLOWING  DATA  HOST  BE  SOPP00000200 
1 LIFE:  * ,/,  ' SENINAJOP  AXIS*,/,'  ECCENTRICITY' ,/, • RIGHT  A SCENSION  ' ,00000210 
1/,'  A’GUHENT  CF  PERIGEE',/,*  TIHE  OF  PERIGEE*,/, • INCLINATION ',//)  00000220 
WPITe(6,125)  00000230 

125  "ORNAT  ( 1 H , • FOR  EACH  G FOUNT  STATION  THE  FOLLOWING  DATA  HOST  B00000240 

IE  SUPPLTE*': ' ,/,  * COLAT ITUD'1 ' , /,  ' LONGITUDE*,//)  00  000250 

W°ITE(P,  14  2)  00000260 

EOF  HA  I (1 H , 'ADDITIONAL  DATA  WHICH  HOST  BE  S UPPL IFP: ' , /, • UPLIN00000270 
IK  FREQUENCY',/,'  DOWNLINK  FREQU  ENCY  • , /,  • CROSSLINK  FREQUENCY  • ,//)  00000280 
WRITE(6,  13  C)  00000240 

'ORNA  T ( 1 H , *THF  US'?  IS  ASKED  TO  ANSWER  S EVE"  AL  QUESTION  S .',/,  00000  300 
i*  IF  ”HE  USE®  DO'S  NOT  KNOW  THE  VALUE  OE  AN  INPUT  VARIABLE,',/,'  H00000310 

1'  SHOULD  ENT'?  THE  NUNBEP  40?. 0.  A DEFAULT  VALUE  WILL  THEN  BE  USF00000320 

ID.',/,'  THE  US'F  IS  EXPECTED,  HOWEVER, TO  SPECIFY  THE  NUHBEP  OF  SATE00000330 
1LLTTES  AND  GROUND  STATIONS  TO  BE  CONSIDERED.',////)  O0000335 

WRI-StS,  14')  00000340 

145  'ORHAT (1 H ,'HOW  HANY  SATELLITES  ARE  TO  PS  CONSIDERED?  ENTER  IN00000350 

1TEGEP  WITH  RIGH^HOST  DIGIT  IN  COLUHN  2.')  00000360 

READ  150,  H 00000  3^0 

WRITE  (6, 1 3')  H 00000380 

135  FORH  AT  (1 H , 13)  00000390 

WPITE(6,155)  00000400 

155  'ORHAT (1H  ,'HOW  HANY  GROUND  STATIONS  AR'  TO  BE  CONSIDERED?  ENT00000410 

1 ER  INTEGER  IN  COLUHN  2.')  0C00O420 

READ  150,  N 00000430 

WRT TF  (6,135)  N 000C0440 

150  FORH  AT  ( 12 ) 00000450 

WRITE  (6,  25)  00000460 

WRITE  (6 ,25)  00000465 

WR  ITE  (6  ,95)  00000470 

95  FORK  AT  ( 1H  , 'ENTER  FOLLOWING  DATA  IN  DFCIHAL  FORH  ANYWHERE  IN  T00000480 

1 HE  HIRST  20  COLUHNS.  ' ,//)  00000490 

WP  ITE  (6,140)  00000500 

140  FORH  AT  ( 1H  ,'WHA"  IS  THF  SEHIHAJOR  AXIS  VALUE  IN  NAUTICAL  HILES00000510 

1 FOR  SATELLITE  1?')  00000520 

00  101  1=1, H 00000530 

K= I* 1 00000540 

READ  170,  A (I)  00000550 

IF  (A  (I)  . NF.  400. 0)  GO  TO  605  00000560 

A ( I)  = A (K)  0C000570 

J1=1  00000580 
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13 


i 


5 WRITE(6,  305)  »(I) 

5 FORHAT ( 1H  , F20. 8) 

IF(I.EQ.H)GO  TO  101 
PRINT  310, K 

G FORNAT ( 1H  , *FOR  SATELLITE*, 
1 CONTINUE 

HR  ITS  ( 6,  315) 

5 FORHAT ( 1H  , /,  • HRAT  TS  “HE 

DO  103  1-1, N 
K*I*1 

READ  110,  E(I) 

IF(E(I)  . NF.  400. 0) GO  TO  610 
F(I)«l(K) 

J1-1 

0 W°ITE(6,335)E(I) 

TF(I.EQ.H)  GO  TO  103 
P5IRT  310,  K 

3 CONTINUE 
WRITE  (6,  320) 

FORRAT(1H  ,/, • HHAT  IS  THE 
1 R SATELLITE  1?*) 

DO  104  1-1, R 
K-T*  1 

READ  110,  H (I) 

IF(H(I)  .BE.  400. 0)  GO  TO  615 

ji-i 

5 HPITF(6,305)H(I) 

I p (I . EQ.  N)  GO  TO  1*4 
PRINT  310,  K 

4 CONTINUE 
HRITT  (6,32*') 

5 FORR  AT  (IN  ,/,*  WHA’"  IS  THT 

IS  FOR  SATELLITE  1?') 

DO  106  1-1, R 
K-TRl 

READ  110,  HP(I) 

IF  (HP(I)  . RE.UO0.C)  GO  TO  620 
WP(T)  * HP  ( K) 

J1-1 

1 HRTT5(6,  305)  HP(I) 

Ip  (I.EQ.  N)  G 0 TO  1"6 
PRINT  310,  K 

5 CONTI NUP 

WRITE  (6,330) 

3 RORR AT ( 1H  , /, • WHAT  IS  THE 

1 1T?  1?') 

DO  101  1=1, H 
K-IR1 

PEAD  no,  T P (I) 

IF  (TP(I).NE.4*O.P)GO  TO  625 
TP  ( I)  -’■P  (N) 

J1-1 

5 WRITE  (6,3*5)TP(I) 

IF(I.FQ.H)GO  TO  101 
PRINT  310,  K 
7 CONTINUE 

WRITE ( 6,  335) 

> pORNAT(1H  ,/,•  WHAT  IS  THE 

1T5L1ITE  1?') 

DO  108  1-1, H 


?CCENTRICIlrY  OF  SATELLITE  1?') 


C0000580 

00000600 

00000610 

00000620 

13,'?')  00000630 

00000640 

00000650 

ECCENTRICITY  OF  SATELLITE  1?')  00000660 

00000670 
00000680 
0000069C 
00000700 
0000071  * 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 

SIGHT  ASCENSION  YALOE  IN  DEGREES  FO0000078C 

00000790 

0000*600 

00000810 

00000820 

00000830 

oooooewo 

00000850 

0C000860 

00000870 

00000880 

00000890 

00000900 

ARGUH'NT  OF  PERIGEE  YALOE  IN  DEGRSE0000091 0 

00000920 

00000930 

00000°40 

00000950 

00000960 

000*0970 

00000980 

00000990 

00001C00 

*0001010 

00001C20 

00001030 

TIHF  OF  PERIGEE  IN  HOURS  FOR  SATELL09*01  040 

00001050 
00001060 
00001070 
00001080 
00001  *90 
0C0C1100 
00001110 
00001120 
00001130 
00001140 
00001150 
00001160 

INCLINATION  YALOE  IN  DEGREES  FOP  SA00001170 

OOOOUBO 

00001190 
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r r 


i 


63r 


ice 


mo 


K-I  + 1 

READ  IPO,  XI(I) 

IF(XI  (I)  .NE.4FP.C)G0  TO  0 30 
XI  (I)  = XI(K) 

J1=1 

WRITE  (6  ,3  05  ) XI  (I) 

I F ( I.  EQ  . N)  00  TO  1 09 
PPINT  310,  K 
CONTINUE 
HP  ITE  (6,  340) 


00001200 

00001210 

00001220 

00001230 

00001240 

000C1250 

00001260 

00001270 

00001280 

00001290 


F0PNAT(1H  ,//,•  WHAT  IS  THE  COLATITUDE  IN  DEGREES  FOB  GBOOND  S00001300 
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20  5 
1 r2 

345 


1 TATION  1?') 

00  102  1=1,  N 
K = I*  1 

®EAD  170,  C(I) 

IF(C(I)  .NE.  40''.)  GO  TO  6?5 
C(I)=C(K) 

J1  = 1 

HRITE(6,  305)  C (I) 

TF(I.FQ.W)G0  TO  102 
PRINT  2"5,  K 

FORH  AT  ( 1 H , • FOB  GPOUNO  STATION*  ,13  ,•  7 •) 
CONTI  NUr 
HPITE  (6,345) 


00001310 
00001320 
000C1330 
00001340 
00001350 
00001360 
00001370 
00001380 
00001390 
00001400 
0000141C 
OCOO  1420 
000C1430 


®ORNAT(1H  ,/,*  WHAT  IS  THP  LONGITUDE  IN  DEGBEES  FOB  GROUND  STA00001440 
1TICN  1?’)  000C1450 

DO  1 "9  1=1, N 00001460 

K = I ♦ 1 00001470 

FPAD  1 7 P , Pt(I)  00001480 

170  FOPfl  AT  (F2 C . R ) 00001490 

IF(?L  (I)  .NE. 400.0)  GO  TC  640  00001500 

PL  (I)  = PL(K)  00001510 

J1  = 1 00001520 

640  WPITr(6, 305)  PL  (T)  00001530 

IF(I.FQ.N)GO  TO  10°  00001540 

PRINT  2n5,  K 00001550 

10°  CONTINUE  00001560 

HPI  TE  (6  ,84?  ) 00001570 

940  FORNAT  (1 H ,////,'  ENTER  THE  FOLLOWING  FATA  IN  EXPONENTIAL  FOBH00001580 

2 (FG.  300. 0*10. **6  IF  WPITTFN  AS  300. 006) • ,/, * WITH  THF  FIGHTHOST  00001590 
2FIGTT  APPEARING  IN  COLUFN  10.',//)  00001600 

WDTTE  (6,845)  000C1610 

845  FORN  AT ( 1 H ,’WHAT  IS  THF  UPLINK  FREQUENCY  IN  HZ?*)  00001620 

P=AD  95C,  rR1  00001630 

850  FORNAT  (DIO. 2)  00001640 

IF  ( FP1.EQ.  400.0)  FBI  =300. 0D6  00001650 

WPITE(6,  855)  FP1  00001660 

e55  FORNAT  (1 H ,D20.5)  00001670 

NEITE(6,860)  00001680 

P60  FORNAT  ( 1 H , /,  • WHAT  IS  TH®  DOWNLINK  FPFOOFNCT  IN  HZ?*)  00001690 

READ  85P,  FB2  00001700 

IF(FP2.  FQ.400.0)  PP  2=2  4 5. 0D6  0000171  0 

WRTmE(6,855)  FR2  00001720 

WRITF  (6 ,86F)  00001730 

865  FORN  AT  ( 1H  ,/,*  WHAT  IS  THE  CROSSLINK  FREQUENCY  IN  HZ?')  00001740 

9°AD  950,  FR 3 00001750 

IF  (FB3.FQ.uOO.F)FP3  = 3P.OD°  00001760 

W?  ITF  (6 ,955)  PR  7 00001770 

T'tJI.EQ.O)  GO  TO  660  00001780 

WRITE  (6  ,67P  ) 000C1790 
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FORN  AT ( 1 H ,//,'  YOU  HAFE  NOT  USED  A COHPLETE  DATA  SET.',/,*  TH00001800 
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* EPFFORE,  THE  RESULTS  ARE  FOP  D7N0NSTRATT0N  PURPOSES  ONLY.')  C00C1810 

GO  TO  675  00001820 

5 BP  ITF  (6,665)  00001830 

5 FORM  AT  ( 1 H , //,  • YOUR  DATA  SET  IS  CO  HPLi'TE  • • ) 00001840 

5 BRITT  (6,105)  00001850 

5 FOPH  AT  (*  H ,///)  00001860 

P*  3.  1415926  00001870 

RE*3449.  00001880 

PTD*57. 29577951  00001890 

BPAD*15.0/PTD  00001900 

CL=1. 61984*10. **5  00001910 

CON=6. 987*10.**  (-6)  00001920 

BEITE  (6, 190)  00001930 

3 EORHAT (1H  , 'RESULTS  LISTED  in  ORDER  OF  APPEARANCE  ARE:',/.'  RA00001940 

2RGE  BFTBEEN  GROUND  STATION  AND  SATELLITE  IN  NAUTICAL  NILES*, /,*  EL00001950 
2 FT  A""  ION  ANGLE  IN  DEGREES  FPO!l  GROUND  STATION  TO  SATELLITE' ,/,  • UPL0000196C 
21 NK  DOPPLER  SHIET  IN  HZ*,/,*  DOBNLINR  DOPPLER  SHIFT  IN  HZ*,/,*  SUB000C1970 
2 SATILL ITE  COLATITUDE  IE  DEGREES*,/,  • SUBSATELLITE  LONGITUDE  IN  DEG00001980 
2REES',//A  0000 1990 

DP  890  IN=1,N  00002000 

B(IN)=B(IN)/5TD  00002010 

BP(IN)  =BP(IN) /PTD  0 0002020 

YI  (IN)  = XI  (I  N)  /PTD  000C20  30 

3 CONTINUE  00002040 

DO  900  J = 1,N  00002050 

TAU=CON*A (J) ***.5  00002060 

DO  7CD  ,N  O00P2070 

WKI^Efe,  10)  00002080 

FOPNAT  (1 H , •SAT*,^X,*HPS*,12X,*P»NGE*,'>X,*ANGLE*,RX,*0PDOPPLER00002090 
1 • , 8X,  • PNDOPPLFR  * ,e  X , • SUBC  • , °X , *SUBL*,6X,  'STATION*,//)  00002100 

T*C(K)/?TD  00002110 

G = RL  (K) /RTD  00002120 

uo  500  1*1,25  P0C02130 

*>T*I  00002140 

CALL  ELLIP(PI,7(J)  , B P ( J)  ,TP  (J)  , A (J)  ,T  AU , PSP  ,P)  000C2150 

CALL  PRIHE(*5P,B(J)  ,XI(J),R,0.,XS,YS,ZS)  00002160 

?S*  (DAT AN  (YR/XS)-BEAr*PI)  *7TD  00002170 

IF(XS.LT.O.)  FS*FS*P*RTD  00002180 

TS*PTD*DAPCOS  (ZS/P)  00002190 

F=G*BP  AD*P  I 0C002200 

X*RE*DSIN(T)  *DCO*(E)  00002210 

T»R**DSIN  (T)*DSIN  (*)  00002220 

Z = PE*  DCOS  (T)  00002230 

ACC*DAPCOS(  (X*XS*Y*YS*Z*ZS)  /(RE*R))  00002240 

RGE*  ( (XS-X)  **2*  (YS-Y)**2*  (ZS-Z)  *«2)  •*.  5 00002250 

ARG*DSTN  (4CC)*P/FGE  00002260 

I* (ARG.GT. 1 .OOrOORO) ARG*1 . 00000000  00002270 

D*  (DAPSIN  (ARG)  “P/2. ) *PTD  00002280 

PT=DSQRT(R*R-RE*PE)  00002290 

lE(RGE.LE.RT)  D*-D  00002300 

CALL  DOPE  (R , F (J)  , A (J)  ,T  A 0 , PS  P,B  P (J)  ,H(J)  ,XI(J),XT,YV,ZT)  00002310 

XT*-BRAD*RE*DSIN(T)«DSIN(P)/3ffl0.  00002320 

YT»tRAD*RP*DSIN  (T)  *DCOS(P)  /360O.  00002330 

PUT*  (XT-XT)  * (XS-X)  ♦ (YT-YT)  • (TS-Y)  *ZT*  (ZS-Z)  00002340 

DOT*DOT/RGE  00002350 

OPDOP  — DOT*pR  1/CL  00002360 

DNDOP*-DOT*FP2/CL  000023',0 

BPIT*(6,20)J,PT,PGE,D,U  PDOR,  DNPOP,  TS,  PS,K  00002  380 

FOPNAT  (1H  ,I3,F8.1,F1U.3,F14.3,F17.3,F17.3,F13.3,T13.3,I10)  00002  390 

3 CONTINUE  00002400 


n 


i i 


WRITE  (6,25) 

00002410 

25 

FOPNAT(1H  ,//) 

0C002420 

IF(JI.EQ.O)  GO  TO  70p 

00002430 

WRITE  (6,690) 

00002440 

680 

FORNAT(1H  , ' DO  YOU  WANT  THIS 

PART  OF  THE  PPOGHAH  TO  CONTINUE?' 00002450 

1,/,  • ENTER  UN  IN TEGEF  0 IF  TH"  ANSWER  IS  NO  OR  AN  INTEGER  1 IF  THE00002460 


1 ANSW"P  IS  TFS  IN  CCLUNN  2.  ') 


9000247C 


"FAD  15T,  NO 

00002480 

IR(NO.EQ.C)  GO  "0  685 

00002490 

WPITE  (6,25) 

00002500 

"’CO 

CONTINUE 

00002510 

909 

CONTINUE 

00002520 

685 

WPITE  (5  ,3") 

00002530 

3° 

®ORHAT(1H  ,/////) 

00002540 

L = H-  1 

00002550 

IE  (L.EQ.O)GO  TO  69'' 

00002560 

WP  ITE  (6 ,250) 

C00C2570 

260 

EORH  AT ( 1H  ,'THIS  PA®T  OF 

THF  PROGRAH  COHPUTES  CROSSLINK  VALUES0000258G 

40 


50 

300 


7 FOR  ALL  PAIRS  OF  SATELLITES.',///,'  RESULTS  LISTED  IN  ORDER  OF  AP00002590 
1PFAFANCE  AFE:',/,'  RANGE  BETWEEN  TWO  SATELLITES  IN  NAUTICAL  NILES'  C0002600 
1,/,  • POINTING  ANGLES  B"TW"RN  SA  TELLI  IE  S (C 1 26L12 ) IN  DEGREES',/,'  D00002610 
2opPLt’R  SHI1""  EXPERIENCED  AT  SECOND  SATELLITE  DUE  TO  SIGNAL  TRANSHI00002620 
2TTED  "RON  THE  "IRST  SATELLITE  IN  HZ.',/,'  CROSSLINK  VISIBILITY  WHE00002630 
2FE  C 1EA.NS  THEFE  IE  NO  VISIBILITY  BE"WE"N  SATELLITES  AND  1 
2HEPE  IS  VISIBILITY',///} 

DO  3RD  H1  = 1 ,L 
TAU  = CON*A.  (N  1)  **1.  « 

N1=H1 ♦' 

DO  250  H2  = N 1 , N 
W®  ITF  (6 , 40) 

®ORNAT  (1H  , 2X  , 'HRS  • , 1 2X  , ' RANG"'  , 1 1X  , • Cl  2'  , 1 2X  , • L 12  • , 1 1 X , 

2 ' CRDOPPLE®'  ,6X, 'CVTS'  ,FX, 'S  ATI',  6X,  'SAT2',//) 

"AU2=CON*A  ( H 2)  **1.c 
DO  300  1=1,25 
RI=I 
NOVIZ=1 

CALL  ELLIP(RI,®(H1)  ,WP(M)  ,TP(H1)  ,A  (Hi)  ,TAU,PSP,R) 

CALL  ELLIP(RI,®  (H2)  ,WP  (H2)  ,TP(N2)  , A(N2)  ,TAU2,  F2P.F2) 

CALL  PPINE  (FS®,  W (HI)  , XI  (HI)  ,R,C.  ,XS  ,YS,ZS) 

CALL  PRIHE  (E2P.W  (H2)  , XI  (H2)  ,R2,0.  ,X2,Y2  ,Z2) 

TC1=X2-XS 
TC2  =Y2 - YS 
*C  3=Z  2-ZS 

R1 2= (TC1**2*TC2**2*TC3**2) **. 5 
"1  2 = PTD*  DATAN  (TC2/TC1 ) 

I F (TC 1 . LT . 0 . ) F 1 2=F  1 2+P*RTD 
T12  = DARC0S(TC3/R12)  *RT  D 
VIZ  = (-XS*TC1-YS  = TC2-ZS*TC3)/  (R12*R) 

IE ( V IZ. GT. . 4999409) VI Z=. 9999499 
DEL  = DARCOS  (VIZ) 

DELH  = DARSIN  (RE/R) 

IF  ( DEL . LE.  DELH)  NOVIZ  = 0 
CALL  DOPE  (P,E(H1)  , A ( Hi)  ,T  AU  , FS  P,  WP  (H  1 ) , W (H 1 ) , XI  (HI ) , X V,  Y V , ZV)  0900  29  30 
CALL  DOPr  (R2  , E (H2)  , A (H2)  , T AU 2,  F2P,  WP  (H 2)  , W (H 2)  , XI  (H  2)  , S2,  02 , V2)  0 0002940 
DOT*  ( S2-XV)  *TC1  + (D2-YV)  *TC2»  (V2-ZV)  *TC3 
DOT=DOT/P 1 2 
DOP12  = -DOT*  ®P3/CL 

WRITE  (6, 50)  RI,R  12,T12,F  12.DOP12  ,NOVIZ,H1  ,H2 
F0RHAT(1H  ,F5.  1,F  17. 3,  FI  a.  3,  FI  5.3,  F 20. 3, 19, 110,  11°) 

COVTl^lVf 


BEANS  "00092640 
00002650 
00002660 
00002670 
00002680 
00002690 
00002700 
00002710 
00002720 
00002R30 
00002740 
00002750 
009C2760 
00002770 
C0002780 
00002790 
00002800 
00002810 
00002820 
00002830 
C00C2840 
00002850 
00002860 
00002870 
00002880 
00002890 
00002900 
00002910 
00002920 


00002950 

00002960 

00002970 

00002980 

00002990 

00CP3000 
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NB  ITr  (6  ,6C ) C00C301  0 

EORB  A""  ( 1 H ,//)  D0003020 

IF  (J1  . EQ.  r)  GO  TO  2?r'  00003030 

WB ITF  (6  ,689 ) 00003040 

BEAD  150,  HO  00003050 

TF  (NO.  EQ.  0)  GO  TO  6°0  00003060 

NPITE(6,25)  00003070 

CONTINUE  00003080 

CONTINUE  00003C90 

CONTINUE  00003100 

END  0000  31 10 

SUBEOUTINE  PPT  BE  (FS  P,  NS  , XIS  , XPS  , YFS  , XS,  TS  ,ZS)  00003120 

THPLICI?  PEAL*0  (A-H.O-Z)  00003125 

A 1 1 = OCOS  ( FSP)  *DCOS(NS)  -DCOS  (XIS)  *PSIN  (NS)  *0  SI  N (PSP)  0000  3130 

A12=-DSIN(FSP)  *DCOS  (NS)  -DCOS  (XIS)  *DS  IN  (N  S)  * DCOS  (FSP)  0 0003140 

A21  =ECOS  (FSP)  *DSIN  (NS)  ♦ DCOS  (XTS)  *DCOS  (NS)  *DSIN  (FSP)  00003150 

A22=-DSIN  (FSP)  * OSIN  (NS)  ♦ DCOS(XI°)  *PCOS(NS)  *DCOS(FSP)  0 000  3160 

A31  = C°IN(XIS)  *DSTN(*SP)  00003170 

A32=DSIN  (XIS)  *DCCS  (FSP)  00003180 

XS=A11*XPS*A12*YPS  C 0003 1 90 

YS=A21*XPS+A22*YPS  00003200 

ZS=A31*XPS*A32*YPS  00003210 

BETUBN  00003220 

FNO  00003230 

SUBBOUTINF  ELL  IP  (T  ,E,  NP  ,T  P,  A ,TAU  ,FS  P,  P)  00003240 

IBPLICIT  PEAL*9  (A-H,0-Z)  00003245 

P*3 . 141 5926  00003250 

Z=2 .*?*  (T-TP) /TAU  00003260 

SeB=DSIN(5.*Z)  0000  3270 

S6F=r'SIN  (6.  *Z)  00003280 

S7b  = 0SIN  C’.*Z)  C0003290 

C'B  = DCOS(5.*Z)  00003300 

C6B=DCOS(6.*Z)  00003310 

C7B=EC0S(i.*Z)  00003320 

SB  = DSIN(Z)  00003330 

S2B=DSIN  (2.  *Z)  00003340 

S3B  = DSIN  (3.*Z)  00003350 

S4B=DSIN  (4. *2)  00003360 

CF=DCOS(Z)  00003370 

PSP=Z*2.*E*SB*1.2  5*  (E**2)*S2B*  ((E**3)/12.)  * ( 1 3.  *S  3B- 3.  *SB)  00003  380 

FSP=FSP*NP*((F**4)  /96.)  *(10?.*S4B-44.*s2H)  00003  390 

p5=  ( (p**5)  /960.)  ‘(I00"1.  ‘S5B-6  4 5.  «S3B*50.*SB)  00003400 

F6=  ( (E**6)  /°60.  )*  (122  3.  *S 6B-°0 2. *S 4B *85 . *S2B)  00003410 

p->=  ( (?**i)  /32256.J  * (4  727  3 .*  S7N-  41699.  *S5B*  5995 . *S3B*"»4  9 . *CB)  0 0003420 

FSP=FSP*F5*F6*F"’  00003430 

"=  (A*  (1.-F**2))/(1.*E* DCOS  (FSP-NP) ) 00003440 

PFTUPN  00003450 

END  000C3460 

SUBROUTINE  DOPE  (P  , E , A , TAU  , ES  P,  N P,  NS  , XIS  , X D,  Y D,  ZD)  00003470 

IBPLICIT  REAL*R(A-P,0-Z)  00003475 

P«3. 1415926  00003480 

T»»SP-NP  00003490 

RD*(A*E*2.*P/('AU*(1.-E**2)  ** . 5) ) *D  SI N (T)  0 000  3500 

TD*  (2  • *P/TA0)  *(  ( 1.  -P**2)  **(-1 .5) ) • (1 . *S*DCOS  (T)  )**2  00003510 

XPC *PD*DC OS  (FSP)  -P*TD*USIN  (FSP)  00003520 

XPC«XPV3600.  00003530 

YPC  = PD*DSIN(PSP)  *R*TD*DCOS(FSP)  00003540 

YPC*YPC  /360r*.  00003550 

CALL  PR  IS  F (0 . ,NS,XIS,  XPC,  YPC , XD  , YD , ZD)  00003560 

PETUPN  00003570 


m*+* 


APPENDIX  5 
PROGRAM  PERTP 

PERTP  calculates  the  secular  changes  In  argument  of  perigee, 
right  ascension,  eccentricity,  and  semi-major  axis  as  a function 
of  lunar  perturbations.  A normalized  apogee  distance  (V5>  is 
plotted  on  the  CALCOMP  plotter  as  a function  of  time. 

Unlike  Appendices  1 through  4 which  require  distances  in  nauti- 
cal miles,  PERTP  requires  distances  in  kilometers. 


Ill 


//PTRTP  JOE  (6’6P,DQ1 .DESK) , 'EDELMAN  F • , NOTIFY  = TSO 1 41 , 

//  CLASS=F,TI?|s  = 2,TYPPUN=HOLD 
/*  SETUP  TAPE  PLXXXX 
//  FXFC  FOFTGCLG 
//FORT. SYFIN  IP  * 

C THIS  PPOGRAW  PRODUCES  THE  SECULAR  VARIANCE  OF  , AN  INITIALLY, 

C HIGH  ALTI?Urr  ORBIT,  AMD  AT  THE  SAME  TIME  GENERATES  STABILITY  OF 
C ECCENTRICITY,  EIGHT  ASCENSION  AND  ARGUMENT  OF  PERIGEE.  IT  ALSO 
C PPOPUCFS  THE  '■VEP  CHANGING  SEMI-MAJOR  AXIS. 

C AND  THEN, USING  THE  CALCOMP  PLOTTER 
C PRODOCrc  PIOTS  FOE  DIFFEFENT  ANGLES 
C OF  INCLINATION, GRAPHING  VS  VS.  DAYS. 

C ALL  VALUES  IN  THIS  PROGRAM  HAVE  BFFN 
C CONVERTED  to  THE  B'TPIC  systfh, 

C WHERE  1 NAUTICAL  NILF=1.852  KM. 

DIMENSION  T3US(2090)  .XAPPAY  (42)  ,YAR°AY(42) 

CON  = 2.  T72A«*’o.  **(-6) 

PI  = 3. 1 41  592  e 
oTr=s-’>?cFnosi 
PHO= 393364. 0009 
DO  2C ? J=“S,1?5,4S 
EJ=J 
XI  = PJ 

AO=256452 . 9nQ5 
TAU=CON*AO**1 . 5 
V=(24./taU) 

EO=. 15 

VcO-AO*  (1  . A'O) 

XI  = XI/P'"D 

w=c . 

WP=45 . 

W=W/’TD 

WP=WP/RTD 

CALI  APGPEP  (AO  , XI  , WP  , WS  , V 1 ) 

IFJV.GT.O.)  WP  = WS 
XID=XI*?"D 
NF  ITr  (6 , 1 0) 

10  FORM  A"  ( 1 H ,9X,»?',r,y,.nr.,  15X,,MPD',1',X,'DF»,11X,*V5',15X, 

1 1 D WPD  • , c X , ' XT  • ,6  X , ' D A.  Y S • ) 

MPITE  (6,15)  XID 
15  form  A*"  ( 1H  , 1 0 5X  , F6 . 2) 

OO  1 Of.  T;1 ,40 

PT  = I 

IF  ( I.  FQ  . ’)  E = EO 
DAYS=PI*10. 

X APR  AY  (I)  =DA  YS 
IF  (I.  ST.  1)  AO  = A 
CALI  MOON  (AO,  XT,  A.) 

UPA3=PI*  (1  . /PI.  )*  ( (A./PHO)  «*3) 

A2=  ( (A/P HO)  **2) 

A4=  ( (A/RHO)  **4) 

CT2  = (COS  (XI)  )*»2 
Cl 4=  (COS  (XI)  ) **4 
SP2=  (SIN(WP)  ) **2 
SI2  = (SIN  (XI))*»2 

PAP1=((-135./12«.)  ♦(315./12S.)*CI2) 

PAR2  = ( (2625./20  48.)  -(■’H75./10  24.)  *CI2 ♦ ( 1 1 325  . /204 8 . ) *CT4) 

PAR  3=  ((3 16. /I  2°.)- (3 15. /1 6.  ) *CI2*  (2  20  5. /I  28.)  *CI4) 

PAP 4=  (3.-  (1S./2.)  *SP2*EI2) 

PA»s,(  (22S./32.)  -(315./P.)  *SP2) 

PA?6  = ( (-45.  /’2.)  ♦ (315./64  .)  •SP2*PAR5*CI2*  (2205./64. ) *SP2*CI4) 


0000001 0 
00000020 
00000030 
00000040 
00000050 
0000006C 
00009C70 
00000080 
00000090 
000C01 00 
00000110 
00000120 
00000130 
00000140 
00000150 
00000160 
00000170 
00000180 
000001 90 
00000200 
00000210 
00000220 
00000230 
00000240 
00000250 
00000260 
00000270 
00000280 
0000029C 
00000300 
0000031 0 
0000032C 
00000330 
00000340 
00000350 
00000360 
0900C370 
00000380 
00000300 
0C0C040C 
00000410 
00000420 
00000430 
00000440 
00000450 
00000460 
00000470 
00000480 
0 0000490 
00000500 
90000510 
00000520 
00000530 
00000540 
00000550 
00000560 
00000570 
00000580 
00000590 
00000600 
00000610 
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DW=-2*UPA?*C0S  (XI)  *((3./4.)*A2*PAS1*(AU*PAR2))*V 

DE  = -0PA3*E*SIN  (2*WP)*  ( (-15. /4. ) *SI2*A2*PAR3)  *¥ 

DWP=UPA3*(PAR4*A2*PAP6)  *V 

W=W*PW 

E=**DE 

V 5=  (A*(1.  *E)  ) /Vc0 
WP=liP*DWP 
HD=W*P7D 
WPD=WP*RTD 
DWD*DW*RTD 
DWpr»DWP«RTD 
YA?»AT  (I)  =Ve 

WRITE  (6,20)  E,WD,WPD,DE,  T5,DW  PD,  CATS 
2?  TOPHAT(1H  ,6E17.7,pr>.  2) 

100  CONTINUE 

WRIT'  (*  ,2 ') 

25  EOEH  AT  (1  R , / /) 

CALL  PLOTS  (IBOp,2CCP,6) 

CALL  PLOT  (n  . 0 , 1 . A ,-3) 

CALL  SCALv(XAPRAY,10.rt,tt0,1) 

CALL  SCALE(YARP  AY,8  .0,  40,  1) 

CALL  AXIS  (0 .0,0  ,l2HTIRr  II  OATS , - 1 2,  1 0 . 0 
1 ,0.o,XAFRAY(41)  ,TAPPAY  (42)) 

C DELTA  APOGEE,  IS  THE  NEW  SEHI-MAJOR  AXIS 
C TI"ES,  THE  NEW  ECCENTRICITY  PLUS  ONE,  DIVIDED 
C BY  THr  INITIAL  VALUE  OP  THE  SE.NI-NAJOR  AXIS 
C T IH ES , THE  INITIAL  ECCENTRICITY  PL«S  ONE; 

C OR  IP  EQUATION  PCI,  (A*  (1 . *E)  ) / (AO*  (1 . *10)  ) 

CALL  AXIS  (0.0, 0.0,1 2HTELT  A A POGEE,  ♦12,0.0 
1 ,9C.o,  YAR»AY(41)  ,YAFFA  Y ( «2)  ) 

CALL  LINE(XARRAY,YARPAY.4C,  1,0,C) 

CALL  PLOT  (1R.o#,-3C.o,-7) 

200  CONTINUE 

CALL  PLOT(12. 0,0. 0,999) 

PE  TURN 
END 

SUBROUTINE  A RG PER  (AO,  XI  ,WP,WS,V1) 

PI=3. 14159265 
RHO*  3P3364.0000 

PTD=57. 29577°ei 
»1  =0. 

WS  = o. 

A2=  (AO/RHO)  **2 

V>  = xi 

IF  (V 4.  EQ.  0.)  V 4 = . 00A 1 

NUN=-3.*A2*  ( (45./32.)  - (225.  /32.  ) * (COS  (T4)  ) **  2) 

DEN  = -i.  5*  (SIN  (V4)  ) **2  + A2*  ( (3  1 5 . /64  . ) - (3  15. /B. ) * (COS  (»4)  ) **2 
1 ♦ (2205./64.)  *(COS(V4)  ) **2) 

PAT2  = NU(1/DEN 
IF  (RAT2  .LE.o  . ) RETURN 
JE  ( P AT 2.  GT.  1 .)  RECURS 
A RG  = SQPT( RAT2) 

WS=  ARS  IN  (APG) 

V 1=  1. 

RETURN 

END 

SUBROUTINE  BOON  (AO,  XI,  A.) 

PI*3. 14159265 
P2=PI/2 
T*».  3 

NUE=. 3986*10. **6 


00000620 
00000630 
00000640 
00000650 
00000660 
00000670 
00000680 
00000690 
00000700 
00000710 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 
00000780 
00000790 
00000900 
00000810 
00000P2C 
00000830 
0C00U8U0 
00000850 
00000860 
C 000087  0 
000C0880 
00000890 
00000900 
000C0°10 
0C000920 
00C  0093  0 
00000940 
00000950 
00000960 
00000970 
000C0980 
00000990 
00001000 
00001010 
00001020 
00001030 
00001040 
00001050 
00001060 
00001070 
00001 080 
00001090 
00001100 
00001  110 
0000112C 
00001130 
00001140 
00001150 
00001160 
00001170 
00001180 
00001190 
00001200 
00001210 
00001220 
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pm=383364.9900 

C2=(4.  8998*10.**3)/((PM-AO)  **2) 

WS*63n  . 7828  * (AO**  (-1 . 5)  ) 

WN=2. 66381*  (19.**  (-6)  ) 

W*  (WS*COS  (XI)  -MP) 

IF (M.LT.P.) K=-W 

W?BL  = ( (W)  **  2 + (MS* SIN  (XT)  ) **2) 

NRFI=SQR?  (MPEL) 

PAT=AO/  (P «- AO) 

RAP=.3?333333*STN  (TE*PAT)»  ( (COS  (TF*RAT)  ) **2+2.) 
DVR=2.  *C2*PAP/(WR?L*RAT) 

DA  = ( (DVR)  **2)  * (A 0**2)  /NOE 
V=SIN  (XI) 

TF(V.EQ.C.)  V=.OOCnOOOC1 
ANG=TE/V 

IF ( ANG . GT. P2)  ANG=P2 

DAD=96490.  * (MR  EL* ANG/ (PI**  2)  ) *9  A 

DADY=365. *DAD 

P AT  A=  DADY/AO 

A*A0*’9. *DAD 

RETURN 

EN  D 

//LKFD.SYSL1B  DD  DISP=SHR 
//  DD  rSN  = SYS1  .CALCC«PrDTSP=SHR 
//GO.  PLOTTAPE  DD  DSN  = PL0'",DISP=  (.KEEP)  , 

//  0NTT=(','APF7,  .DEFEP)  ,DCB=DEN=1, 

//  VOL=SE»=PLTXXX, LABEL*  (,NL) 


00001230 
000012110 
00001250 
00001 260 
00001270 
00001280 
00001290 
00001300 
00001310 
00001320 
00001330 
900C1340 
0C001350 
00001360 
00001370 
00001380 
00001390 
00001400 
0000141 0 
00001420 
00001430 
000C1  440 
C0001450 
00001460 
00001  470 
00001480 
00001490 


APPENDIX  6 


PROGRAM  SATLUNAE 


SATLUNAE  has  the  same  purpose  as  SCOREE,  (App.  7),  but  approx- 
imate integrations  rather  than  the  time  iterations  of  SCOREE  are 
used.  Its  utility,  then,  is  for  time  >5  years  or  when  very  strin- 
gent CPU  requirements  are  imposed. 

The  future  time  of  interest  for  the  perturbed  elements  is 
entered  on  line  350  in  days. 

SATLUNAE  also  Includes  a more  convenient  crosslink  pointing 
angle  coordinate  system  then  SATE.  See  Figure  11.  TU12  and  FU12 
are  chosen  in  a local  satellite  coordinate  system. 
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//T  S 0 4 2 n A JO  3 (6360, D^l, DESK)  ,• CHRIST OPHE B P,,NOTIFY=TS342C, 

//  TIME=2 
//  EXEC  FOPTGCG 
//FOFT.SYSIN  Df  * 

C THIS  S AT  ELL  ITT  VISIBILITY  PROGFAM  IS  WRITTE N *OR  A FA  XI HOH  OF 
C IFF  SATELLITES  AND  TWO  GROUNr  STATIONS. 

C THE  APR  AYS  CONTAIN  THE  FOLLOWING  I NFOPK  ATION  FOR  EACH  SATELLITE: 

C A — SEHIMAJOR  AXIS  IN  KM 
C F--ECCENTPICITY 
C W--  FIGHT  ASCENSION  IN  DEGREES 
C WP--APGUMENT  OF  PERIGEE  IN  DEGREES 
C TP--TIME  OF  PEPIGEE  IN  HOUFS 
C XI  — INCLINATION  IN  DEGREES 

C ARRAY  E CONTAINS  THE  COLATI'UTE  AND  LONGITUDE  rOa  EACH  GROUND  ST  AT 
DIMENSION  A (10)  ,E0  (10)  ,»0  (10)  ,WP0  (10)  ,TP(10)  ,XI  (1  0)  ,B  (4) 

DATA  A/2656  1 . , 2656 1 . , 4 2 1 6 4.  , 4 2 1 64.  , 4 2 16 4.  , 42 1 64. , 

1 106240., 106247., 1C62 4"., 1C6 247./ 

DATA  EC/.  -’2  5,.  125,0.  ,'\,0.,0.,D.,C.,r.,0./ 

DA” A HO/0.,270.,0.,0.,0.,0., 

4 9.,0.,0.,0./ 

DATA  WP0/-9A.,-QD.  ' ,0.,n.,0.,”.,C.,C.,0./ 

DATA  TP/0., -6. ,-3.,-°., -IS. ,-21., 

1 0. ,-24. ,-49. ,-72./ 

DATA  XI/6  3. 4?S,«3.  UAS,''.  ,n.  ,n.  , 3.  ,QC.  ,qo.  ,90.  ,?«./ 

DATA  B/47. 54,  57. 23,  298.73,242.9/ 

WPIT'7  (6,4) 

4 "OFMATCH  ,•  ABTG*  ,13X, 'El',  13X,*W*  , 13X,'HPS’ , 13X,*  V5*  ,/) 

DIMENSION  E (10)  ,WP  (l^)  ,ABIG  (10)  ,E1  (I'')  , «(  1")  ,WPS  ( 10)  , V5  ( 1 0) 
DIMENSION  XM  ( 1 0) 

RTD=57. 29573 
DO  P JR  = ’ ,1" 

NO  (JP)  =V»P  (JF)  /RTD 
WP"  (JP.)  = WPA  (JR)  /PTD 
XM  (JP)  = (XI  (JR)  - 23. 4)  / FTP 
DAYS=4Q0. 

CALL  LUNA.  (A  (JR)  ( JP)  ,W'  (JR)  , W P"  (JR)  ,XM  (JP)  , DATS,  ABIG  (JR) 

1 ,ri  (JP)  ,W(JF)  , WPS  (JP)  ,V5  (JP)) 

WD=W  (JP)  *4TD 
HPSD  = WPS ( JP) *PTD 

WPTT1r(6,0)ABIG(JR)  ,F1  (J3)  , WD,  W PSD,  V 5 ( JP) 

9 FOPM  A” ( 1 H , 1F14. 1, 1F14.9, IF  14 . 6 , IF  1 4 . 6, 1 FI  4 .7) 

0 CCNTINUr 

C NUMS  IS  THE  NUM3FP  OF  SATELLITES  TO  BE  CONSIDERED. 

C NUNG  IS  THE  NUMBER  OF  GROUND  STATIONS  ”0  BE  CONSIDERED. 

NUMS= 1 n 
NUM  G=2 
P=  3.  1 41  5926 
RE  = 6 370 . 8 

°”D=57.P95“t7951 
W4  AA=  1 5 . 0/P  TD 

C CL  IS  THE  VELOCITY  OF  LIGHT  IN  NAUTiefti  MILES  PEP  SECOND. 

CL=2 .90703*10 . **5 

C FP  IS  THF  UPLINK  FREQUENCY  IN  HZ. 

C F»1  IS  THE  DOWNLINK  FPEQUFNCY  IN  HZ. 

FPsACO.PEIO.*** 

'5  1 = 2 45. 0*1 0 . **6 
TON  = 2.7721P*1r.**(-6) 

DO  900  J=1,NUMS 
XI  (J)  =X I ( J)  /PTD 
E (J) *E1  (J) 

A (J)  =ABIG  (J) 


WP(J)  = WPS(J)  00000620 

C TAU  IS  -HE  PERIOD  O'  PO'ATION  O'  SATELLITE  J.  00000630 

T»U=CON*A(J)  **1.5  00000640 

DO  TOO  K=1 , RUNG  00000650 

WRIT'  (6 , 10)  00000660 

10  FOP' AT (1H  , *SA7*,5X, • H' S' , 1 2X , • PA NGE* ,9 X,' ANGIE' ,R X, • UPDOPPLEROOOOO 67C 

1 ' ,8X,'  ONDOPOLER'  ,9 X,  'SUBL',°X,  'SUBC'  ,6X,'  STATION'  ,//)  00000680 

T=B(K)/PTD  00000690 

G=  B (K*  2)  /PTD  00000700 

C I IS  TH'  HOOF.  00000710 

DO  500  1=3,51,3  00000720 

PI=I  00000730 

C SUBROUTINE  ELLIP  COMPUTES  THE  PAHGE  'ROH  G'OCE'TER  TO  SATELLITE  J AND  00000740 
C THE  ANGLE  MEASURED  IH  ORBIT  PLANE ; REFERRED  10  RODE  OF  ORBIT  PLARE  00000750 

C ARD  EQUATORIAL  PLARE.  00000760 

CALL  ELLIP('I,'(J)  ,WP(J),TP(J),A(J),TAU,*SP,P)  00000770 

C SUBROOTIRF  PRIME  COMPUTES  THE  INERTIAL  CARTESIAN  COORDINATES  OF  THE  00000180 

C SAT'LLIT'  J (WITH  GFOCENTFP  AS  THE  ORIGIN).  00000790 

CALL  PPIHE  (FSP.H  (J)  ,XI  (J)  , ?,0.  ,XS,TS,ZS)  00000800 

C FS  IS  THE  SUBSATELLITE  LONGITUDE  IN  DEGREES.  00000810 

'S=  (AT AN (TS/XS)  -WPAD'PI)  'RTF  000009  20 

IF  (XS.LT.O.  )FS='S*P*PTD  0000083C 

C TS  IS  THE  SUBSAT'LL ITZ  COLATITUDF  IN  DEGREES.  000C0840 

TS  = RTD*APCOS(ZS/P)  00000850 

F=G*  WRAD*  RT  00000860 

C (X,Y,Z)-  ARE  T HE  INERTIAL  CARTESIAN  COORDINATES  OF  GROUND  STATION  K.  00000870 

X=R"SIH(T)  »COS  (F)  00000880 

Y = RE*STN(T)*SIN(F)  00000890 

Z = P E* CO S (T)  00000900 

C ACC  IS  THE  ANGL'  BE'W  F'N  GROUND  STATION  K AND  SATELLITE  J FROH  the  00000°10 

C CENTER  OF  THE  EARTH.  00000920 

C RGE  TS  THE  RANGE  FROM  GROUND  S”A'ION  K TO  SATELLIT'  J IN  NAUTICAL  00000930 

C NILES.  00000940 

ACC=A  RCOS ( (X*XS*Y*YS+Z»ZS)/ (RE'R) ) 00000950 

RG'=  ( (XS-X)»'?  + (YS-I)  “2+  (ZS-Z)  "2)  *».  5 00000960 

ARG  = SIN  (ACC)  »R/PG'  00000970 

IF  (ARG.GT.  I.-PDOOCO)  ARG=  1 . 0D',nDr  11  00000980 

C D IS  THE  ELEVATION  ANGLE  IN  DEGREES  FRO*  GRCU'P  STATION  K TO  C0000990 

C SATELLIT'  3.  00001000 

D=  (ARSIN(APG) -P/2.)  *PTD  00001010 

PT=SQPT  (F*R-RE*PE)  00001020 

IF(RGE.LE.RT)  0=-D  00001030 

C SUBPOUTI NE  DOPE  COMPUTES  THE  COMPONENTS  OF  SATELLITE  VELOCITY.  00001040 

CALL  DOP  F (R  , E ( J)  , A ( J)  ,?  AO  , ES  P,  W P (J)  , S ( J)  , X I ( J)  , X V,  Y V , ZV)  0 0001050 

C (XT, YT)  A'E  THE  COnPON'NTS  OF  THE  RELATIVE  VELOCITY  OF  GROUND  00001C60 

C STATION  K.  00001070 

X7=-H?AD*RE*SIN  (T)  *S I N (F) /3600.  C0001C80 

Y”'  = NPAD'RE*SIN  (T)  *COS  (F)  /3600.  3 0001090 

DOT=  (XV- YT)  * (XS-X)  ♦ (YV-TT)  * (YS-Y)  *ZV*  (ZS-Z)  0 0001100 

DO  T=- DO'/P  GE  00001110 

C UPDOP  IS  THE  UPLINK  DOPPLER  SHIFT  IN  HZ.  00001120 

C DND"P  IS  THE  DOWNLINK  DOPPL'R  SHIFT  IN  HZ.  00001130 

UPDOP  = DOT*FP/CL  00001140 

DNDOP  = DOT*FP1/CL  00001150 

WRITE  (6,20)  J,RI,  RGE,  D,  UPDOP,  ONDOP,  FS,  TS , K 00001160 

20  cOPM  AT  (1  H , I3.F8.  1,'17.  3,F14.  3,'17.  3,F1N.  3,  F13.  3,F13.  3,110)  00001170 

50"  CONTINUE  00001180 

NRITF  (6 ,2N)  00001190 

25  FORM  AT  (1 H ,//)  00001200 

700  CONTINUE  00001210 

°0O  CONTINUE  00001220 
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WRITE  (6 ,3C) 

30  FORM  AT (1H  ,/////) 

L=N'1"S-  1 

IF  (L.',0.o)GO  10  350 

C FP  IS  THE  CROSSLINK  fpFQO'NCY  m HZ. 
rP=60.*1P.**9 
HI  IS  THE  FI 3 ST  SATELLITE. 

H2  IS  TRF  SECOND  SATELLITE. 

TO  350  11=1, L 
TA  H=CON*A  (Ml)  **1. 5 
N = M1*1 

DO  250  M2=M,NUKS 
WRITE  (6  ,tin) 

in  POPS  A-1’  ( 1 H . 2X, ' HF E ’ , 1 2X,  ‘RANGE  ' , 5X , ■ LI  2 • ,5X , *C  12 • , 31 , ' FU1 2 < , 

2 8X, ‘TO  1 2 • , 4X , •CPDOPPLEP • ,6X,  'CVIS' , 6 X , ' SAT1 • , 6 X , • S AT2 • ,//) 

TA02  = CON*A(M2)**1.5 
DO  300  1 = 3,150,3 
?I=I 

C NOVIZ  DETERMI NES  CP^SSLINK  VISIBILITY  WHERE  0 MEANS  NO  VISIBILITY 
C AND  1 MEANS  VISIBILITY. 

NOV IZ= 1 

CALL  ELLIP  (FI,F(M1)  ,H  P ( Ml)  ,T  P (M  1)  , A (K  1)  ,TAU,FSP,R) 

CALL  ELL  TD  (?I,F  (M2)  ,HP(M2)  ,TP(M2)  , A(M2)  ,TA02,  F2P.P2) 

CALL  PRIME(FSP,N(M1)  ,XI(M1)  ,R,0.,XS,YS,ZS) 

CALL  PRIME  (F2P,W(r2)  ,XI  (M2)  , R2,  ,X2,Y2,Z2) 

TC1  = X2-XS 
TC2= Y2-YS 
TC3=Z2-ZS 

C F12  IS  THE  RANGE  BETWEEN  TWO  SATELLITES  IN  KM. 

C El  2 AND  T 1 2 ARE  THE  POINTING  ANGLES  BETWEEN  TWO  SATELLITES  IN  DEGREES. 
P12=(TC1»*2+TC2**2+TC3**2)  **.5 
CALL  ONPPIM  (TCI  ,TC2  ,TC3  , FS  P,  W (N  1 ) , X I (N  1)  , R 12,  F012,  ?U  12) 
pi  2 = RTD*A,"AN  (TC  2/TC  1) 

IF  (TCI  . LT.O.  ) F12=F12*P*RTr' 

VIZ=  (-XS*TCi-YS*TC2-ZS*TC3)  /(R12*R) 

T r ( V I Z . GT . . 90909000)  VIZ  = . 09900990 
IF  ( VI  Z.  LT.  -.  90099909)  vi  Z=-.  Q990999Q 
rrL=ARCos  (viz) 

DELM  = AR  SIN  (RF/P) 

IF  (DEL.  LE.  DELM)  NOVIZ=D 
IP3=TC3/P 1 2 

IF  (TR  3 . GE.  1 .)  TR  3= . 999 oc 000 
IF(TP3.LE.-1.)  TK3  = -.  09909990 
T1 2 = ARCOS  (TP3)*PTD 

CALL  D0PE(R,F(M1)  ,A  (Ml)  , TA  0 ,F  SP  , WP  (Ml ) ,W(H1)  , XI  (HI)  , XV  , YV,ZV) 
CALL  DOPE  (R2,E(M2)  , A ( M2)  , T AU2 , F2  P,  W P (M2 ) , W (M2 ) , XI  (M2)  ,S  2 ,0  2 ,V  2) 
D0T=  (S2-XV)  *TC1  ♦ (D2-YV)  *TC2*  (V2-ZV)  *TC3 
COT=DOT/R  12 

C TOP  1 2 IS  THE  CROSSLINK  DOPPLER  SHIFT  IN  HZ. 

DOP12  = -DOT»rF/CL 

WRI"E(6,50)  PI,R  12,F12,T12,FU12,T012,D0P12,N0VIZ,H1,M2 
50  “OPM  AT ( 1 H ,F5.1,F17.3,E9.1,F7.1,F10.1,F6.1,F17.3,I9,I10,I10) 

300  CONTINO? 

WPI”F  (6,60) 

60  FORMAT  ( 1H  , //) 

250  CONTINOF 

350  CONTINOF 

END 

SOBPOUTINF  PRIME  (ESP,  WS  , X IS  , X PS  , YPS  , XS,  YS,ZS) 

A1  i = COS  ( FSP)  *COS  (WS)  -COS(XIS)  *SIN  (WS)  *SIN  (PSP) 

A12  = -SIN  (FSP)  »COS(WS)  -COS  (XI S)  *SIN  (WS)*COS  (FSP) 
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A21*C0S(FSP)  *SIN  (NS) +COS  (XIS)  *COS(NS)  *SIN(FSP) 
A22=-SIN  (PSP)* SIN  (US)  *COS  (XIS)  *COS  (NS)  *COS(FSP) 
A31-SIN  (XIS)  * SIX  (FSP) 

X 32  = ST  N (XI 5)  *COS  (FSP) 

XS=A11*XPS*A12»YPS 
YS=A21*XPS*A22*YPS 
ZS=A  31*XPS*A32*YPS 
RE TOFU 
END 

SUBROUTINE  ELLIP(T,E,NP,TP, A,TAU,FSP,P) 

P=3. 1415926 
Z = 2.*P*(T-"P)  /TAU 
P2*2.*P 

2 IF  (Z.GT.P2)Z*Z-P2 

IF  (Z . GT.  P2)  GO  TO  2 
E1  = 7*E*SIN(Z) 

c2  = (Z*E*  (SIN(*1)  )-  (E*COS  (El)  ) *E1)/(1.  -E*COS  (El)  ) 

Q*0. 

4 E3=(Z*F*(SIN(E2)  ) -(?*COS(E2)  ) *F2) /(I  .-E*COS  (E2) ) 

Q»Q*i  • 

DE=E3-Z2 

DE2=DS**2 

E2=E3 

IF(TE2. G”. . 00*00901) GO  TO  « 

TH*AECOS  ( (COS  (E2)  -F)  / ( 1 . - r*CO  S ( E2)  ) ) 

IF  (Z.GT.P)  TH=2.*P-TH 
F5P=NP*TH 

P=  (A*<  1.-E**2)  ) /(1.+E*COS  (FSP-WP)  ) 

RETURN 
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END 

SUBROUTINE  DOPE  (R,  E,  A , T At) , FSP  , HP,  NS  , XIS  , XD  , YD  , ZD) 

P=3. 1415926 
T=  FSP -HP 

®DS  (A*E*2.*P/(TAU*(1.-E**2)  **.5))  *SIN(T) 

TD=  (2  .*P/TAU)»  ((1.-E**2)»*  (- 1 . 5) ) * ( 1 . +E*COS  (T)  ) **  2 
XPC=R**COS  (FSP)  -P*TD*SIN  (r*=P) 

XPC=XPC/360P. 

Y PC  - PD*  SIN  (ESP)  + P*TD*COS  (FSP) 

YPC=YPC/3  600. 

CALL  PPIBE(C.  ,NS,XIS,XPC,YPC,X*  ,YD,ZD) 

PETURN 

END 

SUBROUTINE  DRPRI H ( TCI,  TC2,  "-C3 , 'SP , NS , XI  S, PI  2 , FU1 2 ,TU  1 2) 
A11=COS(FSP)*C*S(NS)-COS  (XIS)  *SIN(WS)  *SIN('I'SP) 

A12  = -SIN  (FSP)  * COS  (NS)  -CO  S (XIS)  * SIN  ( NS)  *COS  (FSP) 

A 21  =COS  (FSP)  *SIN  (NS)  ♦COS(XIS)  *COS  (NS)  *SIN  (®* SP) 

A22  =-SI  V (FSP)  *SIN(NS)  *COS  (XIS)  *COS  (NS)  *COS  (FSP) 

A3 1 = SIN  (XIS)  *STN  (FSP) 

A32*SIN  (XIS)  *COS  (FSP) 

A 1 3 *SI N (XIS)  PSIN(NS) 

A23  = -SIN  (XIS)*COS  (NS) 

A33*COS  (X  IS) 

Z1»A11*TC1*A21*TC2*A31*TC3 

X1*A12*TC1*A22*TC2*A32*?C3 

Y1«A1?*TCUA23*TC2»A33»TC3 

XU*X1 

YU*-Y1 

ZU--Z1 

PTD*57. 29577951 
P=3. 14159265 


00002130 
00002140 
00002150 
00002160 
00002170 
00002180 
00002190 
00002200 
00002210 
00002220 
00002230 
00002240 
00002250 
00002260 
00002270 
00002280 
00002290 
00002300 
000C231 0 
00002320 
00002330 
00002340 
00002350 
00002360 
00002370 
00002380 
00002390 
00002400 
00002410 
00002420 
000C2430 


50 


100 


1 


FU12=R?D*ATAN  (YU/XU) 

IF(X0.1T.  0.)  PU12=EU12*P*PTD 
ZF=Z0/P12 

IB (ZF.  GE. 1 . ) ZP= . 999999°9 
IF  (ZR.  L E.  - 1 . ) ZP  =-  . °9°99999 
T012  = APCOS(ZR)  * FTD 
RETOPN 
END 

SUBROUTINE  LUNA(A,E9,WA,WP,',XI,T,ABIG,E1,WfWPS,V5) 

P=3 . 1 41 5°26  5 
T=24.  *T 
P2=P/2. 

P^D* 57. 29  57  8 

WH  = 2.  6638*10.  **  (-6) 

WS=6?P.',828*A**(-1.6) 

38=393368. 9 
RUE  = . 3996*10. *»6 
TB  = .3 

IE (XI . LT.  .0001)  XT=.00A1 
IF  (XT.  EQ.  P)  XI  =P-  .90  01 
BAT=A/ (F8-A) 

P AR=.  3 333333  3*  SIN  (TE*cA'r)  * ( (COS  (TE*P  AT) ) **2*2.) 

HPEL=S2?T(  (W'*COS(XT)  -99)  **2*(WS*SIN  (XI)  )**2) 
C2=(9.9Q<50*io.**3)/((5*_j)**2) 

DVP  = 2.*C2*9AP/(WbEL*FAT) 

V7  = SIW  (XT) 

IF  (V.  EQ.  C.  ) V7=.ACC0PA''01 
ANG=TE/V 

IF  (ANG.GT.  P2)  ANG=P2 

PAH=36A0.»  (W?EL*ANG/P**2)*  ( (PV  P**2)  * A*  * 2) /EU  E 
C5=DAH 

DA  Y = °7  60. * D A H 
AB  = »*  (C5*T)/2. 

PHO=RM 

C6=P*(1./91.)  *(1./RHO)  **  3 
C7=2. 77219*10.** (-6) 

S2I=  (SIN  (XT)  ) **  2 
C 21=  (COS  (XI) ) **2 
CUT =C 21**2 
AB2=  (AE/RHO)  **  2 

C9=3.*AP2*(  (-u'./32.)  ♦(225-/32.)  *C2T) 

C10  = -7.5*S2T*A?2*(  (316./6U. ) -(315./9.)  *C2I ♦ ( 220e ./6 u. ) *cui ) 
RATTO=-C°/C10 

C11=-  ( 16./4.)  *S2I*AP  2*  ((315. /I  28.)~  (31  5 . /1 6 . ) *C2 1*  (2205. /1 28. ) 

1 I) 

ups*''. 

IF(PATTO.LT.O.)  GO  *0  5P 
IF  (RATIO. GT.1.)  GO  TO  50 
WPS=AFSTN  (SQPT  (°ATTO)  ) 

E1=F0*EXP  (-C6*(AB**1.5)  * (SIN  (2.* WPS)  )*  (Cl  1/C7)  *T) 

GC  TO  100 
WPC*(P/U.)-.A01 

C9  = 3.-  (15./U.)  *S2I*AR2*(1. 054688-22. 5* C2I+17. 226 56 *C4 I) 
WANG=C6*C8*(AB**1 . 5)  /C7 
C21=2.*C6*  (AB**1.5)  *C9/C7 
C20=  2.  *WPC 

El  =E0*"XP  (-C6*  (A  8**1. 5)  *(C1  VC7)  * (-  ( 1 . /C2 1)  * COS  (C20)  ♦ (1 . /C2 1 ) • 

2 COS(C20*C2  1*T)  ) ) 

rw=-2.  *C6*(AB**3)  *COS(XI)  *(.76*AR2*  (-  (135./128.)*  (315./128.)  * 

3 ) ♦ (AR2**2)  *(  (262e./204  8.)  -C»97  9./102  4.)*C2I*  (173  25./2048.)*C4I 
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1 )) 

B«dO*r«H*T 
ABIG=A+  <C**T) 

V5C*A*  (1  .*BC) 
V5*ABIG*(1.»EH  /V5r 
B1T09N 
BSD 


APPENDIX  7 


PROGRAM  SCOREE 


The  lunar  perturbation  results  of  PERTP  are  combined  with  the 
link  program  SATE  to  calculate  new  crosslink  relations  after  a long 
period  (e.g. , a few  years).  The  future  time  of  interest  is  entered 
on  line  260  as  EPOCH,  in  days. 

Semi-major  axes  are  in  kilometers. 
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//SCORE  JOB  (6360, D°1, DESK)  ,'EDELBAN  E ' , NOTI*  Y*TS  0 1 4 1 , 

//  TIBE=2 
//  EXEC  FORTGCG 
//FORT.SYSIN  OD  * 

C TEN  SATELLITES  AND  TNO  GPOUND  STATIONS. 

C THE  ARRAYS  CONTAIN  THE  FOLLONING  INFORBATION  FOR  EACH  SATELLITE: 

C A — SEHIHAJOR  AXIS  IN  NAUTICAL  BILES 
C E--FCCE*'T  PICITY 
C H— RIGHT  ASCERSION  IV  DEGREES 
C W P-  - ARGHHEN'F  OF  PERIGEE  IN  DEGREES 
C TP — TINE  OF  PERIGEE  IN  HOOES 
C XI -- INC LIN ATT  ON  IN  CEGREES 

C APRAY  B CONTAINS  "HE  CO  LATITUDE  AND  LONGITUDE  FOR  EACH  GROUND  STATION. 
DI  BENS  I ON  AO  ( ID)  , EO  (10)  ,40(10)  , NPO  (1  0)  , TP  (1  0)  , XI  (13  ) , E (4) 

DATA  f o/200n0".  ,20',9,'P.  ,200000.  , 2™ ODD. , 2r'('01n .,  200000 . , 

1 25000C . , 25C0CC.,25OOCC.,25COD0./ 

DAT  A f0/.  1,  . 1,  . i, . 1, .1,  .1,  . 1,.1,.1 ,.1/ 

DA NO/0. ,C. ,C . ,0. ,0. ,0. , 

4 0., 72., 1UN. ,216./ 

DATA  HPO/0.  ,0.,0.,O.,C.,0.,0.,0.,0.,0./ 

DATA  TP/F.  ,-80.  ,-160.  , -40  . , - 1 20 . , - 2C0. , 

1 0., -57. 6, - 115.  2,-  172.5/ 

DATA  XI/23.  4,  23. 4,23. 4,20 3. 4, 233. 4, 203. 4, 113. 4, 113. 4, 113. 4, 113.  4/ 
DATA  B/47. 54, 5". 2 3,259.  24 2.8/ 

C EPOCH  ,IS  THr  TIBS  IN  DAYS. 

EPOCH  = 400 . 

C NUBS  IS  TH r NUBBER  OF  SATELLITES  TO  BE  CONSIDERED. 

C NUHG  IS  "Hr  NUBEFP  CF  GROUND  STATIONS  TO  BE  CONSIDERED. 

NU  FS=  1 C 

NUMG=2 

P=  3.  1 415926 

PE=6370. 8800 

FTD=57. 205-7051 

NPAr=  15.0/RTD 

WPI"F(6,15) 

15  FOPBAT  (1H  ,"2X,'  V5*  ,2UX,*  A*  ,24X,  • S'  , 24X  ,'  H'  , 23X,*  HP') 

DIB  ENS  ION  A (10)  ,F  (10)  ,V  (1  0)  ,NP(10)  , VS(10> 

DO  3 JP«1, 10 
XI  (JP)  =XI  ( JP)  /°TD 
NO  (JP)  =WO  (JP)  /FTP 
NPO  (JP)  =NPO  (JP)  /P "D 

CALL  PERT  (EPOCH  , AO  (JP)  ,EO  (JP)  ,NO(JP(  ,NPO  (JP)  ,XI  (JP)  , 

1 A (JP)  , E (JP)  , N (JP)  ,NP  (JP)  , V 5 (JP)  ) 

NP  = N (JP)  *P"D 
NPD=NP (JP) *RTD 

NETT*  (6, 25)  Y5  (JP)  ,A  (JP)  ,E  (JP)  ,»D,NPD 
25  FORH  AT  ( 1H  ,5E24.9) 

3 CONTINU* 

C CL  IS  THE  VFLOCITY  OF  LIGHT  IN  KILOBETERS  PER  SECOND. 

CL*2. 99994*10.  *«5 

C FP  IS  THE  UPLINK  FRF QUINCY  IN  HZ. 

C FR1  IS  TH*  DO NN LINK  FREQUENCY  IN  HZ. 

FP*?00.0*1C.**6 
RR1=245.0*10.**6 
CON=2. 77218*10. «*  (-6) 

NP ITE  ( 6,  27) 

27  FORI  AT  (1  H ,//////) 

DO  900  J»1 , NUBS 

C TAU  IS  THE  PERIOD  OF  POTATION  OF  SATELLITE  J. 

TAU  = CON*A  ( J)  **1.5 
DO  700  K=1 , NUBG 


( 
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WRITE(6,10)  00000620 

PO  B N AT  ( 1 H , • SAT'  , 5X,'  HP  S'  , 12X,'  RANGE'  , 0 X , • ANG  LF*  , 9 X , • U PDOPP  LER  00000630 

1 • ,8X,'DNDOPPLFR'  , 9X,  ' SU8L  • , °X  , • SUBC'  , 6X  , • STATION'  .//)  00000640 

T=B(K)/BTt)  00000650 

G=  B (K  + 2)  /PTD  00000660 

IS  TH*  HOOB.  00000670 

DO  500  1=3,150,3  00000680 

BI=I  00000690 

SUBROUTINE  ELLIP  COHPUTrS  THE  BANGS  FPOH  GEOCEHTEB  TO  SATELLITE  J AND  00000700 
THE  ANGLF  HEASOBED  IN  OBBIT  PLANE; REFERRED  TO  NODE  OF  OBBTT  PLANE  00000710 

AND  EQOATOB IAL  PLANE.  0000072" 

CALL  ELLIP(BI,E(J)  ,HP(J)  ,TP(J)  ,A(J)  ,'-A0,FSP,B)  00000730 

SOBBOOTINE  PBIFE  COEPUTES  THE  INEBTIAL  CABTESIAN  COOBDINATES  0*  THF  0000P7B" 
SATELLITE  J (WITH  GFOC'NTEP  AS  "HE  OBIGIN)  . 00O007,‘0 

CALL  PPIHE  (PSP,  N (J)  ,XI(J)  , F , 0 . , XS  , T S , ZS ) 00000760 

C PS  IS  THE  SOBSATELLITE  LONGITODE  IN  DEGBEES.  0000077(5 

FS=  ( ATAN  (YS/XS)-WFAD*PI)*RTD  0 000078  0 

IF(IS.LT.O.)  FS  = FS*P*RTD  OOOCP7«0 

C TS  IS  TH*  SOBSATELLITE  COLATITODE  IN  DEGREES.  00000800 

TS  = PTD*ARCOS(ZS/R>  0000081'' 

F=G*WRAD*RI  00000820 

C (X, Y , Z)  ABE  THE  INEBTIAL  CABTESIAN  COOBDINATES  OF  GBOOND  STATION  K.  00007030 

X = RF*SIN  (T)  *COS  (F)  00000840 

Y=RE*SIN(T)  *SIN  (F)  00000880 

Z=p**COS(T)  0 0000  860 

C ACC  IS  THF  ANGLE  BFTNEEN  GROUND  STATION  K AND  SATELLITE  J FBOH  THE  0000C870 

C CENTEB  OF  THE  EARTH.  00000880 

C BGE  IS  THE  BANGE  PBOH  GROUND  STATION  K TO  SATELLITE  J IN  NAUTICAL  00000890 

C HILES.  00000900 

ACC=APCOS  ( (X*XS*Y*YS+Z*ZS)  / (R  E*R)  ) 00000°10 

PGE=  ( (XS-X)  **2*  (YS-Y)  **2*  (ZS-Z1  **2)  **.  S O000O920 

APG=SIN  (ACC)  *B/PGE  C00C0930 

IF (A  EG. GT. 1.0 CO  7" 70) AFG=1. 00001000  000009  90 

C D IS  THE  ELEVATION  ANGLE  IN  DrGRFES  PPO*  GROUND  STATION  K TO  0000095C 

C SATELLITE  J.  "0010960 

D=  (ARSIN  (ARG) -P/2.)  *RTD  000C0°70 

BT  = SQFT  (R*B-PE*RP)  00000980 

IF  C GE.  L E . BT)  D=-C  00000990 

C SUP  BOUT!  NR  DOPE  COHPUTFS  TH*  CONPONENTE  OF  SATrLLIiE  VELOCITY.  00001000 

CALL  DOPE  (R,E(J),A(J),TAU,FSP,WP(J),W(J),XI(J),XV,YV,ZV)  00001C10 

C (XT,  Y T)  ABE  THE  COHPONFNTS  OF  THE  RPLATI  VF  VELOCITY  OF  GPOUND  00001C20 

C STA"ION  K.  00001030 

XT=-NPA7*RE*SIN  (T)»SIN  (») /361D.  00001040 

YT=  WR  A D*R  E*  SIN  ( T)  * COS  (F)/?600.  0C0H05C 

"OT=  (XV-XT)  *(XS-X)  ♦ (YV-YTJ  * (YS-  Y)  *ZV*  (ZS-Z)  0 000  1060 

DOT=-DOT/RGr  0000107" 

C UPDOP  IS  THP  UPLINK  DOPPLER  SHIFT  IN  HZ.  00001080 

C DNDOP  IS  THE  DOWNLINK  DOPPLER  SHIP"  IN  HZ.  00001090 

UPDOP=DOT*FR/CL  00001100 

DNDOP=  DOT*  FB  1/CL  00001110 

WRIT*(6,20)  J,  PI,  BGF,D,  UPDOP,  DNDOP,  FS,TS,K  00001120 

FOBHAT  (1 H ,I3,E8. 1 ,*17. 3,  FI  9.  3,  F 17.  3,  FI  7.  3,  FI  3.  3,  FI  3. 3, 1 10)  00001130 

CONTINU*  00001140 

WRITE  (6,26)  00001150 

FOBHAT  ( 1 H ,//)  00001160 

CONTTNUF  OOOOII’O 

CONTINUE  00001180 

WRITE  (6,30)  00001190 

FOBHA"(1H  ,/////)  000C1200 

t-NUHS-1  00001210 

IP (l . rQ.  0 ) GO  TO  35"  00001220 
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C FP  IS  THE  CROSSLINK  ''P'QUFNCY  IN  HZ. 

FR=6P.  *10.**9 

C HI  IS  THE  FIPST  SATELLITE. 

C H2  IS  THE  S^CON n SATELLITE. 

DO  350  M 1 = 1 , L 
TAU=COE*A  (H1 ) **1  . 5 
N=M  1 +1 

DO  250  M2  = N,NUHS 
WRTTr  (6  ,«0) 

l»0  FORH  AT  (1  H , 2 X , • HR  S • ,1  2X  , * R A NGS  • , 5X  , • L 1 2 • , 5X , • Cl  2*  , 2 3X , 

2 'CPDOPPLER  •,  6X,  TVIS<  , 6X,  • SA  ?1  • , 6X  , • SA T2  • ,//) 

TA02  = C0N*A(M2)  **1.S 
DO  3 p*  1=3,150,3 
PI  = I 

C NOVIZ  DETERMINES  CROSSLINK  VISIBILITY  WHERE  0 MFANS  NO  VISIBILITY 
C AND  1 MFANS  VISIBILITY. 

NOV IZ=  1 

CALL  ELLIP(PI,e  (Mi)  ,WP(M1)  ,TP(M1)  ,A(M1)  ,TAU,FSP,R) 

CALL  ELLIP(PI,F  (M2)  ,K  P ( 12  ) , T P ( M2)  , A (M  2)  , T AU  2,  F2P,  R 2) 

CALL  PRIME(ESF,  W (Ml), XI  (Ml)  , P , ” . , XE  , YS,  ZS) 

CALL  PRIME(F2P,  W (M2)  , XI  (M2)  , P 2 , n . , X2  , Y2  , Z2) 

TC1=X2-XS 
"■C2=Y2-YE 
TC  3 = Z2- ZS 

C R’2  IS  THE  RANGE  BETWEEN  TWO  SATELLITES  IN  NAUTICAL  MILES. 

C FI  2 AND  ”12  ARE  THE  POINTING  ANGLES  BETWEEN  TWO  SATELLITES  IN  DEGREES. 
R 1 2=  (”C1**2*TC2**2*TC  3*«2)  **.5 

CA  LL  UN  PRIM  (TC'  ,TC2  ,TC3  , FS  P,  W (Ml)  ,XT  (Ml)  , PI  2,  E0 12,  T0 12) 

F 1 2=F  TD*  AT  AN  (”C2/TC1) 

IE  (TC  1 . LT . 0 . ) F 1 2=F  1 2*  P*PTD 

VT  Z=  (-XS*TC1-YS*TC2-ZS*TC3)/  (HI  2*P) 

IF (VIZ. GT . . QQooaoooj VTZ=. Q99Q9°°Q 
IF (VI7.LT. - . 0990000 O)  vi Z = - . °99  999°9 
DEL  = A?COS  (V”Z) 

DELM=ARSIN  <?S/F ) 

IF(DFL.LE.rEIM)  NOVTZ  = C 
TP  3 = TC  3 /R1 2 

Tf  |TD 1 . GE.  1 . ) TP  3=  . CQQCQQQQ 
JpjTp1.LE.-1.)  pp  Is—  . OOOOOQOG 

T12=ARCOS(”E‘>)  *P~D 

CALL  DOPE  (R  ,E  (Nl),‘(M’|,rAU,FSP,WP(M1),W(M1),XI(H1),XV,YV,ZV) 

CALL  DOPE  (R  2,  e (M2)  ,«  (M2)  , TA  U2  , E 2P  , WP  (M2  ) ,W(M2)  , XI  ( H 2)  ,S2,U2,V2) 

POT=(S2-XV)  *TC1*  (U2-YV)  *TC2*  (V2-ZV)  *TC3 

DOT=DOT/R12 

C DOP 1 2 IS  THE  CPOS'LINK  DOPPLEP  SHIet  IN  HZ. 

DOP12=-DOT*ep/CI 

WFTTE  (6,50)FI,R12,F12,T’2,EU12,TU12,DOP12,NOVIZ,M1,M2 
so  FORM  A” ( 1 H , E5. 1, FI”. 3,F°. 1 ,FP. 1 ,P1C. 1 , = 6.1 ,F11. 3,19 ,110,1 10) 

3Pr  CONTINUE 

WPITE  (6,60) 

6C  FORM  A”  ( 1 H , //) 

25C  CONTINUE 

35n  CONTINUF 
END 

SUBFOUTINE  PRIME(FSP, WS,XIS,XPS,YPS,XS,YS,ZS) 

A11=COS  (FSP)*COE  (WS)-COS  (XIS)  *S  IN  (WS)  *S  IN  (FSP) 

A12  = -SIN  (FSP)*  COS  (WS)  -COS  (XI S)  * SI  N ( W S)  * COS  ( 'SP) 

A 2 1 = COc  (ESP)  *S IN  (WS)  ♦ COS(XIS)  *COS  (WS)  *SIN  (FSP) 

A22  = -SIN  (FSP)  *SIN  (WS)  *COS  (XIS)  *COS  (WS)  »COS  (FSP) 

A31  =SIN  (XIS)  *SIN  (FSP) 

A 32=SIN  (XIS)  «COS  (FSP) 
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XS=A11*XPS+A12*YPS 
YS=A21*XPS*A22*  YPS 
ZS  = A31*XPS  + A32*YPS 
RFTUFN 
END 

SUBROUTINE  ELI  IP  (T,  E , NP  r TP,  A , TA  U ,F  SP  , F) 

P=  3.  1 4 1 5°26 
Z =2  . * P*  (T-TP)  /TAD 
P2=2.*P 

2 IF(Z.GT.P2)  Z=Z-P2 

IF  (Z.GT.  P2)  GO  TO  2 
E 1 = Z ♦E*S  IN  (Z) 

£2  = (Z*E*(SIN  (El)  ) - (E*COS  (El)  ) *E  1 ) / (1  . -E*COS  (El)  ) 

Q=0 . 

4 F3=  (2 ♦ E*  (SIN  (E2)  ) - (E*COS(E2)  ) *E 2)  / ( 1 . -E*COS  ( E2) ) 

Q'O* 1 . 

DE=E3-E2 
0E2=DE**2 
E2=E  3 

IF  (DE2  .GT.  .OOOOO'K' 1)  GO  TO  4 
TH=APCOS  ( (COS  (E2)  - E)  / (I  . -E*COS  (E2)  ) ) 

IF  (Z.GT.  P)  TH=2.  *P-TH 
FSP  = N®+TH 

F=  (A*  (1.-E**  2)  )/(’.♦!*  COS  (FSP-NP) ) 

FETOFN 

END 

SUBROUTINE  DOPc  (P,  E,A  ,TAU,FSP,NP,NS,XIS  ,XD,  YD,ZD) 

?=  3.  1415926 
T = F SP- NP 

P5=(A*F*2.»P/(TA0*(1.-E**2)**.5))*SIN('r) 

TD*  ( 2.  * P/TAU)  *(  ( 1.-e**2)  **(-1 .5)  ) *(’  . *E*C  OS  (T)  ) **2 
XPC=PD*COS  (FSP)  -F*TD*SIN  (FS  P) 

XPC=XPC/3600. 

YPC  = °D*SIN(rSP)  ♦ P*TD*COS  (FS®) 

YDC  = YPC/36n,>. 

CALL  PRIME  (0.,  NS,  XIS,XPC,YPC,XD,  YD,  ZD) 

pr-Upfl 

END 

SU  BPOUT  IN  r UNPPIN  ('C 1 , TC  2,  TC  3 ,F  SP  , NS  , XI  S ,P1  2 , FU  1 2 ,TU  1 2) 
A1  1=C0S(FSP)  ‘COS  (NS)  -COS(XIS)  *S  IN  (NS)  * SIN  (FSP) 

* 1 2 =-  SI  N (®s  P)  *COS  (NS)  -COS  (X  IS  ) *S  IN  (NS)  * COS  ( FSP) 

A 21  = COS  (FSP)  * SIN  (NS)  ♦ COS(XTS)  *COS(NS)  *SIN(FSP) 

A2  2 = -SIN  (FSP)  *SIN  (NS)  +COS  (XI  S)  *COS  (NS)  * COS  (FSP) 

A31  =STN  (XI  S)  »SIN  (®S®) 

A32  = SIN  (XIS)*COS  (FSP) 

A1  3 = SIN  (XIS)  *STN  (NS) 

A2  3 =- SI  N ( XT  S)  *C"S  ( NS) 

A 3 3 = CO S (XIS) 

Zl  = A11*TC1*A21»TC2*A31*TC3 

X1=A12*TC1*A22*TC2**32*TC3 

Y 1 = A1 3*TC1*A23»TC2*A33*TC3 

XD=  X 1 

Y0=-Y1 

ZU=-Z1 

PTP=57.  2°S77<»S1 
P = 3 . 1 415926s 
FU12=P?D*ATAN  (YO/XD) 

IF  (XU.  LT.O.)  FU12=FU12*P*PTD 
ZB  =ZU/F1  2 

IF(ZP.GE.  1.)  ZP  = .°?900<*°9 


00001840 
00001850 
00001860 
00001870 
00001  880 
">0001890 
00001900 
00001910 
00001920 
00001930 
00001  940 
01001950 
00001960 
C0001970 
00001980 
00001990 
000C2000 
00002010 
00002020 
00002030 
00002040 
00002050 
00002060 
0 000207C 
00002080 
000020°0 
00002100 
0C00211C 
00002121 
">00021  30 
000021U0 
000021 50 
00002160 
00002170 
00002180 
00002190 
0 3 OC  2200 
00002210 
00002220 
00002230 
00002240 
00002250 
00002260 
0 00  0 227  0 
00002280 
000C2290 
00002300 
00002310 
00002320 
00002330 
00002340 
00002350 
00002360 
00002370 
00002380 
00002390 
00002400 
00002410 
00002420 
00002430 


IF  (ZP  . LF.  - 1 . ) ZP  =-  . °999OQC*° 

TU  1 2=  A PCOS  ( Z P)  * pTP 
PPTUPM 

end 

SUBROUTINE  PERT (EPOCH, AO, EO , WO , WPO , XI , 

1 A ,E,W,WP,V6) 

C THIS  SUBROUTINE  ?PODUCrS  t-he  SECULAR  VARIANCF  OF  , A N INITIALLY, 

C HIGH  A LTI'rUDP  ERBI”,  AND  AT  THE  3 A HE  TIME  GENERATES  STAPILITY  OF 
C ECCENTPICITY,  FIGHT  ASCENSION  AND  AH  GO  MINT  OF  PERIG".  IT  ALSO 
C PRODUCES  THE  EV  PF  CHANGING  S^MI-MAJOP  AXIS. 

C ALL  VALUES  IN  THIS  PROGRAM  HAVE  E'EN 
C CONVERTED  TO  THE  MrTRTC  SYSTEM, 

C WHERE  1 NAUTICAL  MILE=1.8S2  KM. 

CON=  2.  77  218*1°. »*(- 6) 

PI  = 3 . 1 41 592 

29S"’7?S1 
RHO=  383R6U.  °00  C 
W = WO 
V P= WPO 

IF  ( vp.  '0.  o.)  HF=.norci 
V6=EPOCH/4°  ? . 

DO  2°0  1=1, 4C 
TAU=C0N*»0** 1 . 5 
V=  (24. /TA  U)  *1  " . 

V=V*V6 

V50=  AO*  ( 1 . ♦ EO) 

CALL  ARGPEH  (A0,XI,WS,V1) 

IF  ( V 1 . GT.  0 .)  WP-WB 
IF  (I.rQ.  1)  E = T" 

C MOON  REQUIRES  ; NEW  S^MI-MAJOP  AXIS 

C AT  EACH  NFW  ITEP  AT  ION . ,"H':'PT'FOFr'  FOP  EVERY  ITFPA'TON  AFTER  THE 
C FIRST,  AO  IS  REALLY  THE  NEW  'A  ■ AS  COMPUTED  BY  THr  PPEVTOOS  ITERATION 
C SO  FOB  THE  SECOND  ITERATION  AC  REALLY  BQUALS  A1,FC>R  "He  THIRD 
C A2  AND  SO  ON.  THUS  “H'  LOGIC  STAT^M^NT  IS  REQUXRED. 

IF  ( I.  GT . 1)  AO  = A 
CALL  MOON  (AO,  XI  , A) 

UP  A?  =PI*  (1  . /PI . ) * ( (A/PHOi)  **  2) 

A2=<  (A/RHO)  **2) 

A4=  ( (A/RHO)  **4) 

CI2  = (COS  (XI)  ) *-*2 
CI4  = (COS  (XT)  ) **4 
SP2=  (SIN(WP)  ) *«2 
ST2  = (SIN  (XI) ) **2 

PAP  1=  ( (-  i?5./12°.)  ♦(31s./128.)*CT2) 

PAR2=(  (2625. /2°4e.)  -(■’875./1C  24.)  *CI2*  (17325. /2049.  )*CI4) 

?AP3»(  (315./128. ) - (31E./16. ) *CI2*  (22rS.  /12F.)  *CI4) 

PAP 4=  (3.-  (IS.  /2.)  *SP2*SI2) 

PAP  5= ( ( 2 2C  . / 32. ) *(Jl1e*./8.)  *SP2) 

PAP6=  ( (-U5./32. ) ♦ (315./F 4.  )»SP2*PAR5*CI2*  (22CS./64.)  *SP2*CI4) 

DW  = -2*UPA3«COS  (XI)*  ((3  . /4  .)  ♦ A 2* PAR  1*  (A4«PAP2)  ) *V 

DE=-U'>A3*F*STN(2*WP)*((-15./4.)»SI2*A2*PAR',)*V 

DWP=UPA3* (PAR4+*2*PARS) »v 

W=W*DW 

E=E*DE 

Vc  = (A*  (1  . ♦ F)  ) /VK0 
WP=  WF  + DWP 
200  CONTINUE 

C DEL'1' • APOGEE,  IS  THE  NEW  SFNI-MAJOP  AXIS 
C TIMES,  THE  NEW  ECCENTRICITY  PLUS  ONE,  DIVIDED 
C EY  ”Hr  INITIAL  VALUE  OF  THE  SEMI-MAJOP  AXIS 


00002440 
00002450 
00002460 
00002470 
00002480 
00002490 
00002500 
00002510 
0°002520 
00002  53  0 
0 00 C 2540 
00002550 
OC002560 
00002570 
O^OOTSSO 
00002590 
00002600 
00002610 
000C2620 
0 00  C 26  30 
00002640 
00002650 
00002660 
00002670 
00002680 
00002690 
C00027CC 
0 000271  0 
00002720 
0 ■'00273  0 
00002740 
00002750 
UC0C2760 
0 000277  0 
00002780 
OOOC279C 
00002800 
00002810 
00002820 
00002830 
00002840 
00002850 
00002860 
000C287C 
00002880 
00002890 
00002900 
00002910 
00002920 
00002930 
00002940 
000029  50 
0C0C2960 
00002970 
00002980 
00002990 
00003000 
C00030  10 
0O003020 
0 COO  30  30 


? 

f 

■ 

j 


1 
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c 

c 


TIMES,  "HE  INITIAL  ECCENTRICITY  PLOS  ONF ; 

CP  IN  pQUATION  FOPM,  (A»(1.=E)  )/(AO*(1  . ♦ ?0 ) ) 

EFTOFN 

END 

SUBPOUTI Np  ARGPEP(A9,XI,NS,V1) 

PI  = 3. 14159265 
RHO=383364.000C 
PTD=57. 29597951 
V1=''. 

WS=0. 

A2=  (AO/PHO)  **2 
VU  = XI 

I=(74.E0.0.)  74=. 0901 

N0M=-3.+A2*  ((45./32.J  *(2  25./32.  ) • (COS  (7  4)  ) **2) 

DEN=-7.  5=  (SIN  (74)  ) **2  + A2»  ( (315./64.)  -(315./8.)  * (COS  (7  4)  ) 
1 ♦ (22C5./64.)  *(C0S(V4)  ) **4) 

EAT2=NnN/DEN 
IF  (F  AT2  . L F.  9 . ) °ETUF  N 
IF  (PAT2.  GT.  1.)  HETOPN 
AEG=SQRT  ( PAT2) 

WS=APSIN  (ARG) 

71=’. 

PETOFS 

END 

SOPPOUTI  NF  MOON  ( AO,  XI  , A) 

PI  = 3.1415°2  85 
P2=PT/2 
TE  = . 3 

N0*=.?986»’C.*=6 
PM =383364.9000 

C2=  (4.  8998  = 10.  » = 3)  /( (FM-AO)  **2) 

HS=6?''.7828»  (AO*»  (-1.  c) ) 

WM=2.  6638  1=  (10. *»  (-6)) 

W=  (WS*C0S  (XI)  -NM) 

Tp ( M . LT. C . ) H=-W 
W R EL  = ( (W)  **  2 = (NS*  SIN  (XI)  ) **2) 

NREL=SQPT  ( NPFL) 

PA  T = AO/ (PM  - AO) 

PAP  = . 3333333 3*  SIN  (TF=FA")  * ( (COS  (TE*PAT)»  *=2  = 2 .) 

D7R  = 2.*C2*PAR/(  WPFL*FAT) 

DA  = ( (D7R) *»2) * (AO»*2) /MOE 
V = S IN  (XI) 

IF  (V.EQ. '’.)  7 = . 9000090 0 1 
A VG=TE/7 

IF  ( AN  G.  G"  . P 2)  AN  G = P 2 

D?D  = 86400.  *(W  PEL* ANG/( PI* *2) ) *DA 

DAPY=36E.*DAD 

EAT  A=DACY/AO 

A=AO+10.*DAD 

PFTOPN 

END 


I ••••••••! 


00003040 
00003050 
00003060 
00003070 
00003080 
00003090 
C90031CC 
00003110 
00003120 
00003130 
00003140 
00003150 
00003160 
00003170 
**2  009C31 80 

00003190 
000C32C0 
00003210 
00003220 
90003230 
09003240 
00003250 
00903260 
00003270 
00003280 
90003290 
90003300 
C 000  3310 
00093320 
00003330 
00003340 
00003350 
000"336C 
00003370 
0°003380 
00003390 
00003400 
00003410 
00003420 
00003430 
00003440 
00003450 
00003460 
00003470 
00003480 
00003490 
00003500 
00003510 
00003520 
00003530 
00003540 
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APPENDIX  8 


PROGRAM  AZ1 


A ground  station  azimuth  angle  calcualtlon  is  added  (lines 
870-1025)  to  program  SATE. 
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ePKECgQIfG  rlfiE  BLANK- NOT  FILMED 


//TS0420A  JOB  (63*0  ,D°1  ,DESK)  , 'CHRISTOPHER  P',  OOOOOOIO 

//  TIME=1  00000020 

//  EXEC  FnRTGCG  00000030 

//FOPT.SYSTN  DD  * 0C000C40 

C AZ 1 (APR  "5)  GIVES  AZIMUTH  ANGLE  000000«5 

C THIS  SATELLITE  VISIBILITY  PROGRAM  IS  WRITTEN  FOE  A MAXIMUM  OF  00000C50 

C TEN  SATELLITES  AND  TWO  GROUND  STATIONS.  00000060 

THE  ARRAYS  CONTAIN  THP  FOLLOWING  INFORMATION  "OP  EACH  SATFLLITF:  0000007C 

A — SEMI MAJOR  AXIS  IN  nautical  MILES  00000080 

E — ECCENTRICITY  00000090 

W — P IGHm  ASCENSION  IN  DEGREES  00000100 

WP — APGUMENT  OF  PEP IG EE  IN  DEGREES  00000110 

TP— TIME  0"  PERIGEE  IN  HOURS  00000120 

XI  — INCLINATION  IN  DEGREES  00000130 

APDA Y B CONTAINS  THE  COLATI"UDE  AND  LONGTTODE  FOR  EACH  GROUND  STATION. 00000 140 
DIMENSION  A (10)  ,E  (1")  ,W  (10)  ,WP(1C)  ,TP  (1")  ,XI(  10)  ,R(4)  00000150 

DATA  A/14342.,14342.,14342.,61421.,61421.,61421.,  00000160 

3 61421.  , 1 1°402.  , 1 1<»402.,  1 19402./  00000170 

DATA  S/.V25,.72?,."'2S.4,.4,.4,.4,.%.5,.  V 00000180 

DATA  »/n/.2’tl.  0,  63. 435, 63. 436,63. 435,  S3. 435,  00000190 

4 24  3.  733,  243.733,243.733,  243.733/  00000200 

DATA  W P/- ° 0 . 0 , - 9 " . D / 0000021 0 

DATA  TP/0. ,-6.0, -12. ,0., -26. 587, -53. 173, -79. 76,d.  , 0000022" 

1 -96. -'0,-192.  16/  00000230 

DATA  XI/63.4i6,63.435,27.4,9D.,o9.,o9wq9.,0.,0.,"./  00000240 

DATA  P/45. 0,0.0/  00000250 

C NtlMS  IS  THE  NUMt>x'R  0F  SATELLITES  ”9  be  CONSIDpFED.  00000260 

C NUMG  IS  'H7  NUMBER  OF  GROUND  STATIONS  to  EE  CONSIDERED.  O0000270 

NHMS=2  00000280 

NU“G=1  000C029C 

P-3 .141 s°26  00000300 

Pi  2=. 5*°  00000302 

P 3 2= 1 . 5*?  0000 C 30  4 

P2=2.‘P  00000306 

88=344".  0 000031 C 

FTD  = 57.  29577051  00000320 

WR  AD= 1 5 , D/RT"  00000330 

S  THE  VELOCITY  OF  LIGH”  IN  NAUTICAL  MILES  PER  SECOND.  00000340 

CL=1.61“8U*10.**5  00000350 

FP  IS  THE  UPLINK  FREQUENCY  IN  HZ.  00000360 

ER 1 IS  THE  DOWNLINK  FREQUENCY  IN  HZ.  00000370 

E?=3"0.r*10.**6  00000380 

p?1=245.0»10.»»6  "0000390 

C0N=6.987*  10.**  (-6)  00000400 

DO  900  J=1, NUMS  00000410 

W(J)=W(J)/RTn  00000420 

WP  (J)  =WP  (J) /P7D  00000430 

XI  (J)  = XI  (J)  /°TD  C0000440 

TAU  TS  THE  ospiOD  Op  ROTATION  OF  SATELLITE  J.  00000450 

TAU=CON*A (J)** 1. 5 00000460 

DO  700  K = 1 , NUMG, 2 00000470 

WRI TE  (6 , 1 0 ) CCCO04e0 

1"  FORMAT!  1H  ,’SAT>  , SX,  ■ HRS*  ,12X,*  RANGE*  ,OX,*  ANGLE  • ,8X,  'AZIMUTH',  00000490 

1 9 X , •UPDOPPL'-R',  RX,  ' DN  TOPP  LEP  • ,9X,  • SOBL*  ,°X,'  SUBC'  ,6  X , • ST  ATION  • ,//)  0O00050C 
T=B  (K)  /RTD  00000510 

K 1 =K ♦ 1 00000520 

0=B(K1)/PTD  00000530 

I TS  THE  HOUR.  00000540 

DO  500  1=1,25  00000550 

°I=I  00000560 

SUBROUTINE  EL L IP  COMPUTES  THE  RANGE  FROM  GEOCENTER  TO  SATELLITE  J AND  00000570 


C CL 
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C THE  ANGLE  BEASURED  IN  ORBI7  PLANS ; REFER RE  D TO  NODE  OF  OREI'  PL  ME 
C AND  EQUATORIAL  PLANE. 

CALL  ELLTP  (PT,F  (J)  , HP  (J)  , TP  (J)  ,A(J)  ,TAU,FSP,P) 

C SUBROUTINE  PPTNE  COBPO-ES  THE  INERTIAL  CARTESIAN  COORDINATES  OF  THE 
C SATELLITE  J (WITH  GEOCENTER  AS  TH»  0°IGIN)  . 

CALL  PR  IB  F (FSP  , W (J)  ,XI  (J)  , R , 0 . , IS , YS  , Zc) 

C FS  IS  '■HE  SOBSATELLITE  LONGITUDE  IN  DEGREES. 

FS=  (ATA  N (TS/XS)  -NRAO*RI)  *P'D 
IF(XS.LT.O.)PS*,,S*F*R?D 
11  IF  (FS.LT. -360.)  FS=FS*360. 

IF  (ES.  LT.-36".)  GO  TO  11 

C TS  IS  THE  SUBSATFLL ITE  CCLA.TITUDF  IN  DEGREES. 

TS  = RTD*ARCOS  (ZS/P) 

f=g*hpad*ri 

C (X,Y,Z)  ARE  THE  INFFTIAL  CARTESIAN  COORDINATES  OF  GROUND  STATION  K. 
X = °E*SI  N (■)  *CO  S ( F) 

Y = RE*SIN  (T)*SIN  (F) 

Z = RF*  CO  S (■) 

C ACC  IS  -H'1  ANGLE  RETWEEN  GERUND  STATION  K AND  SATELLITE  J 7ROH  THE 
C CENTER  OF  THE  EARTH. 

C PGF  IS  THE  RANGE  *ROB  GROUND  STATION  K ■0  SATrLLITE  J IN  NAUTICAL 
C BILES. 

ACC=A  RCOS  ((X*X$*Y*YS*Z*ZS)/  (RE*P)  ) 
PGF*((XS-X)**2*(YS-Y)**2*(ZS-Z)  **2)  •*.' 

A~G=SIN  (ACC)  *P/PG? 

IP  (ARG.GT.  I.d^ddodD)  ARG*1.TCOnODD'' 

C D IS  THE  ELEVATION  ANGL'  IN  ’'FGEE'S  EFOB  GROUND  STATION  K TO 
C SATELLITE  J. 

D=  (APSIN(ARG)  -P/2.)  *PTD 
RT=SQR?(P*P-Pr*RE) 

TF(?GF.LT.?T)  D=-D 
C AZIFUT  IS  THE  AZINUTH  ANGLE. 

SL  AT=  TS/PTD 
SLON**S/RTD 
APG=G-SION 
ALPHA*  (SLAT-T)/2.D 
BETA*  (SLAT  + T)  /2.F 
GAHBA*ABS  (APG)  /7.D 

IF  (GAB  HA.  GT.P1  2.  AND.  GAB  BA.  LT.  P)  G ABB  A=  P-G  ABB  A 
ETNA- SIN  (ALPHA) 

SINR  = STN  (BETA) 

COSA  = COS  (ALPHA) 

COSB=COS(BETA) 

COTG*COTAN(GAB"A) 

U1  = SINA*C0TG/SINB 
ALPHA  = ATAN  (U’J 
U2  =COSA  *COTG/CCS  B 
PFTA*  ATAN (02 ) 

IF  (BE~A.GT.P32  .AND.  BETA  .LT.P2)  BFTA*BFTA*P 
IF  (BETA  .GT. -P12  .AND.  EFT  A . LT . 0 . ) BET  A«  BET  A>P 
AZ  IB  UT=  (ALPHA  ♦BETA)  *RTP 
AFG*S IN (APG) 

IP  (APG.  GT.  0 . C)  A ZI  BUT* 360.  AZ IHUT 

C SUBPOUTINF  DOPE  COHPUTrS  THP  COHPONFNTS  OF  SATELLITE  VELOCITY. 

CALL  DOPE(R,E(J)  ,A  (J)  , TA  U ,F  SP  , WP  (J)  ,i(J)  ,XI(J)  ,XV,YV,ZV) 

C (XT,  YT)  A°F  THE  COHPONFNTS  OF  THE  REL  AT  IF  I VELOCITY  OF  GROUND 
C STATION  K. 

XT*-WRAD*PE*SIN  (T)  • SIN  (F)  /3600. 

YT*WPAD*PF*SIN(T)  *COS(F)  /3600. 

DOT*  (XV-XT)  • (XS-X)  ♦ (YV-YT)  * (IS-Y)  ♦ (ZV)  * (ZS-Z) 

DO?*-DOT/PGF 


C0000580 
00000590 
0C0CP6D0 
o CDC06 1 0 
0O00062D 
00D0D630 
00000640 
00000650 
000DC660 
00000662 
00000664 
00000670 
00000680 
00000690 
0C000700 
000C07  10 
00000720 
000G0730 
0C0C074C 
00000750 
00000760 
0 C000770 
OC0C0780 
00000790 
0C0C0P0C 
00000610 
0000082  0 
00000830 
00000840 
0COOOP5O 
000C086C 
00000870 
00000880 
00000890 
00000°00 
00000910 
00000920 
00000930 
00000932 
0000094C 
00000950 
O00C0960 
00000970 
00000980 
00000090 
00000995 
00001000 
C0001  005 
00001007 
00001009 
00O01015 
00001020 
00001025 
00001030 
ODOO  1040 
00001050 
00001060 
00001070 
00001080 
00001090 
00001100 


133 


C UFDOP  IS  "HE  UPLINK  DOPPLFF  ?HIFT  IN  KZ.  0000111* 

C DUDOP  IS  THE  DOWNLINK  POFPIER  SHIFT  IN  HZ.  00001120 

UPDOP  = DOT*  FP/CL  00001130 

DN  DOP  = DCT*  Fp  1/CL  00001140 

WPTTE  (6  ,20)  J,RI,RGP,P,  AZINUT,UPPOP,  DNPOP,  FS.TS,  K 00001150 

2°  FOpf1  AT  (1  H , T3,F8. 1,Fi7. 3,F14.  3,F15.  3,F1P.  3.F17.  3,F13.3,F13. 3,  00*01160 

1110)  000C1170 

600  CONTINUE  0000118C 

WRITE  (6  ,2P)  00001190 

26  FOPMA”(1H  ,//)  00001200 

-’00  CONTINUE  00001210 

90*  CONTINUE  0*001220 

WP  I'T,p  ( 6 , i*)  00001230 

1*  FOPN A T ( 1 K ,/////)  000*1240 

L = NUM  S- 1 00001250 

IF (L . FQ . 0)  GO  TO  350  00001260 

C FP  IS  THE  CP*  SSI ItfK  FPrQU  ENCT  IN  HZ.  00001270 

F==39 . »1 0 . «»9  *0001280 

C Ml  IS  T HE  FIPST  SA’"-LLI1’E.  00001290 

0 F2  IS  THE  SECOND  S A.TPLLTT p.  00001300 

00  360  M 1 = 1 ,L  00001310 

TAU  = CON*A  (M  1)  **1.  6 00001320 

N=  Ml ♦ 1 000*1330 

00  26*  R2=N,NUKS  0000134C 

WPIT"(6,4*)  00001350 

U * FO?MAT(1H  ,2X,*HPS  • ,12X,  ' PANGF*  , 1 1X,*L12*  ,12X,'C12' ,11X,  00001360 

2 • CFDOPPLEP'  ,6X, 'CVIS*  ,6X, 'SA.TI  • ,6X,  '5AT2' ,//)  00001370 

TAU2  = CON*A  (M2)  »*  1 . p 00001380 

00  30*  1 = 1, 2P  00001390 

PI =T  00001400 

C NOVIZ  DFTPRMINES  CROSSLINK  VISIBILITY  WHFFE  0 MEANS  NO  VISIBILITY  00001410 

C AND  1 MEANS  VISIBILITY.  00001420 

NOVI  Z=  1 0*001430 

CALL  FLLIP  (pI,  F (M  1)  ,WP  (Ml)  ,TP  (Ml)  ,A  (Ml)  ,TAU,FSP,P)  0C0C1440 

CALL  ELIIP(PI,F(M2)  ,WF(N2)  ,TP(M2)  ,A  (M2)  ,TAU2,P2P,P2)  00001450 

CALL  PRIME  (ESP,  W (Ml ) ,XI  (MU  , P,0.  ,XS,  YS.ZS)  00001460 

CALL  FBTME(F2F,W(M2),XT  (M2)  ,R2,C.,X2,Y2,Z2)  0000147C 

TT1=X2-XS  00001480 

TC  2 = V 2 - Y S 0*001490 

TC3=Z  2-  ZS  00001600 

C F 1 2 IS  tHf  FANGF  BETWEEN  TWO  SATELLITES  IN  NAUTICAL  MILES.  00001510 

C FI  2 AND  T12  A PF  THP  POINTING  ANGLES  BETWFEN  TWO  SATELLITES  IN  DEGREES. 000C1 520 
? 12=  (TCI**  2*TC2»*  2*”C3»*2)**.5  *000153  0 

rl2  = PTr*A~AN  (TC2/TC1)  00001540 

IF  (TC1.LT.  0.)  FI  2=  FI  2*  P*  FT  D 00001550 

VIZ=(-XS*TC1-YS*TC2-ZS*TC3)  /(F12*R)  00001560 

I F (VIZ. GT. . *9990000) viZs. 0q9990O9  00001570 

DSL=APCOS  (VIZ)  0000158C 

DFLM=APSIN  (PF/P)  000*1590 

IF  ( DEL  . LE . D ELM)  NOVIZ  = 0 00001600 

T12=APC OS  (TC3/P12)  *PTP  00001610 

CALL  DOPE  (P  ,*  (Ml)  ,A  (Ml)  ,T  AU,PSP,  WP(M1)  , W (Ml)  , XI  (HI)  , XV,YV,ZV)  00001620 
CALL  DOPE  (R 2,  E (M2)  , A (M 2)  , TA  U2 ,F2P  , WP  ("2 ) , W ( M2)  , XI  ( M2)  ,S 2 ,02  ,V 2)  000016  30 
DO?=  ( S2-XV)  *TC1*(02-TV)  *?C2*(V2-ZN)  *TC3  00001640 

DOT=DOT/F1 2 00001650 

C DOP 1 2 IS  THE  CROSSLINK  DOPPLE?  SHIFT  IN  HZ.  00001660 

DOP12=-DOT*FP/CL  00001670 

WRITS  (6,5*)  RT,P12,p12,T12,DOP12,NOVIZ,H1,H2  00001680 

6*  FORM  AT ( 1 H , F' . 1 , F 17 . 3 , F 1 4 . 3 ,F  1 5 . 3 ,F 20 . 3 , 19 ,1 10 , 1 10)  00001690 

30*  CONTI NUr  00001700 

WPI TE  (6  ,60 ) 00001710 
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60  *OSNAT(1H  ,//) 

250  CONTTNOr 

350  CONTISUF 

END 

SU8P09TT NF  P?IH5(FSP,MS,XIS,XPS,YPS,XS#YS,ZS) 
A11=COS(FSP)*COS(MS)-CCS  (XIS)  *SIN  (MS)  *SIN(*SP) 
a’2=-sin  (fsp)*cos  (ms) -cos  (xis)  *sih(»s)  *cos(Fsp) 

A21  = COS(FSP)  *SIN(H?)*COS(XIS)  *COS(MS)*SIN  (FSP) 

A2  2=-SIN  (FSP)  *SI”  (MS)*COS  (X  IS)  *COS  (MS)* COS  (FSP) 

A3 1 = S IN  (XIS ) *S  IN  (*SP) 

A32=SIN  (XIS)  *COS  (FSP) 

XSsA 1 1*XPS*A12*YPS 

yS=A21*XPS*A22*YPS 

ZSsA?1*XDt:*A32*YPF 

FFTUPN 

END 

SOBFOU’'INB  FLI I P ( T,  =■,  MP,  TP,  A,TA0,?SP,3) 

P=3. 1415926 
Z=2.*P*  (T-TP)  /TAD 
t>2=2.  *P 

2 I*(Z.ST.P2)  Z=Z-P2 

IP(Z.GT.P2)  GO  *0  2 
Pl=Z*F*SIN(F) 

F2  = (Z*F*  (SI«  (FI)  ) - (F*COS  (FI)  ) *E1)  /(  1 . -E«COS  ( E 1)  ) 

0*'’. 

a r3*  (Z*F*(STN(f  2) ) - (F*COS  (F2) ) *F2)  / (I  .-F*COS  (F2)  ) 

0=0* 1 . 

D*=F3-F2 

DS2*nE**2 

F2=F3 

TF(0F2.GT.  . :CCnCOr  ’(GO  TO  a 
7H*>PC05  ( (COS  (F2)  -F)  / ( 1 . - F'COS  ( F2 ) ) ) 

I*  (Z.GT.  P)  TH=2.  *P-TH 
F9P  = M?4-~H 

P=  (A*(  3.-F**2)  ) /( 1.*F*C0S (FSP-MP)  ) 

FF~o»N 

ENO 

SPBRCU'INF  OCof  (P,  ",  A,  T AO,  FSP  , MP  , MS  ,XIS  , Xt>,  YD,  ZD) 

P=3.1415°26 

T= FSP-MP 

PD=(A*E*2.  *P/("‘An*(1.-E**2)  **.5)  ) *SIN  (■") 

*7  = (2  . * P/TA  U)  * ((1.  -~»*2(**(-1. 5)  ( *(1.*F*C0S(T))  **2 
X?C=?  D*  COS  (FSP)  -R*  TO*  SIN  (FCP| 

X PC *X PC /3Pftrt. 

TPC  = PD*SIN  (FSP)  ♦ F*TD*COS  (*SP) 

YPC*YPC/3699. 

CALL  PPTHF(P.,MS,XIS,XPC,YPC,XD,YD,ZO) 

P'TOPN 
FN  n 


00001720 
0ni]01730 
00001740 
00001750 
00001760 
00001770 
00001780 
00001790 
00001800 
0000 1 81 0 
000C182C 
00001830 
00001840 
00001850 
00001860 
00001870 
00001880 
00001890 
00001400 
000019  1C 
00001*>20 
0CCC1 °30 
00001940 
CC001950 
00001960 
00001970 
000G19e0 
00001990 
00002000 
90002010 
0C0C2020 
C 0002C  30 
00002040 
00002C50 
00002C6C 
000C207? 
00002080 
09002090 
00002100 
00002110 
00002120 
00002130 
00002140 
C0002150 
90002160 
00002170 
00002180 
0C0021 90 
00002200 
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APPENDIX  9 


PROGRAM  NOLINKE 

i ; 

i 

NOLINKE  searches  for  any  visible  satellite  at  regular  intervals 
around  the  earth.  If  no  satellite  is  visible,  the  coordinates  of 
that  unfortunate  location  are  printed  out  in  the  foreground.  A good 
10- satellite  system  can  require  thousands  of  elevation  angle  cal- 
culations for  very  little  printout.  Care  should  be  used  with  NOLINKE 
or  large  CPU  time  and  expense  will  result. 

The  orbital  elements  (A  in  nautical  miles)  are  entered  on  lines 
20-70.  Elevation  angle  requirements  are  entered  on  line  210  as  EM, 
in  degrees.  Time  Intervals  of  one  hour  are  examined  (RI  - time  in 
hours  on  line  230). 


i 

w 


DINING  I ON  A (10)  ,'(103  , H(1C)  , HP  ( 1 0)  , TP  (1  0)  , XI  (1 0) 

00000010 

DATA  A /1 4342. ,14342.  ,229'?., 227 67., 22767., 22767., 

00000020 

3 V’369., '7369. , 5736°.  ,51360./ 

00000025 

DATA  5/.  7 25,  .725,0.,  0 . , 0. , C . , 0 . ,0.  , 0.  ,0 ./ 

00000030 

DAT*  H/9. ,270. ,0. ,0. ,0. ,". ,9.,0.,0.  ,0./ 

00000040 

DATA  WP/-9C.  ,-9C.  ,0.,  C.  , 0.  ,0.  ,0.  ,O.,0./ 

00000050 

PA ’’A  Tp/0.  ,-6.  ,-3.,-?.,-15.  ,-21 . , 

0000006C 

1 0.  ,-24.  ,-49.  ,-92./ 

000C0065 

DATA  XI/6  3.  43  5,  6 3.  4 35,0.  ,0.  ,0.,C.,?C.  ,9  0.  ,9  0.  ,97./ 

00000070 

HPT  IF  ( 6 , 7) 

00000100 

•J 

FORNAT  ( 1 H ,5X,'T1*  ,«X,  'FI',  13X,  *NO.  • , 13X,*T*  ,/) 

OOODOIIO 

PTD  = 57. 295H95 

00000120 

DO  11  (1=’,10 

00000130 

H (T ) = H (N)  /PTD 

000001 40 

HP(K)=HP(B)  /PTD 

000C0150 

XT  (1)  =XI  (N)  /PTD 

00000160 

1 1 

CONTEND' 

0C0C01 70 

P=3.  1415926 

00000100 

= 3UUD . 

00000190 

HP  A D= 1 5 . 9 /PTD 

00000200 

rr  = ’5. 

0000021 0 

DO  699  -=i,6 

00000  220 

RT=T- 1 . 

00000230 

DO  5 DO  J = 1 , 13 

0000024C 

PJ=J 

0000025C 

(P.1-1  . ) *1  5.  /PTD 

0000C260 

DO  400  K = 1 , 25 

09000290 

PK  = K 

0 0000  280 

FTM=((RK-1.)*16./?TD)+3T*HPAD 

C 0000290 

F=  FIN- P I*HP  AD 

00000300 

TON  =6 .909*10,  ••  (-fj 

00000  310 

SA?NO=r . 

00000320 

DO  300  t=1,  ID 

00000330 

PL  = L 

00000  340 

TA«=CON»A  (L ) **  1 .6 

00000350 

CALI  FLLIP(PI,F  (L)  ,HP(L)  ,T(L)  ,A(L)  ,TA0,F5P,R) 

00000360 

CALL  PRT(1E(?SP,H(L)  ,»T(L)  ,P,d.,XS,XS,Z5( 

00000370 

T = PT*SIN  (T)  *COS  (FIN) 

00000380 

Y = RF«?IN  (?)  *S  IN  ('IN) 

00000390 

Z=R'*COS (T) 

00000400 

ACC=  APCOS  ( (Y*XS  + Y*YS*Z*?S)/  (R'*P)  ) 

00000410 

RGF=  ((XS-X)  **2*  (YS- Y)  *«2*  (ZS-Z)  **2)  **.5 

00000420 

A°G  = SI  N ( ACC)  *P/RGE 

00000430 

1=  ( A PG.GT..9QU9°99)A?G=. 09990Qq 

00000440 

0=  ( APS  IN  ( ARG3  - P/2.3  *°TD 

00000450 

RT  = SCPT(S*P-RF*RF) 

00000460 

IF  (pGr.  LF.  RT)  D=-D 

00000470 

IF(D.GT.  SB)  SATNO  = SATNO*1. 

00000480 

IF  (SATNO.GT.  . 1)  GO  TO  400 

00000490 

IF  (PL. LT.  9. 9)  GO  TO  30  0 

00000500 

T1  = ’t’*°TD 

00000510 

PI  =F*PTD 

00000520 

HRITF(6,20)T1 ,*1 ,SATNO,RI 

00000530 

20 

'ORB  AT  ( 1H  , 4F 1 2 . 3) 

00000540 

300 

CONTINOF 

00000550 

IF  (PJ  .LS.  1 . ) GO  TO  500 

00000560 

IF ( P J . GT.  12.)  GO  TO  500 

00000570 

UCO 

CORTINOF 

00000580 

c np 

CONTINnF 

00000590 

600 

CONTINOS 

00000600 

'MD 

00000610 
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SOBROUTIN  p PPIN  3(FSP,BSfXIS,XPS,YPS,XS,  YS,ZS) 

A1 1 = COS(FSP)  *COS(NS)  -COS(XIS)  *SIN  (BS)  *SIN(FS?) 

A 1 2 --  SI  N (FS  P)  »COS  (BS)  -COS  (XIS)  »SIN  (BS)  *COS  (FSP) 
A2 1 =COS  (FSP ) * STN  (BS)  *COS  (XIS  ) *COS  (BS)  »S  IN  ( FSP) 
A22=-SIN  (FSP)  *SIN  (BS)  ♦COS(XIS)  *COS  (BS)*COS  (FSP) 
A31  = STN(XIS)  *SIN(FSP) 

A32  = SIN  (XIS)*COS  (FSP) 

XS=A11»XPS*A12*YPS 
YS=f.  21*XPS*A22*YPS 
ZS  = A3i*XPS*A32*YPS 
RFTOPN 
END 

SUBROUTINE  PLL  TP  (T , E,  BP  ,T  P,  A,  T AO  , FS  P,  R) 

P=  3.  1415926 
Z = 2.*P*  (T-TP)  /TAO 
P2=2.*P 

2 IF(7  .GT.P2)  Z = Z-P2 

IF  (2.GT.  P2)  GO  TO  2 
E1  = Z SIN  ( Z) 

E2=  <Z*E*(SIN  (El)  )-  (e*COS(E')  ) *31)  / (1  .-E*COS  (El)) 

Q=n. 

tt  *3  = <Z**»  (SIN  (F2)  )-  (r*COS  (E2))*I2)/(1.-Z*COS(E2)  ) 

Q=Q  + 1. 

DF=E3-F2 
FE2  = OF**2 
r2=f  3 

IE  (DE2.GT..r''00“C?1)GO  TO  4 

TH  = AFCOS  ( (COS  (E2)  - E ) / ( 1 . -F*C  OS  (E2)  ) ) 

ir(?.G?.P)  TH«2.*P-TH 

*SP=NP*TH 

°*  (A*  ( 1.  -S**2)  ) /(  1 . *F*COF  (PSP-BP)  ) 

9E“0PN 

END 


C00C0620 
00000630 
00000640 
00000650 
00000660 
CC000670 
00000680 
00000690 
00000700 
00000710 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 
00000780 
00000790 
CC000800 
00000810 
00000820 
0'1000830 
00000°40 
Of 00085C 
00000852 
00000860 
00000870 
0 DO 00 880 
0000089C 
00000900 
00000910 
00000«20 
00000930 
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